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SECTION  I 
INTRODUCTION 

In  this  report  user's  information  for  the  ADDRESS  program 
(Automated  Design  of  Damage  Resistant  Structures)  is  presented. 

The  first  volume  of  this  report  discusses  the  theory  behind  the 
ADDRESS  program  and  presents  representative  applications.  In 
this  introductory  section  an  overview  of  the  user's  instruction  rttehual 
is  given.  Also  included  is  a  discussion  of  the  supplemental 
data  found  in  the  appendices. 

1.1  OVERVIEW 

Section  2  of  this  report  details  the  data  input  requirements. 
Particular  attention  should  be  paid  to  the  user  notes  at  the  end 
of  that  section. 

Program  output  is  discussed  in  Section  3.  This  includes 
printed  output  as  well  as  informaton  to  local  files  for  further 
processing.  Formats  for  these  files  are  specified  in  the  code. 

Control  cards  and  program  updates  are  given  in  Section  4. 

Several  distinct  versions  of  ADDRESS  can  be  easily  obtained  to 
handle  vibration  problems  and  models  which  use  beam  finite  elements. 
Segmentation  directives  are  also  given  which  help  to  reduce  the 
core  requirements  for  analysis  problems. 

The  organization  of  the  subroutines  is  detailed  in  Section  5. 
This  information  is  helpful  in  adding  new  features  to  the  program. 

1.2  ORGANIZATION  OF  THE  SUPPLEMENTAL  DATA 

In  the  appendices  other  useful  programing  information  are 

presented.  Appendices  A  and  B  were  generated  using  the  UDRI  cross- 

2 

reference  program  XREF  .  The  first  two  pages  of  Appendix  A  con¬ 
tain  a  brief  explanation  of  the  use  of  XREF. 

The  variable  name  -  subroutine  name  cross-reference  table  of 
Appendix  A  gives  an  alphabetical  listing  of  all  the  program 
variables.  After  the  colon  following  the  variable  name  is  a  listing 
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of  all  the  subroutines  where  this  variable  is  used.  An  asterisk 
before  the  subroutine  name  indicates  that  the  variable  is  defined 
or  redefined  in  that  routine.  Appendix  B  is  a  list  of  the  common 
blocks  and  the  routines  in  which  they  appear. 

The  subroutine  descriptions  of  Appendix  C  provide  the 
following  information: 

* 

•  subroutine  name  {the  descriptions  are  in  alphabetical 
order  by  name) ; 

4 

•  a  short  statement  of  the  purpose  of  the  routine: 

•  the  call  sequence  with  the  arguments  as  they  appear  in 
the  subroutine ; 

•  a  description  of  the  arguments; 

•  a  list  of  the  subroutines  which  call  the  given  subroutine; 

•  a  list  of  the  externals,  including  any  intrinsic 
functions ; 

•  the  files  used  by  this  routine; 

•  special  notes  or  user  information. 

The  updates  in  Appendix  D  show  the  changes  that  are  needed 
to  put  the  mass  and  stiffness  matrices  out  of  core.  The  PUTSK, 

PRNTSK,  and  SK  function  routines  are  listed  with  extensive  user 
comments. 

Appendix  E  is  a  listing  of  the  TRUSS  program  discussed  in  . 

Reference  1  which  compares  various  iterative  reanalysis  techniques 
for  the  calculation  of  stresses  and  deflections  in  a  damaged  truss  *■' 

structure.  The  program  is  written  for  the  CDC  Fortran  V  compiler. 
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SECTION  2 

INPUT  INSTRUCTIONS  FOR  THE  ADDRESS  PROGRAM 

In  this  section  the  data  requirements  of  the  ADDRESS  pro¬ 
gram  are  discussed.  The  input  consists  of  the  following  blocks 
of  data: 

A.  General  Input 

B.  Coordinate  Data 

C.  General  Element  Data 

D.  Materials  Data 

E.  Connectivity  Data 

F.  Fiber  Orientation  Data  (optional) 

G.  Boundary  Conditions 

H.  Loads  Data 

I.  Displacement  Data  (optional) 

J.  Damage  Data  (optional) 

K.  Lumped  Mass  Data  (optional) 

The  following  pages  describe  each  of  these  data  blocks  in  detail. 
Each  data  block  consists  of  one  or  more  card  sets,  and  each  card 
set  consists  of  one  or  more  cards.  The  notes  inform  the  user  as 
to  the  number  of  cards  in  each  set  as  well  as  explanatory 
information , 
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A.  GENERAL  INPUT- READ  IN  SUBROUTINE  INITIAL  (CARD  SET  1)  AND 
INGNRL  {CARD  SETS  2-5) 


CARD  SET 
{FORMAT) 

VARIABLE 

NAME 

DESCRIPTION 

NOTES 

I 

(115) 

NSTR 

Number  of  problems  to  be  solved 

2 

{ 8A10 ) 

TITLE 

Any  alphanumeric  description  of 
the  problem 

- 

3 

(1615) 

MM 

=2  for  a  two-dimensional  problem 
=3  for  a  three-dimensional  problem 

- 

LMTDSP 

=0  if  there  is  no  displacement 
constraint 

=1  if  displacement  constraint  is 
the  same  for  all  nodes 
=2  if  displacement  constraint  is 
not  the  same  for  all  nodes 

LMTSTR 

=0  if  no  stress  constraints  exist 
=1  if  stress  constraints  exist 

- 

NDMGCAS 

Number  of  damage  cases 

- 

NLMPMSS 

Number  of  lumped  masses 

- 

MAXDCCL 

Number  of  cycles  of  iteration 
using  the  recursion  relation 
based  on  displacement  gradients 

“ 

MAXECCL 

Number  of  cycles  of  iteration 
using  the  recursion  relation 
based  on  energy  gradients 

— 

IAREAS 

=0  if  the  initial  design 

variables  of  the  elements  are 
set  to  1.0  inches 
=1  if  you  wish  to  input  the 

initial  design  variables  for 
each  element 

INDMIN 

=0  if  the  minimum  allowable  size 
is  the  same  for  all  elements 
=1  if  the  minimum  sizes  of  the 
elements  are  to  be  input 
directly 
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CARD  SET 

VARIABLE 

(FORMAT) 

NAME 

DESCRIPTION 

NOTES 

KANLYZE 

=0 

if  you  wish  to  use  the  program 
for  structural  optimization 

- 

=1 

if  you  wish  to  use  the  program 
for  structural  analysis  only 
(no  resizing) 

MAXSZE 

=0 

if  no  maximum  size  is  to  be 
specified  for  the  elements 

- 

=1 

if  the  maximum  allowable 
sizes  of  the  elements  are  to 
be  input  directly 

LAYERD 

=0 

if  problem  contians  no  layered 
composite  elements 

- 

=1 

if  problem  contains  layered 
composite  elements 

INDANG 

=0 

if  no  composite  elements  or  0° 
fibers  are  defined  per  elements 
with  respect  to  the  global 
coordinate  system 

=1 

if  the  0°  fibers  are  defined 
per  element  with  respect  to 
the  local  coordinate  system 

-2 

if  0°  fiber  orientation  is 
the  same  for  all  composite 
elements  with  respect  to  the 
global  coordinate  system 

MNLAYR 

=0 

if  same  for  composite  elements 

— 

=1 

if  the  minimum  proportions  of 

0°,  90°  ±45°  layers  will  be 
input  for  each  member 

IPP 

=0 

if  no  postprocessor  file  present 

— 

=1 

if  postprocessor  file  on  DPOST 
created 

4 

(8F10.3) 


AEMNMM 

DINCR 


THKLAM 

SPRDF 


Minimum  allowable  design  value 

A  parameter  to  determine  the 
active  set  of  displacement 
constraints 

Minimum  composite  layer  thickness 
Shear  panel  reduction  factor 
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CARD  SET 
(FORMAT) 

VARIABLE 

NAME 

DESCRIPTION 

NOTES 

5 

(215,  F10.0) 

MODES 

Number  of  eigenvalues  to  compute 

- 

NO  I 

Maximum  number  of  iterations 

- 

TOLVEC 

Eigenvector  tolerance 
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B.  COORDINATE  DATA  -  READ  IN  SUBROUTINE  INXYZ 


CARD  SET 
(FORMAT) 

VARIABLE 

NAME 

DESCRIPTION 

NOTES 

6 

{ 1 5X , 215 ) 

NODES 

Total  number  of  nodes  in  the  finite 
element  model.  Limit  is  90. 

- 

MGEN 

Mesh  generator  flag 
=0  if  you  wish  to  read  in  the 
coordinates 

=1  if  the  coordinates  are  to  be 
read  or  generated  in  the  user 
subroutine  MESG 

7 

(15, Al, 14 
3E10.0) 

N 

ISYS 

Node  Number 

Reference  coordinate  system 
=_:  Cartesian  system  x,y,z, 

=A:  Cylindrical  R,0,Z 

=B:  Spherical  R,6,<j> 

=C:  (user  defined) 

=D:  (user  defined) 

=E:  (user  defined) 

(1) 

NINCR 

Increment  for  node  point  generation 

(2) 

XI 

X-node  coordinate 

- 

X2 

Y-node  coordinate 

- 

X3 

Z-node  coordinate 

(3) 

•4 
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C.  GENERAL  ELEMENT  DATA  -  READ  IN  SUBROUTINE  ELEMIN 


CARD  SET 
(FORMAT) 

VARIABLE 

NAME 

DESCRIPTION 

NOTES 

8 

(415) 

ITYPE 

Always=3  (fixed  variable) 

{4} 

NMT 

Total  number  of  material  property 
sets.  Limit  is  20. 

- 

NELEM 

Total  number  of  elements 

- 

NC0MP 

Number  of  composite  material 
property  sets.  Must  always  be 
less  than  or  equal  to  NMT 

- 
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D. 


MATERIALS  DATA  -  READ  IN  SUBROUTINE  INP03 


CARD  SET 
(FORMAT) 


VARIABLE 

NAME 


DESCRIPTION 


9  Y0UNGM(I)  If  material  is  isotropic  (same 

(5E10.0)  elasticity  in  all  directions) 

then  it  is  the  Young's  modulus 
for  the  I-th  material  (psi) 

If  material  is  composite  then  it  is 
the  elastic  modulus  in  the  0° 
fiber  direction  for  the  I-th 
material 

P0IS0N ( I )  If  the  material  is  isotropic  then 
it  is  Poisson's  ratio  for  the 
I-th  material 

If  the  material  is  composite  then 
it  is  Poisson's  ratio  for  the 
transverse  strain  due  to  stress 
along  the  0°  fiber  direction 
for  the  I-th  composite  material 

RH01(I)  Density  of  the  I-th  material 

ELASM(I)  If  the  material  isotropic  then  =0 

If  material  is  composite  then  it  is 
the  elastic  modulus  transverse  to 
the  0°  fiber  direction  of  the 
I-th  composite  material  (psi) 

SHEARM(I)  If  material  is  isotropic  then  =0 

If  material  is  composite  then  it 
is  the  shear  modulus  for  the 
I-th  composite  material 


10  ALSTRS (1)  Tension  stress  allowable  in  psi 

(5E10.0)  in  the  0°  fiber  direction 

ALSTRS (2)  Compression  stress  allowable  in 
psi  in  the  0°  fiber  direction 

ALSTRS (3)  Tension  stress  allowable  in  psi 
transverse  to  the  0°  fiber 
direction 


NOTES 

(5) 
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CARD  SET 
{FORMAT) 


VARIABLE 

NAME 


DESCRIPTION 


ALSTRS ( 4 )  Compression  stress  allowable  in 
psi  transverse  to  the  0° 
fiber  direction 


ALSTRS (5)  Shear  stress  allowable  in  psi 
transverse  to  the  0°  fiber 
direction 


If  LAYERED  {see  Card  Set  2)  =  0,  go  to  Card  Set  18 
If  INDANG  (see  Card  Set  2)  =  1,  go  to  Card  Set  13 
If  INDANG  {see  Card  Set  2)  =  2,  go  to  Card  Set  14 
If  ITYPE  ji  1,  omit  Card  Sets  11A  and  11B. 


NOTES 


E 


CONNECTIVITY  DATA  -  READ  IN  SUBROUTINE  INC0NN 


CARD  SET  VARIABLE 

(FORMAT)  NAME  DESCRIPTION  NOTES 


11  IELNO 

(915, 

3E10.0)  I TYPE 


IPR 


Element  number 

Element  type 

=1  for  beam  element 

=2  for  axial  bar 

=3  for  triangular  membrane 

=4  for  quadrilateral  membrane 

=5  for  shear  panel 

Material  property  set  for  this 
element 


(6) 


(7) 


LAMF0  Fiber  orientation  parameter 

=0  for  isotropic  element 
=1  for  fiber  orientations 
0°,  90°,  ±45°  in  the 
proportions  .25,  .25, 

. 50  respectively. 

=2  for  fiber  orientations 

0°,  90°  in  the  proportions 
.50,  .50  respectively. 

=3  for  fiber  orientations 

±45°  in  the  proportion  1.00 
=4  for  fiber  orientations 

0°,  ±45°  in  the  proportions 
1/3,  2/3  respectively. 

=5  for  fiber  orientations 

90°,  ±45°  in  the  proportions 
1/3,  2/3  respectively,' 


KGEN  Node  increment  for  element  (8) 

generation 

N0D ( 1)  Local  node  number  1  (9) 

N0D(2)  Local  node  number  2 

N0D ( 3)  Local  node  number  3 

N0D(4)  Local  node  number  4 

THICK  Member  size  (10) 

THICKMN  Minimum  element  thickness  (10) 

THICKMX  Maximum  element  thickness  (10) 
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CARD  SET 
(FORMAT) 

VARIABLE 

NAME 

DESCRIPTION 

NOTES 

11A 

(5X, 

7E10.0) 

VII 

Moment  of  inertia  about  the  y— axis 
at  node  1 

- 

YI2 

Moment  of  inertia  about  the  y— axis 
at  node  2 

- 

211 

Moment  of  inertia  about  the  z-axis 
at  node  1 

- 

ZI2 

Moment  of  inertia  about  the  z-axis 
at  node  2 

- 

XJ1 

Polar  moment  of  inertia  at  node  1 

— 

XJ2 

Polar  moment  of  inertia  at  node  2 

* 

YK1 

y- component  of  the  mass  moment  of 
inertia  at  node  1 

- 

11B 

(5X, 

7E10.0) 

VK2 

y-component  of  the  mass  moment  of 
inertia  at  node  2 

- 

ZK1 

z- component  of  the  mass  moment  of 
inertia  at  node  1 

- 

ZK2 

z-component  of  the  mass  moment  of 
inertia  at  node  2 

- 

YM 

ZM 

Mass  centroid  offsets 

YG 

ZG 

Geometric  centroid  offsets 

12 


F. 


FIBER  ORIENTATION  DATA  -  READ  IN  SUBROUTINE  INLAYR 


CARD  SET 
(FORMAT) 


12 

( 3F10 . 3) 


GO  TO  CARD 


13 

(6F10. 3) 


GO  TO  CARD 


14 

(3F10. 3) 


IF  I AREAS 


15 

(3F10.3) 


16 

(3F10.3) 


VARIABLE 


NAME 

DESCRIPTION 

NOTES 

XANG(I) 

The  angle  in  degrees  that  the  0° 
fibers  of  the  I-th  element 
makes  with  the  global  x-axis 

Ul) 

YANG (I) 

The  angle  in  degrees  that  the  0° 
fibers  of  the  I-th  elelment 
makes  with  the  global  y-axis 

ZANG(I) 

The  angle  in  degrees  that  the  0° 
fibers  of  the  I-th  element 
makes  with  the  global  z-axis 

SET  16 

XANG(I) 

The  angle  in  degrees  that  the  0° 
fibers  of  the  I-th  element 
makes  with  the  local  element 
coordinate  system 

(12) 

SET  16 

XA 

The  angle  the  0°  fibers  make  with 
the  global  x-axis 

(13) 

YA 

The  angle  the  0°  fibers  make  with 
the  global  y-axis 

- 

ZA 

The  angle  the  0°  fibers  make  with 
the  global  z-axis 

- 

0  SKIP  CARD  SETS  15  and  16 


AEX(I)  Proportion  of  fibers  in  the  0° 

direction  for  the  I-th  element 


AEY(I)  Proportion  of  fibers  in  the  90° 

direction  for  the  I-th  element 
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IF  MNLAYR  =  0  SKIP  CARD  SET  17 


17  AEXMIN ( I ) 

{ 3F10 . 3) 

Minimum  proportion  of  0°  layers 
for  the  I-th  element 

AEYMIN (I) 

Minimum  proportion  of  90°  layers 
for  the  I-th  element 

AEXYMIN ( I ) 

Minimum  proportion  of  ±45°  layers 
for  the  I-th  element 

14 
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BOUNDARY  CONDITIONS  -  READ  IN  SUBROUTINE  GETBC 


CARD  SET 

VARIABLE 

(FORMAT) 

NAME 

DESCRIPTION 

NOTES 

18 

(215) 

NBC1 

Number  of  type  1  constraints 

(14) 

NBC  2 

Number  of  type  2  constraints 

(15) 

IF  NBC1  =  0 

GO  TO  CARD 

SET  21 

19 

N1 

Beginning  node  number  to  be 

(16) 

(315) 

constrained 

N2 

Ending  node  number  to  be 

constrained 

INCR 

Node  number  increment 

20 

JD(I) 

Nodal  components  constrained 

(16) - (17) 

(315) 

(type  1) 

IF  NBC 2  =  0 

SKIP  CARD 

SETS  21  and  22 

21 

JD  (I) 

Nodal  components  constrained 

(17) - (18) 

(315) 

(type  2) 

22 

(1015) 

ND  (I) 

Nodes  constrained 

(18) - (19) 

15 
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LOADS  DATA  -  READ  IN  SUBROUTINE  INLOADS 


CARD  SET 
(FORMAT) 


23 

(1615) 


24 

( 3  (F10 . 0  , 
215)  ) 


IF  LMTDSP 


IF  LMTDSP 


VARIABLE 

NAME 


DESCRIPTION 


NOTES 


LOADS 


Number  of  load  conditions 


NJLOADS 


Number  of 
the  Ith 


load  components  in 
loading  condition 


(20) 


TEMP ( I ) 
IM(I) 


JM(  I) 


Magnitude  of  the  Ith  load 

Direction  of  the  load 
=1  if  in  x  direction 
=2  if  in  y  direction 
=3  if  in  z  direction 

Node  number  where  load  is 
applied 


(21) 


0  SKIP  CARD  SETS  25,  26,  and  27 


2  GOT  TO  CARD  SET  26 
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I. 


DISPLACEMENT  CONSTRAINT  DATA  -  READ  IN  SUBROUTINE  INDSPL 


CARD  SET 
(FORMAT) 


VARIABLE 

NAME 


DESCRIPTION 


25 

(3F10.3) 


DEFMAX(I)  Absolute  value  of  the  displacement 
constraint  in  the  Ith  direction 
for  all  nodes 
=1  for  x  direction 
=2  for  y  direction 
=3  for  z  direction 


IF  LMTDSP  =  1  SKIP  CARD  SETS  26  and  27 


26 

(15) 


KH 


Number  of  displacement  constraints 


27 

( 3 (F10 . 0 , 
215)  ) 


TEMP { I ) 
IM  ( I ) 


JM  ( I ) 


Magnitude  of  the  displacement 
constraint 

Direction  in  which  the  constraint 
is  applied 

=1  if  in  x  direction 

=2  if  in  y  direction 

=3  if  in  z  direction 

Number  of  the  node  where  constraint 
is  applied 


IF  NDMGCAS  =  0  SKIP  CARD  SETS  28,  29,  and  30 


NOTES 

(22) 


(23) 
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J.  DAMAGE  DATA  -  READ  IN  SUBROUTINE  INDMGE 


CARD  SET 

VARIABLE 

{FORMAT) 

NAME 

DESCRIPTION 

NOTES 

28 

( 2 (F10 , 0 , 

TOLDMGD 

Damage  displacement  iteration 
tolerance 

- 

15)) 

MAXDDIT 

Maximum  number  of  damage 
displacement  iteration 

— 

TOLDMGF 

Damage  mode  shape  iteration 
tolerance 

- 

MAXDFIT 

Maximum  number  of  damage  mode 
shape  iteration 

- 

29 

(15) 

NDMG(J) 

Number  of  damage  elements  in 
damage  case  J 

(24) 

30 

IDNG(J, I) 

Ith  element  in  damage  case  J 

(15, 

2F10.0) 

DKF 

Stiff  damage  factor 

- 

DMF 

(TR ( J , I ) ) 

Mass  damage  factor  in  the 

Ith  element  and  in  the  Ith 
damage  case 

* 

IF  NLMPMSS  =  0  THEN  STOP 
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K. 


LUMPED  MASS  DATA  -  READ  IN  SUBROUTINE  INPMSS 


CARD  SET 
(FORMAT) 


31 

(15, 

F10.0) 


VARIABLE 

NAME  DESCRIPTION 

LMN0DE ( I )  Node  location  of  Ith  lumped  mass 
WLMPMSS ( I)  Value  of  Ith  lumped  mass  in  lbf. 


NOTES 

(25) 
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NOTES : 


-J 


(1)  Repeat  Card  Set  5  N0DES  times,  then  add  one  blank  card. 


(2)  A  nonzero  value  of  NINCR  causes  nodes  to  be  equally 
spaced  between  the  last  and  current  nodes,  with 
numbering  increment  NICR.  As  an  example,  the  data 


10 

0. 

0. 

0. 

20  2 

10. 

-10. 

0. 

is  equivalent  to 

10 

0. 

0. 

0. 

12 

2. 

-2. 

0. 

14 

4. 

-4. 

0. 

16 

6 . 

-6. 

0. 

18 

8. 

-8. 

0. 

20 

10. 

-10. 

0, 

If  MM  =  2  (i.e. , 

two-dimensional 

problem) 

then  X3  is 

zero 

ITYPE  in  ELEMIN 

is  not  the  same 

as  ITYPE 

in  INP03. 

In 

subroutine  ELEMIN  the  following  variable 
made : 

assignments 

are 

NMAT  =  NMT 
MEMBS  =  NELEM 

IS0TRN  -  NMAT  -  NC0MP  -  Number  of  isotropic 

material  property  sets 

(5)  Repeat  Card  Sets  9  and  10  NMT  times.  The  isotropic 
material  properties  are  ordered  before  the  composite 
properties . 

{6)  Repeat  Card  Set  11  NELEM  times  if  KGEN  =  0  {see  Note  8 
below) .  Add  two  blank  cards  at  the  end. 

(7)  This  is  determined  by  the  ordering  established  in 
Card  Sets  9  and  10. 

(8)  A  nonzero  value  of  KGEN  on  the  second  card  of  a  pair 
causes  elements  between  the  last  and  current  elements 
to  be  automatically  generated.  With  the  exception  of 
node  numbers,  all  elements  generated  are  assigned  the 
same  data  as  the  current  element.  Local  node  numbers 
for  the  previous  element  are  incremented  by  KGEN  to 
generate  each  succeeding  element.  More  than  one 
element  must  be  generated  to  use  this  feature;  node 
numbers  for  the  current  element  need  not  be  given.  As 
an  example,  the  data 
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85 

5 

1 

0 

0 

49 

50 

52 

51 

.0190 

88 

2 

is 

equivalent 

to 

85 

5 

1 

0 

0 

49 

50 

52 

51 

.0190 

86 

5 

1 

0 

0 

51 

52 

54 

53 

.0190 

87 

5 

1 

0 

0 

53 

54 

56 

55 

.0190 

88 

5 

1 

0 

0 

55 

56 

58 

57 

.0190 

(9)  For  beams,  leave  N0D(4)  blank.  Beam  N0D(3)  determines 
principal  axes  of  bending. 

For  bars,  leave  N0D{3)  and  N0D(4)  blank. 

For  triangular  membrane,  leave  N0D(4)  blank. 

For  each  element,  let  N0D(1)  be  the  lowest  node  number 
and  N0D(2)  the  next  lowest. 

For  quadrilaterals  and  shear  panels,  N0D(3)  and  N0D(4) 
are  determined  by  continuing  in  the  direction  defined 
by  N0D (1) ,  N0D ( 2 ) . 

(10)  For  bars,  THICK  is  the  cross-sectional  area.  For  all  others, 

except  beams,  THICK  is  the  element  thickness. 

For  beam  analysis,  set  THICK  =  0;  THICKMN  =  area  at  end  A 
of  beam;  THICKMX  =  area  at  end  B  of  beam. 

(11)  Include  only  if  LAYERD  =  1  and  INDANG  =  0  (see  Card  Set  3) . 

Read  continuously  for  I  =  1  to  I  =  NELEM.  If  an  element 
is  not  composite,  set  XAND(I)  =  YANG  £  I )  =  ZANG(I)  =  0.0 
for  that  element. 

(12)  Include  only  if  LAYERD  =  1  and  INDANG  =  1  (see  Card  Set  3) . 

Read  continuously  for  I  =  1  to  I  =  NELEM.  If  an 
element  is  not  composite  set  XANG(I)  =  0.0  for  that 
element. 

(13)  Include  only  if  LAYERD  =  1  and  INDANG  =  2  (see  Card  Set  3) . 

(14)  A  type  1  constraint  is  used  to  fix  selected  degrees  of 

freedom  over  a  specified  range  of  nodes  with  specified 
increment  in  node  numbers. 

(15)  A  type  2  constraint  is  used  to  fix  selected  degrees  of 

freedom  over  a  list  of  nodes. 

(16)  Type  1  constraint  data  to  be  entered  NBC1  times.  Omit 

if  NBC1  =  0. 
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(17)  Values  of  1,  2,  and  3  can  be  used.  To  fix  u  at  a  node 

(x  component  of  deflection)  include  a  1  value.  To  fix 
v  at  a  node  (y  component  of  deflection)  include  a  2 
value.  To  fix  w  at  a  node  (z  component  of  deflection) 
include  a  3  value.  For  example,  make  Card  Set  13 

12  3 

to  fix  all  degrees  of  freedom.  To  fix  only  u  use 

10  0 


(18) 

Type  2 
if  NB 

(19) 

Each  of 
speci 

(20) 

Repeat 

(21) 

Repeat 

cards 

(22) 

Repeat  ] 

(23) 

Repeat 

(24) 

Repeat  ] 

(25) 

Repeat  1 
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SECTION  3 
PROGRAM  OUTPUT 

Output  of  the  ADDRESS  program  is  very  similar  in  nature  to  the 
OPSTATJ  program.  Formats  have  been  changed  to  some  extent  and 
additional  information  on  damage  cases  is  printed.  Output  for  a 
typical  static  optimization  run  is  as  follows: 

1.  Input  Data  Echo.  Fifty  lines  of  data  are  printed  per 
page  with  header  and  trailer  numbers  to  help  check  the  column 
locations  of  the  data.  The  title  card  information  is  printed  on 
all  pages  of  the  output,  together  with  a  page  number. 

2.  General  Data.  Card  Set  3  information  is  presented  with 
explanations  of  the  numerical  values. 

3.  Nodal  Coordinates.  Data  is  always  printed  in  Cartesian 
coordinates  even  though  another  system  may  be  used  for  input. 

4.  Materials  Data.  Properties  from  Card  Set  9  and  the 
allowables  from  Card  Set  10  are  printed  for  each  material. 

5.  Element  Connectivity.  Element  number  (ELEM) ,  type  {ITYPE) , 
material  code  (MATL) ,  fiber  orientation  parameter  (LAM) ,  local 

nude  number,  and  member  sizes  are  printed,  based  on  the  Card  Set  11 
data. 

6.  Composite  Element  Data.  If  the  model  contains  composite 
elements  the  fiber  orientation  data  from  Card  Sets  12,  15,  16,  and 
17  are  printed. 

7.  Boundary  Conditions .  A  summary  of  the  type  1  and  type  2 
boundary  conditions  (see  Notes  14  and  15  of  Part  G  of  Section  2) 
is  printed.  For  the  current  library  of  elements,  components  4 
through  10  will  always  be  printed  as  zeroes. 

8.  Load  Summary.  For  each  load  case  the  net  loads  (FX,  FY, 

FZ)  in  the  coordinate  directions  are  printed,  together  with  their 
respective  moments  {MX,  MY,  MZ) . 
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9-  Damage  Data.  Data  from  Card  Sets  28,  29,  and  30  are 
printed  for  each  damage  case. 

10.  Population  Information. 3  Output  from  Subroutine  P0P 
concerning  the  distribution  of  elements  in  the  stiffness  matrix. 
This  information  is  generated  before  the  stiffness  matrix  of  the 
structure  is  assembled. 

(a)  Gross  Population  =  total  number  of  elements  in  the 
upper  triangle  of  the  matrix. 

Apparent  Population  =  actual  number  of  elements 
considered  as  non- zero  by  a  given  solution  scheme. 

Thus  the  apparent  population  represents  the  number 
of  storage  locations  required  for  the  stiffness 
matrix . 

(b)  Starting  Row  Numbers  for  each  column  -  the  number  of 
the  row  where  the  first  non— zero  element  occurs  in 
each  column. 

(c)  Numbers  of  Diagonal  Elements  in  Single  Array  Stiffness 
Matrix.  For  each  Column  I  the  actual  number  of 
elements,  ID  (I),  in  the  upper  triangular  matrix  up 

to  and  including  that  column,  i.e., 

ID  (I)  «=  -  l  b. 

j=l  3 

where  b..  is  the  row  number  given  for  Column  I  in  (b) . 
Thus  for  the  last  column,  ILAST, 

ID (I LAST)  -  Apparent  Population 

11-  Relative  Design  Data.  After  the  static  analysis  is 

performed,  the  relative  sizes  of  the  elements  are  printed. 
The  maximum  relative  size  is  1.0  and  all  other  sizes  are 
given  relative  to  this  largest  size.  Next,  a  summary  of 
the  time  for  the  analysis  is  printed.  This  is  broken 
down  into  time  for  the  matrix  assembly,  decomposition,  and 
forward  and  backward  substitution. 
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12 •  Damage  Reanalysis  Error  Summary.  If  damage  cases  are 
included,  then  the  analysis  iteration  error  summary  is  printed 
for  each  case.  The  norm  of  the  vector  error  is  printed  with  the 
load  case  index  in  parentheses.  Error  information  on  the 
extrapolation  process  is  also  printed.  Total  time  required  to 
analyze  each  damage  case  is  indicated.  This  should  be  less  than 
the  time  for  the  analysis  of  Item  11  above  for  the  reanalysis 
to  be  cost-effective, 

13.  Initial  Scale  Factor  (BASEA) .  This  is  defined  in  sub¬ 
routine  SCALE  to  be  .01  times  the  square  root  of  the  quotient  of 
STRMAX  and  ENGCAP.  The  variable  STRMAX  is  the  maximum  over  all 
load  vectors  of  the  inner  product  of  the  load  and  displacement 
vectors.  In  STIFFK ,  ENGCAP  is  defined  as  the  sum  over  all  elements 
of  AE  times  ELENGTH  where  AE  is  the  member  size  and  ELENGTH  is 
length  (for  bars)  or  area  (for  membranes  and  shear  panels) .  The 
initial  BASEA  is  used  in  PREPAR  to  scale  the  stress  allowables. 

These  scaled  allowables,  ALS,  are  used  in  STRC0N  to  compute  the 
element  stress  ratios,  ESRTI0. 

14.  Scaling  Factors.  If  a  stress  ratio  for  an  element  in 
some  load/damage  case  combination  exceeds  1.0,  its  maximum  value 
is  printed.  The  listing  of  corresponding  element  numbers  that  is 
printed  beside  these  exceedance  values  tells  the  user  which  elements 
are  overstressed  for  the  design  whose  relative  member  sizes  are 
scaled  by  the  initial  scaling  factor  of  Item  13.  The  product  of 
the  initial  scaling  factor  and  the  critical  element  scaling  factors 
gives  the  new  BASEA  which  is  printed  next  with  the  header,  "Scaling 
factor  to  satisfy  stress  constr. " 

15.  Summary  of  Current  Cycle.  Current  cycle  number  of  this 
resizing  is  printed,  together  with  a  weights  summary . 

15.  Resizing  Continues.  Items  11  through  15  are  repeated 
until  one  of  the  following  happens: 

(a)  All  energy  and  displacement  cycles  have  been  completed 
as  specified  in  Card  Set  3; 
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(b)  The  weight  goes  up  in  an  energy  cycle; 

(c)  The  weight  doubles  in  a  displacement  cycle. 

17*  Final  Design  Stress  Summary.  Stresses,  member  sizes, 
element  numbers,  and  ESRATI0's  are  printed.  The  column  of  ESRATl0's 
are  the  last  set  of  ESRATI0's  computed  (see  Item  13).  They  are 
ordered  so  that  the  undamaged  load  cases  come  first,  followed  by 
the  damage  case  information  by  load  case.  Example:  2  load  cases, 

2  damages, 


ELMT 

element  no.  5 


LC  =  load  case, 


ESRATI0 


.  91162E-01 

LC1, 

no  damage 

. 74404E-01 

LC2 , 

no  damage 

. 12544E+00 

LC1 , 

DC1 

. 10595E+00 

LC2  , 

DC1 

. 91743E-01 

LC1 , 

DC2 

. 75008E-01 

LC2, 

DC2 

DC  =  damage 

case . 

1 ® •  Final  Design  Deflection  Summary .  De flections  are  printed 
by  node  in  the  three  coordinate  directions.  Applied  loads  are 
also  printed  in  the  same  format  as  the  ESRATI0's  above  in 
Item  17. 


Additional  output  is  written  to  local  files  8  and  99.  The 
subroutines  beginning  with  "PP"  write  information  to  one  or  both 
of  these  files  for  post-processing.  See  the  subroutine  descrip¬ 
tions  in  Appendix  C  for  further  information. 
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SECTION  4 


CONTROL  CARDS  AND  PROGRAM  UPDATES 

Table  1  lists  the  control  cards  needed  to  run  the  version 
of  ADDRESS  which  is  dimensioned  for  453  elements  {design  variables) 
and  339  degrees  of  freedom.  The  central  memory  requirement  of 
255K  octal  can  be  reduced  to  155K  octal  through  the  use  of  the 
control  cards  and  segmentation  shown  in  Table  2.  When  segmenta¬ 
tion  is  used,  the  program  capacity  is  doubled;  however,  the  1/0 
requirement  also  doubles. 


TABLE  1 

CONTROL  CARD  EXAMPLE  WITHOUT  SEGMENTATION 


UDNT  -  T3000 .  1 02000  >  CM 2 5 50 00 .  D790486 ,  T AYL0R »  30  i 

COMMENT.  INTERCOM  BATCH  JOB  ****  Nu  DECK******** 

AT TACH . LOO ,  DAM0PTBI NARY JENSEN >  MR= 1 . 

ATTACH, TAPES, NEWDATAFORMETALWI NG» CY-3, MR=1 . 

MAP ( PART ) 

LIMIT, 3 100. 

I.  00  v  PL88 10000 ) 

REWIND, DEBUG, 
i  OP Y. DEBUG. 

£  <  I  T  t  S  i 

l-'b'W  J ND-  DEBUG. 

r OPY. DEBUG. 


93672 
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TABLE  2 

CONTROL  CARD  EXAMPLE  WITH  SEGMENTATION 


UD5T ,  T .1  SO ,  1 0800 ,  CM  1 55000 .  D790486 , TAYLOR,  229-30 1 S .  93672 

COMMENT.  INTERCOM  BATCH  JOB  ****  NO  DECK******** 

ATTACH, OLDPL , RFTDAMOPTPL. 

ATTACH,  C,  RFT D AMO P T MODE Z E R OCHAN G E S . 

COMMENT.  THIS  UPDATE  CREATES  THE  M0DES=0  DAMOPT 
UPDATE, F, I*C, L=A 1 . 

FTN, 1 =CQMP I LE , R=3 . 

COMMENT.  ATTACH, INPUT  DATA  FILE  AS  TAPES 
ATTACH, TAPES, W I NGDAMAGECASEECS, CY=2 . 

MAP, ON. 

LDSET, PRESET* ZERO. 

SEGLOAD . 

LOO. 

TREE  DAMOPT- (DOPROB, INITIAL-ECHO) 

LEVEL 

TREE  AVERAGE 
TREE  SAVE 
TREE  EMODE 
TREE  GETBEST 

TREE  INI  TP- (POP, I NPT - <  TREE 1 »  GETBC, INDMGE, INDSPL, INGNRL-PPTTL , 

,  INLAYR,  INLOADS.  INLPMSS.  I  NX  YZ  -  <  CTYPE ,  ME3HG ,  FPNODE ) ')  ) 

T REE  1  TREE  ELEM I N- 1 NP03- 1 NGONN- ( NODCHCK , PPELEM ) 

1  REE  PR  1  NT-  ( LA YFR-LA YCALC ,  PRNTDR-PPDISP ,  F'RNTEL ,  PftNTMDS  ) 

TREE  DMODE-UNTF RC-UN I TEG 
LEVEL 

TREE  CURRENT- ( WEIGHTS, SI VI B2-<  DECOUP, ERROR, ORTHOG, PREMULT, RANDOM) 
, , SCALE -S TROON , DMGFREQ- ( FRQITER-DMGORTH-PRELT ) , 

, ANALYZ- (STIFFK-LLT-PRINTK, B0UND2) ) 

CU R R E  N T  1 NCLU D E  B ACKS U B , F 0 R SU B 
LEVEL 

TREE  LMSIZE 
I  REE  PPMODE-GENMSS 
\ REE  RESTOR 
I  REE  PPHPR 
i REE  REDUCE 

I  REE  QDRLTL- ( PLST I F-CRAMER , CONDNS-CHANGE , SUM ) 
l  .KVEL 

I  REE  ELFORC 

1  REE  DAMAGE- ( DMG IT ER-DX , ERR , GETDK-ASEMDK ) 

LEVEL 

TREE  TRNSFM 

TREE  COuRD 

TREE  ELSTIC 

(REE  ELSTIF-LMPROD 

TREE  LMPMAS 

T RE E  PLMASS 

I  REE  PR E PAR 

TRFE  TRECON 

T  U 1 1  'I  RCi  Li  ST  R-T  R Q  DO  U T-QLSTR  S-S  T  R  E  SS 
i  urlMON 

1 1 1” i  i  l-  ’ 'L  OBAL  DUMMY-SAVE 
i.iIRrENl  GLOBAL  GM-SAVE ,  SK-SAVE 
END  DAMOPT 
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TABLE  2  (continued) 


*  1 D  DELETE! 

RFTDAMOPT. 52 

COMMON  /DM/  GM( 1 ) 
*D  RFTDAMOPT. 236 

COMMON  /OM/  OM( 1 ) 
*D  RFTDAMOPT. 2946 

COMMON  /CM/  GM( 1 ) 


*D  RFTDAMOPT. 3707 

COMMON  /GM/  GM( 
*D  RFTDAMOPT. 555 

COMMON  /EIGVEC/ 
*D  RFTDAMOPT. 646 

COMMON  /EIGVEC/ 
*D  RFTDAMOPT. 713 

COMMON  /EIGVEC/ 
*D  RFTDAMOPT. 799 

COMMON  /EIGVEC/ 
*0  RFTDAMOPT.  1247 

COMMON  /EIGVEC/ 
*D  RFTDAMOPT. 1305 

COMMON  /EIGVEC/ 
*D  RFTDAMOPT. 1340 

COMMON  /EIGVEC/ 
#D  RFTDAMOPT. 1387 
5)40  FORMAT ( 3F10. 3 ) 


) 

U  ( 350 ,  6 ) ,  V  ( 350 »  6 ) 
U  <  350 ,  6 )  »  V  ( 350 ,  6 ) 
U ( 350 1 6 )  »  V  ( 350, 6 ) 
LK350.6) ,  V(350,6) 
U<  350,  6 ) »  V  ( 350 1 6 ) 
11(350.6) ,  VC  350, 6) 
IK 350, 6) , VC 350,  6) 


.M(Mi  1  )  »  XL  dO»  3) 
«U<li 1* 1 ), XL< 10,3) 
»M(1. 1. 1  >,  XLU0.3) 
,W<1, 1, 1 ), XLU0.3) 
.M< 1. 1 . 1 ) i XL( 10.3) 
.MU.  1*  1  )  .  XLtlO.3) 
,W(1. 1. 1) . XL< 10.3) 


*D  RFTDAMOPT. 2371 

COMMON  /El GVEC/  U < 350 , 6 ) , V ( 350, 6 > .  U  U  .  1 .  1  > .  XL UO.  3 > 
*D  RFTDAMOPT. 3262 

COMMON  /EIGVEC/  U( 350. 6 ) . V( 350, 6 > , W ( 1 , 1 , 1 ) , XL( 10, 3 ) 
*D  RFTDAMOPT. 3633 

COMMON  /EIGVEC/  LK 350 , 6 > , V < 350, 6 ) , W ( 1 , 1, 1 ) , XL( 10,3) 


Y(  1.  1) 

Y  <  1 ,  1  ) 

Y  (  1 ,  1  ) 
YU,  1) 
Yd,  1  > 
Y<  1 ,  1  > 
Y(l,  1  ) 

Y  (  1 ,  1  ) 
Y(l,  1  ) 
Yfld  ) 
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SECTION  5 

PROGRAM  ORGANIZATION 

The  ADDRESS  program  contains  101  routines,  counting  the 
main  program.  Figures  1  through  6  show  their  interrelationship. 

In  this  section  the  overall  organization  of  this  calling  sequence 
will  be  outlined.  For  details  of  the  individual  subroutines, 
see  Appendix  C. 

The  main  program  first  calls  INITIAL  which  sets  certain 
variables  and  reads  NSTR  which  is  the  number  of  problems  to  be 
solved  in  this  run.  The  input  data  is  echoed  and  mass  storage  is 
set  up  for  the  global  mass  and  stiffness  matrices.  All  of  the 
analysis  and  optimization  functions  are  controlled  by  D0PR0B, 
which  is  called  within  a  loop  whose  index  ranges  from  1  to  NSTR. 

The  D0PR0B  routine  shown  in  Figure  1  contains  the  program 
logic  for  performing  all  the  aspects  of  the  optimization.  Its 
functions  are  as  follows: 

(1)  Read  the  input  data  (call  to  INITP) 

(2)  Analyze  the  current  design  (call  to  CURRENT) 

(3)  Perform  resizing  for  strength  requirements  (call  to 
EM0DE) 

(4)  Perform  resizing  to  meet  displacement  contraints  (call  to 
DM0DE) 

(5)  Perform  bookkeeping  functions  such  as  getting  the 

best  design,  averaging  with  a  previous  design,  and  saving 
the  results  (calls  to  GETBEST,  AVERAGE,  and  SAVE) 

(6)  Print  the  results  (call  to  PRINT) 

During  the  input  phase  the  routines  of  Figure  2  are  called. 

In  addition  to  input,  a  preliminary  estimate  of  the  storage  space 
required  for  the  mass  and  stiffness  matrices  is  made  by  subroutine 
P0P.  The  input  routines  read  the  following  types  of  data: 
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(1)  General  data  such  as  analysis  and  optimization  options 
to  be  performed  (call  to  INGNRL) ; 

(2)  Coordinate  data  for  the  nodes  in  the  finite  element  model 
(call  to  INXYZ) ; 

(3)  Element  data  including  material  properties,  stress  allow¬ 
ables,  and  connectivity  (call  to  ELEMIN) ; 

(4)  Layer  data  for  any  composite  elements  (call  to  INLAYR)  ; 

(5)  Boundary  condition  information  (call  to  GETBC) ; 

(6)  Lumped  mass  data  (call  to  INLPMSS) ; 

(7)  Applied  load  conditions  (call  to  INL0ADS) ; 

(8)  Information  defining  the  damage  conditions  (call  to 
INDMGE)  and  the  displacement  constraints  (call  to  INDSPL) . 

All  of  the  input  data  are  also  printed  by  these  routines,  and  the 
model  data  is  also  put  on  the  local  file  NFIL=99  for  optional 
post-processing  of  the  model. 

Figure  3  shows  the  relationship  of  the  various  routines  which 
perform  the  analysis  of  the  current  design.  The  functions  performed 
by  CURRENT  are  as  follows: 

(1)  Set  up  the  global  mass  and  stiffness  matrices  and  perform 
a  static  analysis  (call  to  ANALYZ ) ; 

(2)  Perform  a  vibration  analysis  (call  to  SIVIB2) ; 

(3)  For  the  damage  conditions,  compute  the  static  (call  to 
DAMAGE)  and  vibration  (call  to  DMGFREQ)  response 

by  reanalysis  methods; 

(4)  Scale  the  design  to  meet  stress  requirements  (call  to 
SCALE)  and  compute  the  weight  (call  to  WEIGHT) . 

The  routines  shown  in  Figure  4  are  called  when  resizing  for 
displacement  constraints  takes  place.  Virtual  displacements  due 
to  the  dummy  loads  are  determined  by  calls  to  REDUCE ,  F0RSUB ,  and 
BACKSUB.  Calls  to  REST0R  and  UNTFRC  determine  response  to  the 
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loads.  Damage  cases  are  handled  by  the  call  to  DAMAGE.  After 
resizing ,  LMSIZE  is  called  to  determine  the  thickness  distribu¬ 
tions  for  any  composite  elements  that  have  been  resized. 

Figures  5  and  6  show  the  relationship  of  various  utility 
routines.  The  print  routines  of  Figure  5  involve  calls  to 
analysis  routines  to  obtain  results  for  the  final  design.  The 
utility  routines  of  Figure  6  set  up  the  element  mass  and  stiffnesses 
for  both  damaged  and  undamged  elements . 
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MAIN 


D0PR0B 


INITIAL 


INITP  -® 
CURRENT  -® 
SAVE 

EM0DE  -  LMSIZE 
GETBEST 
DM0DE  -  (c) 
AVERAGE 

-  PRINT  -(D) 

-  ECH0 

_  0PENSK 


Note:  See  Figures  2,  3,  4,  and  5  for  details  of 

A,  B,  C,  and  D,  respectively. 


Figure  1.  ADDRESS  Program  Routines:  MAIN,  D0PR0B,  and 
INITIAL  Call  Sequence. 
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-  INPT 

©- 

-  POP 


Figure  2. 


r  ingnrl 


-  INXYZ 


—  PPTTL 
-  MESHG 
CTYPE 
PPHDR 
_ PPNODE 


-  ELEMIN 

-  INLAYR 

-  GETBC 

-  INL0ADS 

-  INDSPL 

-  INDMGE 
L-  INLPMSS 


XNP03 

PPHDR 


PPHDR 
INCONN ■ 


NODCHK 

PPELEM 


Input  Routines:  INITP  Call  Sequence. 
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PUTSK 


-DAMAGE 


-  SCALE 


0- 


ANALYZ 


-WEIGHTS 


“SIVIB2- 


L-dmgfreq 


RANDOM 

0RTH0G 

BACKSUB 

PREMULT 

F0RSUB 

DEC0UP 


ERR0R 

-0 


[  STIFFK 

B0UND2 

- PUTSK 

REDUCE 

rPUTSK 

LLT 

1  DT>’MrT|C 

F0RSUB 

-  BACKSUB 

- GETDK  - 

rv^ 

Lasemdk 

-DX - J 

L 

DMGITER  - 

F0RSUB 

-ERR 

- BACKSUB 

'  REST0R 

-C00RD 

ELF0RC 

•  STRC0N 

PREPAR 

- TRQDSTR 

ASEMBL 

■PUTSK 


STRESS 

-© 


NOTE:  For  details  of  E ,  F,  and  H,  see  Figure  6. 


Figure  3.  Analysis  Routines:  CURRENT  Call  Sequence 
(Overall) . 


CRAMER 
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-  REDUCE 

-  F0RSUB 

-  BACKSUB 

-  DAMAGE  ■ 


GETDK 


DMGITER 

-ERR 


KD 

1-ASEMBK 

-  DX - 

F0RSUB 
- BACKSUB 


REST0R 


-  UNTFRC 


- C00RD 
ELF0RC 
PREPAR 

Q 

PLSTIF 

-UNITEG 


CRAMER 


LMSIZE 


NOTE:  For  details  of  E  and  F,  see  Figure  6. 


Figure  4.  Displacement  Resizing  Routines:  DM0DE 
Call  Sequence. 


REST0R 

REDUCE 
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LAYPR 


PRNTEL 


PRNTDR 


-  PRNTMDS 


PPHDR 


-[ 


LAY CALC 
CURRENT 


“ C00RD 


■BMXSTRS 


ELFORC 

PREPAR 


TRQDSTR 
L-  TRQD0UT 


rQLSTRS 

STRESS 


KE) 


PPHDR 

PPDISP 


BMGE0M 


CRAMER 


"  REST0R 

GENMSS 
-  PPM0DE 


■PREMULT 
■  GETDK  — 

-DX  - - 


f“© 

L  ASEMDK 

r  REST0R 
L  REDUCE 


NOTE:  For  details  of  B,  see  Figure  3,  and  for 

E  and  F,  see  Figure  6. 


Figure  5.  Print  Routines:  PRINT  Call  Sequence. 
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-C00R 

-PREPAR 

— © 

--PLSTIF 
- PLMASS 

-  TRNSFM 

-  LMPMAS 
- ELSTIF 


© 

-CRAMER 


LMPR0D 


" ELSTIC 
-  TREC0N 


QDRLTL 


"PLSTIF 

SUM 

PLMASS 

-C0NDNS 


CRAMER 


CHANGE 


BEAMKM 


[ 


BMSYMM 
BMX  - 


BMGE0M 


REST0R 

REDUCE 

REST0R 

REDUCE 


PRELT 

DX' - — £ 


RESTJ0R 

REDUCE 


Figure  6.  Analysis  Routine  Segments  E,  F,  G,  and  H. 
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APPENDIX  A 


VARIABLE  NAME  /  SUBROUTINE  NAME 
CROSS-REFERENCE  TABLE 
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UDRI  -  Cross  reference  program;  XREF 


This  program  was  written  using  FORTRAN  V  on  the  ASD  CDC  Cyber 
750.  It  is  designed  to  read  a  compiler  listing  of  a  FORTRAN  IV 
program  and  generate  a  complete  list  of  variable  and  labeled  common 
names.  It  also  provides  a  set  of  tables  which  identify  which 
variables  are  used  and  defined  in  which  routines  as  well  as  labeled 
commons . 

Program  access: 

ATTACH, XREF, XREFPLUS, ID=KL00S,SN=ASDAD  (source  code) 

ATTACH, XX, XXPLUS, ID=KL00S,SN=ASDAD  (object  code) 

Compilation  of  source  code: 

FTN5 , I=XREF , B=XX, L0=O , ANSI=0 , 0PT=2 

FTN5,I=XREF,B=O,L0=M/A/R,L=LIST,ANSI=O,0PT=2  (for  listing) 

Preparations  for  execution: 

TAPEl  must  contain  the  program  listing  which  XREF  is  to 
read.  XREF  considers  TAPEl  to  be  multiple  routines;  all 
of  the  same  executable  program.  This  does  not  have  to  be 
a  complete  program  nor  arranged  in  any  particular  order. 

To  develop  a  TAPEl: 

FTN, I=lfn, L- TAPEl, R=3, . . . 

TAPEl  can  be  editted  prior  to  execution  of  XREF. 

Execution : 

XREF  can  be  run  at  the  terminal  or  in  batch  in  the  same 
manner  as  any  other  LGO  file.  (The  presence  of  TAPEl  is 
the  only  requirement.)  As  XREF  reads  the  symbolic  reference 
map,  it  will  display  various  error  messages  related  to  the 
routine  it  is  currently  processing.  This  is  useful  in 
finding  FTN  error  messages  in  large  listings.  For  very 
large  listings,  the  OUTPUT  file  can  be  disconnected  from 
the  terminal  and  PAGEd  after  termination. 
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Results : 


XREF  will  display  a  summary  of  the  number  of  routines  in 
TAPEl ,  the  number  of  variable  names  in  the  cross-reference 
table,  and  the  number  of  labeled  commons  found.  It  will 
then  list  all  of  the  routine,  variable,  and  common  names 
in  alphabetical  order  (listing  is  in  132  columns).  Then 
XREF  will  report  the  number  of  'garbage'  variables  and  list 
their  names.  A  'garbage'  variable  name  is  one  which  was 
found  to  be  unused  in  one  of  the  routines.  Some  of  these 
may  occur  in  the  variable  cross-reference  table,  others  may 
not.  XREF  stores  the  cross  reference  tables  on  TAPE 2 .  This 
file  is  rewound  at  the  end  of  the  program  and  can  be  printed 
by:  COPY, TAPE2, OUTPUT.  TAPE 2  is  constructed  in  two  parts; 

variables  and  commons.  The  names  are  listed  in  alphabetical 
order  (by  internal  CDC  display  code)  with  the  routines  they 
are  found  in  listed  beside  them.  For  variable  names,  an 
asterisk  will  appear  with  the  routines  that  were  found  to 
define  the  variable  in  the  R=3  map.  The  tables  are  designed 
to  be  reproduced  onto  a  8  1/2  X  11  inch  paper  with  regular 
typing  margins  (and  no  reduction) . 
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A  : 

AA  : 

AAE  : 

AAEX 

AAEY 

AAX 

AAXY 

AAY 

AB 

ADR 

AE 

AEMAX  : 
AEMNM  : 
AEMNMM  : 
AEX  : 

AEXMIN  : 
AEXYMIN: 
AEY  : 


ANALYZ 

SIVIB2 

* COORD 
TRECON 

AVERAGE 

F’RNTEL 

AVERAGE 

AVERAGE 

^STRESS 

8-STRESS 

^STRESS 

* COORD 

-m-SCALE 

^AVERAGE 
8GETBEST 
■»  LM  PROD 
SAVE 
TRQDSTR 

EMODE 

AVERAGE 

INCONN 

^AVERAGE 

*LMSIZE 

8-1 NLAYR 

# I NLAYR 

^AVERAGE 


*CRAMER 

STRESS 


ELFGRC 

TRNSFM 


GET BEST 
*SAVE 


GETBEST 


GETBEST 


8-UNI  TEG 

CURRENT 
GETDK 
*LMSIZE 
SCALE 
UNI  TEG 

8  INCONN 


DMODE 


8-INGNRL 

*GETBEST 

PREPAR 

LMSIZE 


LMSIZE 


8GETBEST 


DAMAGE 


ELSTIF 


LAYCALC 
TROD GUT 

8SAVE 

*SAVE 


* DM ODE 
8 1 NCONN 
PREPAR 
STIFFK 
WEIGHTS 

SCALE 

8EMODE 


8 1 NLAYR 
SAVE 


8INLAYR 


DMGFREQ 


LMPMAS 


8LAYPR 


8ELSTIF 

8LAYPR 

PRINT 

STRCON 


8 I NCONN 


LAYCALC 

TRQDOUT 


LAYCALC 


PLSTIF 

8LMPROD 

8PRINT 


8EMODE 

LMPMAS 

PRNTEL 

8TRECON 

8LAYPR 

8LAYPR 
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TRQDOUT 


AEYMIN  : 
AH  : 

AL  : 

ALS  : 
ALSTRS  s 
AM  : 

AM  AX  : 
AMAXAE  : 
AMX  : 
AREA 
AX  : 

AXX  : 
AY  : 

AYY  : 
AZ  : 

B  : 

BA  : 

BAE  : 
BASEA  i 


♦LMSI ZE 

♦INLAYR 

♦LMPROD 

* COORD 
STIFFK 

♦PREPAR 

♦INP03 

♦LMPMAS 

♦LMSIZE 

* AVER AGE 

* SCALE 

♦PLMASS 

* CHANGE 

♦TRECON 

♦LMSIZE 

♦TRECON 

♦LMSIZE 

♦ELSTIF 

*PREPAR 

♦PREPAR 

AVERAGE 

LMSIZE 

WEIGHTS 


PREPAR 

LMSIZE 

ELSTIF 

STRCON 

PRNTEL 

PREPAR 

♦QLSTRS 

♦DMODE 

* PL ST IF 
♦LMSIZE 

♦PREPAR 

♦PREPAR 

♦PREMULT 

♦TRQDSTR 

DMODE 

PRINT 


SAVE 

GETDK 
UN I TEG 

QLSTRS 

♦STRCON 

♦EMODE 

♦QDRLTL 

♦PREPAR 

STRESS 

STRESS 

♦TRNSFM 

EMODE 

PRNTEL 


LMPROD 


STRCON 


STRESS 


♦T RECON 


♦TRECGN 


LAYCALC 

♦SCALE 


PRNTEL 

STRESS 


♦TRECON 


LAYPR 

♦STRCON 
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BB 

• 

* 

aELSTIF 

BND 

m 

aF'RNTDR 

C 

m 

m 

ASEMBL 

ASEMDK 

CARD 

m 

m 

a  EC  HO 

Ct 

m 

■ 

ASEMBL 

aLMPROD 

ASEMDK 

STIFFK 

CCC 

• 

■ 

aTRNSFM 

CTA 

■ 

■ 

aTRECON 

D 

■ 

* 

a DAMAGE 
* REST OR 

aDMGITER 

DEFLMT 

DMODE 

aINDSPL 

DEFMAX 

aINDSPL 

DET 

a CON DNS 

DINCR 

DMODE 

aiNGNRL 

D I SP 

apRNTDR 

DK 

aASEMDK 

DMGITER 

DKF 

a  I NDMGE 

DKFCTR 

ASEMBK 
LIN  I  TEG 

a  I NDMGE 

DM 

aASEMDK 

FRO I TER 

DMFCTR 

ASEMDK 

a  I NDMGE 

DONE 

a INCONN 

aNODCHCK 

DR 

ANALYZ 

PPDISP 

CURRENT 

PRNTDR 

aELSTIF  GETDK  STIFFK 

aELSTIF  GETDK  aLMF'MAS 

aDMGORTH  aBX  aORTHOG 

SCALE 


DMGORTH 

F'RNTEL  ST  ROOM  TRQDSTR 

GEN MSS 

DMODE  aELFQRC  a  INLOADS 

F'RNTEL  SCALE  STRCGN 
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DRAT  10 

DTEXTRA 

DT I  TER 

DX 

DY 

DZ 

E 

El 

E2 

EDDR 

EDR 

EDR1 

EE 

EEK 

EEKK 

EEKM 

EEM 

EEMM 

EFFSTR 

EFSTRS 


LINT  FRO 
*SCALE 
* DAMAGE 
* DAM AGE 
* INXYZ 
■ttINXYZ 
*INXYZ 
#TRQDSTR 

ELSTIC 

UN I  TEG 

ELSTIC 

■aQLSTRS 

*ELFGRC 

UN I  TEG 

UN I TEG 
*ELSTIC 
#GETDK 
#ELSTIF 
*SUM 
*GETDK 
wELSTIF 
*QLSTRS 
wPRNTEL 


* WEIGHTS 


ELSTIF 

■frPREPAR 

PRNTEL 

UNTFRC 

UNTFRC 

PL ST IF 

#PREPAR 

•frPLST IF 

*PREPAR 

*LMPROD 

TRQDSTR 

QLSTRS 


#PREPAR 


*QLSTRS 

STRESS 

*ST IFFK 

*ST IFFK 
*PLMASS 

^STRESS 


PRNTEL 


STRCON 

*TRECGN 

#UNTFRC 


TRQDOUT 


STRCON 


TRQDSTR 


*TRQDSTR 
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Ek  : 

•CONDNS 

•UNTFRC 

•GETDK 

•QDRLTL 

•STIFFK 

UNITES 

EKK  : 

•CONDNS 

QDRLTL 

QLSTRS 

EKM  : 

•CHANGE 

•SUM 

•TRNSFM 

EL 

DECOUP 

•ORTHOG 

•DMGORTH 

•ERROR 

•FORSUB 

•LLT 

ELI  : 

* DECOUP 

ELASM  : 

• INP03 

ELCNST  : 

•  I NP03 

PREPAR 

ELEENG  : 

•PRNTEL 

•QLSTRS 

•STRCON 

TRQDOUT 

•TRQDSTR 

ELENG  : 

•PREPAR 

TRQDOUT 

•TRQDSTR 

ELENG1  : 

•UNITES 

ELENTH  : 

DMODE 

•STIFFK 

EMODE 

WEIGHTS 

•GETDK 

LMF'MAS 

LMSIZE 

ELL  : 

•DECOUP 

EM  : 

* CON DNS 

•GETDK 

•QDRLTL 

•STIFFK 

EMM  : 

* CON DNS 

•LMPMAS 

QDRLTL 

ENG  : 

•STRESS 

EN6CAP  : 

SCALE 

•STIFFK 

ENGG  : 

QLSTRS 

ENGLTA  : 

* DMODE 

•EMODE 

•STIFFK 

ENGMAX  : 

•TRQDSTR 

ENGST 1  : 

# DMODE 

UNI TEG 
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ENGSTR 

♦SCALE 

ENG TOT 

♦F'RNTEL 

ENGX 

LMSIZE 

ENGXY 

LMSIZE 

ENGY 

LMSIZE 

EONE 

♦PLSTIF 

ER 

♦ERROR 

ERR 

* ERROR 

ERRO 

* DAMAGE 

ERR1 

* DAMAGE 

ERREST 

♦DMGFREG 

ERRMAX 

♦ERR 

ERROR 

♦ERR 

ESRTIO 

♦F'REF'AR 

ETA 

* COORD 

ETWO 

♦PLSTIF 

EX 

♦PLSTIF 

EX  Y 

♦STRESS 

EY 

♦STRESS 

F 

♦REDUCE 

FCTR 

♦PRNTEL 

*STRCON  #TRQDSTR 
♦TRQDSTR  *UNITE6 
♦TRQDSTR  #UNITEG 
♦TRQDSTR  *UNITEG 

*  INCONN  *NODCHCK 


♦STRCON  TRQDOUT 
CRAMER 

♦STRESS 


♦STRCQN  * TRQDSTR 


♦TRQDSTR 


♦UNI  TEG 
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FR 

:  ANALYZ 

GENM 

:  *GENMSS 

GES 

:  * RANDOM 

gm 

:  ♦ASEMBL 

I 

:  ♦ASEMBL 

♦CON DNS 
♦DMGFREQ 
♦ELFORC 
♦ERROR 
♦GETDK 
* INLAYR 

♦LLT 

♦PLMASS 

♦PPMODE 

♦PRINTK 

♦QLSTRS 

♦SIVIB2 

♦TRNSFM 

♦WEIGHTS 


10 

• 

a 

♦BACKSUB 

♦STIFFK 

11 

■ 

■ 

♦BACKsub 

♦FQRSUB 

♦PPHDR 

12 

ELEMIN 

I2M1 

♦STRESS 

13 

ELEMIN 

14 

ELEMIN 

PRNTDR 

IA 

m 

u 

ABOUND 2 

IAA 

m 

m 

♦ASEMBL 

* INLOADS  PRNTDR 

PPMODE  PRNTMDS 


*B0UND2 

PREMULT 

♦ASEMDK 

♦BACKSUB 

♦CRAMER 

♦CURRENT 

♦DMGITER 

♦DMGORTH 

♦ELSTIC 

♦EL3TIF 

♦FQRSUB 

♦GENMSS 

*  INCONN 

♦INDMGE 

* INLOADS 

* I NLPMSS 

♦LMPMAS 

♦LMPROD 

♦PLST IF 

♦POP 

♦PPNODE 

♦PRELT 

♦PRNTDR 

♦PRNTEL 

♦RANDOM 

* REDUCE 

♦STIFFK 

♦STRESS 

♦TRQDOUT 

♦TRQDSTR 

♦FORSUB 

* INCONN 

♦B0UND2 

* DECOUP 

♦INDMGE 

♦INDSPL 

PRINT 

♦STIFFK 

INXYZ 

♦PPHDR 

INXYZ  ♦PPHDR 
INP03  INXYZ 

♦TRNSFM 

♦ASEMDK 


SCALE 


♦STIFFK 

♦B0UND2 

♦CHANGE 

* DAMAGE 

♦DECOUP 

♦DMODE 

♦  DX 

♦EMODE 

♦ERR 

♦GETBC 

♦GETBEST 

♦INDSPL 

♦INITIAL 

♦ INP03 

♦INXYZ 

♦LMSI ZE 

♦ORTHOG 

♦PPBISP 

♦PPELEM 

♦PREMULT 

♦PREPAR 

♦PRNTMDS 

♦QDRLTL 

* REST OR 

♦SCALE 

♦SUM 

♦TRECON 

♦UNITES 

♦UNTFRC 

♦POP 

♦PRELT 

ELEMIN 

♦ELSTIF 

♦LMPROD 

♦POP 

PRINT 

♦STRESS 

PRINT  PRNTDR 

♦PPHDR  PRINT 
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I AREAS  : 
IBC  : 
I BND  : 

I  COL  : 

ICTYF'E  : 
IDIAG  : 

I DKMK  ! 
I DMG  : 
IELNO  : 
IELF’R  : 
IELPRMl: 
I FLAG  : 
IH  : 

IHH  i 

I I  : 

IJ 

IL 

IL1  ! 
IM 

I  Ml  : 


I NCQNN 

*GETBC 

BOUNDS 

BACK SUB 
PRELT 

ftINXYZ 

ASEMBL 

*POP 

ftPRELT 

GETDK 

* INCONN 

*  INCONN 

ft  INCONN 

ft RANDOM 

ftBGUNDS 

ftSUM 

» DECOUP 
ftSTIFFK 

ftASEMDK 

ft  INLOADS 

ft DECOUP 

ft  I NDSPL 

ftBACKSUB 


ft INGNRL 

ftGETBC 

*B0UND2 

PREMULT 

BACKSUB 

PRELT 

ftINDMGE 

ftNODCHCK 


ft REDUCE 

ftFORSUB 

ftTRQDSTR 

DX 

ftPPDISP 

ft INLOADS 
*POP 


INLAYR 

PRNTDR 

FORSUB 

ft BOUNDS 
PREMULT 


ftRESTOR 


ftGETDK 

ftUNITEG 


REDUCE 

LLT 

FORSUB 

F'RINTK 


*SUM 


ftPLST IF 
ftUNTFRC 


RESTOR 

ftPOP 

LLT 

STIFFK 


ftPREMULT 
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I  MAX  : 

DAMAGE 

DMGFREQ 

•ERR 

IMODN  s 

•UN I TEG 

I  NCR  : 

•GETBC 

IND  : 

*T RECON 

INDANG  i 

* INGNRL 

INLAYR 

TRECON 

INDEX  : 
INDMIN  ! 

•PREF'AR 

INCONN 

•INGNRL 

INDX  : 

•UN I TEG 

INT  i 

•SI VIB2 

IP  : 

•ELSTIC 

IPi  : 

* DM G FREQ 

•FRQITER 

IPP  ! 

• INGNRL 
PF'TTL 

PPDISP 

PRNTMDS 

PPELEM 

PF'HDR 

PPNODE 

IPR  : 

* INCONN 

IQ 

* BACKS UB 

•FQRSUB 

•LLT 

•PREMULT 

ISQTRN  : 

•ELEMIN 

INP03 

PREF'AR 

I  SYS  : 

* INX YZ 

IT  : 

•BACKSUB 

•DAMAGE 

•DMGFREQ 

ITRI  i 

•POP 

I  TYPE  : 

*CTYPE 

•PPHBR 

•ELEMIN 

•PRNTEL 

•INCONN 

TRQDOUT 

•INXYZ 

•nodchck 

IX  ; 

•CHANGE 

•PLSTIF 

IY  : 

•CHANGE 

50 


J 

■ 

* 

♦ASEMBL 

♦CONDNS 

♦DMGITER 

♦ELSTIC 

♦FRQ ITER 

♦  I NDMGE 

♦LLT 

♦F'LST  IF 

♦PREMULT 

♦QDRLTL 

♦SCALE 

♦TRNSFM 

♦ASEMDK 

♦CRAMER 

♦DMGORTH 

♦ELSTIF 

♦GENMSS 

♦ INDSPL 

♦LMPMAS 

♦POP 

♦PREPAR 

♦QLSTRS 

♦SIVIB2 

♦TRQDSTR 

•JO 

■ 

♦  I  NDSF'L 

♦TRQDSTR 

Ji 

■ 

■> 

♦DMODE 

♦LMPRGD 

♦ELSTIF 

♦PRNTEL 

J3 

♦PRNTMDS 

J3M2 

♦PRNTMDS 

JA 

♦B0UND2 

♦TRNSFM 

JAA 

♦ASEMDK 

♦TRNSFM 

JD 

♦GETBC 

JH 

♦INP03 

♦SUM 

JJ 

♦INXYZ 

JJ1 

♦LLT 

JL 

* INP03 

JL1 

♦DECOUP 

♦ERROR 

JH 

♦ INDSPL 

♦INLOADS 

JNN 

m 

• 

♦PRNTDR 

♦BACKSUB 

♦B0UND2 

♦CHANGE 

♦DAMAGE 

♦DECOUP 

♦DMGFREQ 

♦DM ODE 

*DX 

♦ELFORC 

♦ERR 

♦ERROR 

♦FORSUB 

♦GETBC 

♦OETDK 

♦INCONN 

♦INLOADS 

♦INP03 

* INXYZ 

♦LMPROD 

♦ORTHOG 

♦PLMASS 

♦PPDISP 

♦PPMODE 

♦PRELT 

♦PRNTDR 

♦PRNTEL 

♦PRNTMDS 

♦RANDOM 

♦REDUCE 

♦RESTOR 

♦STXFFK 

♦STRCON 

♦SUM 

♦UN I  TEG 

♦UNTFRC 

♦WEIGHTS 

♦INCONN 

♦  INDSF'L 

♦INXYZ 

♦STRCON 

♦TRQDSTR 

♦UN I TEG 
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JX 

K 


KO 

K1 

K2 

KA 

KAA 

KANLYZE 

KDEFEQ 

KDMG 

KF 

KGEN 

KH 

KHH 
KI I 

KK 

KL 


♦ASEMBL 

♦BACK SUB 

♦DMGORTH 

♦GENMSS 

♦LMPROD 

♦PREMULT 

♦QLSTRS 

♦STRCON 

♦UN I TEG 

♦PRNTEL 

♦GETDK 

♦UNITED 

♦LMPMAS 

♦TRNSFM 

♦TRNSFM 

DOPROB 

* DMO BE 

♦ASEMDK 

♦  UN I  TEG 

* INCONN 

♦ASEMBL 
* I NDSPL 
♦UNITED 

♦ASEMBL 

♦ASEMDK 

♦ELFORC 

♦PRNTEL 


♦BQUND2 

♦B0UND2 

♦DX 

♦  I NDSPL 
♦ORTHGG 
♦PRINTK 
♦RANDOM 
♦STRESS 


♦STRCON 

♦LMPMAS 


* INGNRL 
UNITED 
*GETDK 


♦ASEMDK 

♦INLOADS 


♦ASEMDK 


♦LLT 


* CON DNS 

♦ELFORC 

* INLOADS 

♦PLSTIF 

♦PRNTDR 

♦REDUCE 

♦TRNSFM 


♦UN I TEG 
♦PRNTEL 


LAYPR 


♦BOUND 2 
*  I NP03 


♦PPDISP 


♦PREPAR 


♦DECOUP 

♦FOR SUB 

♦LLT 

♦PPMODE 

♦PRNTEL 

♦RESTOR 

♦TRQDOUT 


♦STRCON 


SCALE 


♦ELFORC 

♦PPDISP 


♦PRNTDR 


♦TRQDOUT 


♦DMGFREQ 

♦FRQITER 

♦LMPMAS 

♦PRELT 

♦PRNTMDS 

♦SCALE 

♦TRQDSTR 


♦TRQDSTR 


♦EMODE 

♦PRNTDR 
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KL1 

♦ DECOUP 

♦PREMULT 

KSAVE 

DOPROB 

♦WEIGHTS 

KSTR 

♦DAMOPT 

WEIGHTS 

KTR 

♦QLSTRS 

TRQDOUT 

KX 

♦ASEMBL 

♦ASEMDK 

♦DMODE 

♦ELFORC 

♦LMPMAS 

♦POP 

3TIFFK 

UNTFRC 

♦SUM 

KXX 

■ 

* 

♦BOUNDS 

KV 

a. 

■ 

♦ASEMBL 

♦ASEMDK 

GETDK 

♦INLOADS 

♦PRINTK 

♦QLSTRS 

UNTFRC 

L 

• 

W 

♦ASEMBL 

♦ASEMDK 

♦DMGITER 

♦DMODE 

♦ERR 

♦GETDK 

♦LMSIZE 

♦PLSTIF 

♦PRINT 

♦PRNTEL 

♦REDUCE 

♦RESTOR 

♦STRCON 

♦TRECON 

♦UNTFRC 

♦WEIGHTS 

LI 

to 

■ 

DECOUP 

♦SIVIB2 

ERROR 

LAM 

• 

A 

AVERAGE 

GETBEST 

LAYPR 

LMSIZE 

LAMFO 

* 

♦INCONN 

LAYERD 

■ 

to 

DMODE 

SAVE 

EMODE 

LDELTA 

to 

a 

♦PRNTMDS 

♦B0UND2 

♦CHANGE 

♦CONDNS 

GETDK 

♦INDSPL 

♦INLOADS 

♦PREPAR 

♦PRINTK 

♦QLSTRS 

TRQDOUT 

TRQDSTR 

♦UN I TEG 

♦BOUNDS 

♦CHANGE 

♦CONDNS 

♦LMPMAS 

♦POP 

♦PREPAR 

STIFFK 

♦SUM 

♦UNITED 

♦AVERAGE 

♦COORD 

♦DAMAGE 

♦DX 

♦ELFORC 

♦EMODE 

♦LAYCALC 

♦LAYPR 

♦LMPMAS 

♦POP 

♦PREMULT 

♦PREPAR 

♦PRNTMDS 

♦QDRLTL 

♦QLSTRS 

♦SAVE 

♦SCALE 

♦STIFFK 

♦TRQDOUT 

♦TRQDSTR 

♦UN I  TEG 

♦ORTHOG 

♦PRNTMDS 

RANDOM 

♦INCONN 

INLAYR 

LAYCALC 

PREPAR 

TRQDOUT 

♦TRQDSTR 

♦INGNRL 

INPT 

PRINT 
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LOSS 

LET 

LFLAG 

LINES 

LINNUM 

LM 

LMNODE 

LMTDSP 

LMTSTR 

LOADS 

LOCK  : 

LT 

M  : 

Ml 

M2  : 

MS  : 

MA 

MAA  : 

MAXDCCL: 


aELFGRC 

* ERROR 

*LAYCALC 

aINLAYR 

a  ECHO 

aLAYCALC 

alNLPMSS 

INDSPL 

aINGNRL 

ANALYZ 

PRNTEL 

WEIGHTS 

TERROR 

aLAYCALC 

aASEMDK 

aBACKSUB 

aASEMBL 

aTRNSFM 

ASEMBL 

INCONN 

ftSTRESS 

^INITIAL 

DOPROB 


apREPAR 

PRNTMDS 

LAYPR 

aLAYPR 

STIFFK 

*  I NGNRL 

SCALE 

CURRENT 

STIFFK 

aGRTHOG 

aDAMAGE 

aFORSUB 

aASEMDK 

ASEMDK 

ftPLSTIF 

aSUM 

QDRLTL 

* I NGNRL 


*SIVIB2 

apRNTDR 


INPT 

GETDK 

STRCON 

aSI V IB2 

aGETDK 

aLLT 

aTRNSFM 

COORD 

POP 

TRQDOUT 

QLSTRS 


a PRNTEL 


SCALE 


a  IN  LOADS 
TRQDSTR 


aPRNTMDS 

aPREMULT 

aCRAMER 

PPELEM 


INPT 

UNITES 


aUNITEG 

aUNITEG 

ELFORC 

PRNTEL 


54 


MAXDDIT: 

DAMAGE 

♦ INDMGE 

MAXDFITi 

DMGFREQ 

♦INDMGE 

MAXECCLi 

DOPROB 

♦ INGNRL 

MAXSK  : 

♦INITIAL 

INI  TP 

MAXSZE  : 

EMQDE 

♦INGNRL 

SCALE 

MB  : 

ASEMBL 

ELFORC 

PPELEM 

TRQDOUT 

ASEMDK 

♦INCONN 

PRNTEL 

♦CONDNS 

LMPMAS 

QDRLTL 

COORD 

♦PLSTIF 

♦STRESS 

♦CRAMER 

POP 

♦SUM 

MSB  : 

♦INITIAL 

QDRLTL 

QLSTRS 

MC  : 

ASEMBL 

ELFORC 

PPELEM 

ASEMDK 

♦INCONN 

QDRLTL 

♦CONDNS 

LMPMAS 

♦STRESS 

COORD 

♦PLSTIF 

♦SUM 

♦CRAMER 

POP 

TRQDOUT 

MCC  : 

♦INITIAL 

QDRLTL 

QLSTRS 

MD  s 

ASEMBL 

♦INCONN 

TRQDOUT 

ASEMDK 

LMPMAS 

♦CONDNS 

POP 

COORD 

PPELEM 

ELFORC 

QDRLTL 

MDEFEQ  : 

♦DMQDE 

UN I  TEG 

MEMBS  : 

AVERAGE 

GETBEST 

LAYPR 

PRNTEL 

UNTFRC 

CURRENT 

INCONN 

LMSIZE 

SAVE 

WEIGHTS 

DMODE 

INDMGE 

POP 

SCALE 

♦ELEMIN 

INLAYR 

PPELEM 

STIFFK 

EMODE 

INP03 

PRINT 

STRCON 

MESS  : 

♦NODCHCK 

MGEN  : 

♦INXYZ 

ML  : 

♦DAMAGE 

MLR*  : 

♦DAMAGE 
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MLR  J  : 
MM  : 

MM2  : 
MNLAYR  : 
MODES  : 

MSK  : 
MXMEMB  : 
M YOUNG  : 

N  : 
N1  : 
N2  : 
N3  : 
NA  : 
NAA  : 
NACTIVE: 
NBC1  : 
NBC2  : 
NBNDRY  : 


* DAMAGE 

ASEMBL 

INDSPL 

POP 

TRNSFM 

•PPMODE 

•INGNRL 

ASEMBL 

DMGFREQ 

GENMSS 

PRINT 

•DX 

•SCALE 

DMODE 

PRNTEL 

*DX 

•GETBC 

•GETBC 

•GETBC 

* ASEMBL 

•ASEMBL 

DMODE 

• GE  T  BC 

•GETBC 

B0UND2 
RES TOR 


ASEMDK 

•INGNRL 

PPDISP 


INLAYR 

ASEMDK 

DMGORTH 

GETDK 

PRNTMDS 


EMODE 

WEIGHTS 

•INXYZ 

•INXYZ 


* ASEMDK 


•INITIAL 


•GETBC 


ELFORC 

INLOADS 

PPMGDE 


BOUND 2 
ELSTIF 
* I NGNRL 
QDRLTL 


•INCONN 


•PRELT 

•LLT 


•SUM 


SCALE 


INPT 


ELSTIF 

INPT 

PRNTDR 


CON DNS 
ERROR 
PPMODE 
SIVIB2 


INLAYR 


• PR E MULT 


PRNTDR 


GETBC 

LMPROD 

STIFFK 


CURRENT 

FRQITER 

PREPAR 

STIFFK 


PREPAR 


REDUCE 
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NCASES  : 


□MODE 

SCALE 

UNTFRC 


*  I NPT 
STRCON 
WEIGHTS 


PRNTDR 

STRESS 


PRNTEL 

TRGDQUT 


QLSTRS 

TRQDSTR 


NCOMP  : 

♦ELEMIN 

NCGN  : 

♦ELFORC 

ND  : 

♦GETBC 

NO CASES: 

♦DMODE 

UN I TEG 

UNTFRC 

NDCYCL  : 

*DOPROB 

♦INITP 

WEIGHTS 

NDEFEQ  : 

♦SCALE 

NDK  s 

♦ASEMDK 

DX 

♦GETDK 

NDMG  : 

GETDK 

♦INDMGE 

NDMGCAS: 

CURRENT 

INDMGE 

PRNTEL 

WEIGHTS 

DAMAGE 

♦INGNRL 

PRNTMD3 

DMGFREQ 

I NLOADS 
STRCON 

DM  ODE 

I  NPT 
TRQDSTR 

GENMSS 
PPMODE 
UNI  TEG 

NDSP  : 

♦ELFORC 

NDUMMY  : 

♦DMODE 

DGPROB 

UN I TEG 

NECYCL  *. 

♦DGPROB 

♦1NITP 

WEIGHTS 

NELEM  : 

♦ELEMIN 

NFDEG  s 

♦LAYCALC 

LAYPR 

NFIL 

♦INITIAL 

F'PTTL 

PPDISP 

PPELEM 

PPHDR 

PPNODE 

NH  : 

♦B0UND2 

♦REDUCE 

♦RESTOR 

NH1  : 

♦REDUCE 
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MINOR 

* INXYZ 

NJLODS 

# IN LOADS 

NL1 

*ORTHOG 

^RANDOM 

NLAM 

*LAYF'R 

NLINES 

*ECHO 

NLMPMSS: 

*  I NGNRL 

INLPMSS 

NM 

BACKSUB 

DAMAGE 

DMGGRTH 

ERR 

GENMSS 

*INPT 

PREMULT 

SIVIB2 

PRINTK 

NMAT 

■ 

m 

#ELEM IN 

INP03 

NMT 

m 

m 

*EL.EMIN 

NN 

: 

BOUNDS 

DMODE 

*  I NPT 

PDF- 

SCALE 

WEIGHTS 

NND 

■ 

* 

*ASEMBL 

#TRNSFM 

fcASEMDK 

NNDEG 

• 

m 

#LAYCALC 

LAYPR 

N NODES 

m 

w 

ASEMBL 

ASEMDK 

EMODE 

GETDK 

POP 

PREPAR 

STIFFK 

ST ROOM 

UN I TEG 

UNTFRC 

NNRM 

m 

B 

*QDRLTL 

NO 

■ 

# CON DNS 

*QDRLTL 

NOD 

m 

m 

#  INCONN 

#NODCHCK 

I  NPT 

STIFFK 

WEIGHTS 

DECOUP 

ERROR 

LLT 

PRNTMDS 

DMGFREQ 

FORSUB 

OR T HOG 
RANDOM 

DMGITER 

FRQITER 

F’RELT 

SCALE 

DX 

F'PMODE 

INDSPL 

REDUCE 

INLOADS 
REST OR 

*ELFORC 

#LMF‘MAS 

* STRESS 

COORD 
# INCONN 
F'RNTEL 
*TRNSFM 
WEIGHTS 

DMODE 

INLAYR 

QDRLTL 

TRQDOUT 

ELFORC 

LMF'MAS 

QLSTRS 

TRQDSTR 
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NODES  : 

INDSPL 

PPNODE 

INPT 

PRNTDR 

NO  I 

•INGNRL 

SIVIB2 

NO NORM  : 

•PLST IF 

NONZRO  : 

INI  TP 

•POP 

NP  : 

•INXYZ 

NP1  J 

•GENMSS 

•PPMODE 

NPAGE  : 

•AVERAGE 
•EMODE 
* INDSPL 
•INLPMSS 
•POP 

•DAMAGE 

•GETBC 

•INGNRL 

•INP03 

•PRINTK 

NSTR  : 

DAMOPT 

•INITIAL 

NTRIAL  : 

DECOUP 

ERROR 

NUFR  : 

•DMGDE 

NWORK  : 

INXYZ 

NX  s 

*DMODE 

NY  : 

•RANDOM 

NZDEG  s 

•LAYCALC 

LAYPR 

PERMBA  : 

•SCALE 

PMU  s 

ELSTIC 

•PREPAR 

PMU1  : 

•ELSTIC 

POISON  ! 

•INP03 

PREPAR 

Q  : 

•DECOUP 

•INXYZ 

PRNTMDS 

PPDISP 

PPMODE 

STIFFK 

•PRNTMDS 

•DMGFREQ 

•GETBEST 

•INITP 

•INXYZ 

•PRNTDR 

•DMGDE 

•INCONN 

•INLAYR 

•LAYPR 

•PRNTEL 

•DOPROB 
•INDMGE 
» INLOADS 
•NODCHCK 
•PRNTMDS 

•ORTHOG 

RANDOM 

•SIVIB2 
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QUAD  : 

TRQDuUT 

TRQDSTR 

RO  : 

* DAMAGE 

R1  : 

♦  DAMAGE 
PRINT 

ELEMIN 

PRNTDR 

RAD  : 

♦INLAYR 

RHD  : 

ELSTIF 

LMPMAS 

RHG 1  : 

DM  DDE 

EMODE 

RHOA  : 

♦PLMASS 

S  : 

♦ERR 

*  INLOADS 

SAGE  : 

♦ CON DNS 

♦LMPMAS 

SHEARM  : 

♦ I NP03 

SK  : 

♦ASEMBL 

PRELT 

BACKSUB 

♦STIFFK 

SKGM  : 

♦PR I NTK 

SM  : 

ELSTIC 

♦PREPAR 

SPRDF  : 

DMODE 

EMODE 

SORT  1 2  : 

♦TRECON 

SSX  : 

♦QLSTRS 

TRQDOUT 

SSXY  : 

♦QLSTRS 

TRQDuUT 

SSY  : 

♦QLSTRS 

TRQDOUT 

STA  : 

♦TRECON 

ST R ENG  : 

DMODE 

♦EMODE 

STRMAX  : 

♦SCALE 

UNTFRC 

INP03  INXYZ  ♦PPHDR 

LMPRGD  PLMASS  ♦PREPAR 
♦INP03  PREPAR  WEIGHTS 

*UNITEG 

♦TRECGN 

♦BQUND2  FORSUB  ♦LLT 


* I NGNRL  STIFFK 


♦STRCGN  ♦TRQDSTR  ♦UN I TEG 
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sx 

#PRNTEL 

QLSTRS 

•STRCON 

•"■STRESS 

SXY 

QLSTRS 

•STRESS 

TRQDOUT 

SY 

QLSTRS 

# STRESS 

TRQDOUT 

T 

#DMGORTH 

*ERR 

TDR1 

*REST DR 

TBR2 

* RE ST OR 

TEMP 

*INDSPL 

•INLOADS 

•LMSIZE 

•NODCHCK 

TEXTRA 

#  DAMAGE 

TFFR 

GETDK 

TRQDSTR 

•LAYCALC 

UNTFRC 

•PREPAR 

STIFFK 

TFR  : 

*PREPAR 

PRNTEL 

TRQDOUT 

THCKMN  : 

* IN CONN 

THCKMX  : 

* INCONN 

THICK  : 

* I NCONN 

THKLAM  : 

*  I NGNRL 

INLAYR 

LAYCALC 

LAYPR 

TIMEO  : 

#ANALYZ 

* DAMAGE 

•DMGFREQ 

•SI VIB2 

TIME1  s 

*ANALYZ 

•DAMAGE 

•DMGFREQ 

•SI VIB2 

TIME2  : 

•ft-ANALYZ 

* DAMAGE 

TIMES  : 

#ANALYZ 

TITER  : 

•DAMAGE 

TITLE  s 

AVERAGE 

EMODE 

INDSPL 

DAMAGE 
GETBC 
*  I NGNRL 

DMGFREQ 

GETBEST 

INLAYR 

DMO  DE 

INCONN 

INLOADS 

TRQDOUT 


TRQDOUT 


TRQDSTR 


DOPROB 

INDMGE 

INLPMSS 


61 


TOLDMGDs 
TOLDMGF: 
TOLVEC  i 
TRANG  : 
TRIANG  : 
TTHK  : 

TWOPII  : 
U  : 

UDR 

UV  : 

v  ; 

w  : 

WBAR  : 
WEIGHT  : 
WLMPMSS: 
WMAX  : 
WMEMB  : 


I NP03 

PPMODE 

PRNTMDS 

DAMAGE 

DMGFREG 

DECOUP 

LMPMAS 

♦CRAMER 

♦GETDK 

♦UNTFRC 

♦PRNTMDS 

♦BACKSUB 

FRQITER 

♦DMOBE 

♦STRESS 

♦ DAMAGE 
GENMSS 

♦BACKSUB 

GENMSS 

SIVIB2 

♦ WEIGHTS 

DGPROB 

♦INLPMSS 

♦PRNTMDS 

♦WEIGHTS 


INXYZ 

PPTTL 

*INDMGE 
* INDMGE 
ERROR 
♦QDRLTL 
QDRLTL 
PLMASS 

* DAMAGE 
GENMSS 

UNTFRC 

♦ DECOUP 
SIVIB2 

* DECOUP 
♦ORTHOG 

SAVE 

STIFFK 


LAYPR 

PRINTK 


* INGNRL 
TRQDSTR 
STRESS 
PLSTIF 

* DECOUP 
♦PLSTIF 


DMGFREQ 


♦DMGFREQ 

PPMODE 

♦WEIGHTS 

WEIGHTS 


NODCHCK 
PR NT DR 


TRQDOUT 

♦STIFFK 


♦DMGORTH 

♦SIVIB2 


♦FOR SUB 


ERROR 

♦PRNTMDS 


POP 

F'RNTEL 


UNTFRC 

♦TRQDSTR 


ERROR 


♦FRQITER 


FRQITER 

♦RANDOM 
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WSHEAR  : 

•WEIGHTS 

WTLAST  : 

DOPRQB 

alNITP 

a£AVE 

WEIGHTS 

WW  : 

•WEIGHTS 

X  s 

COORD 

PPMODE 

aDMGITER 

PPNODE 

aFORSUB 

aPRELT 

INLOADS 
PR NT DR 

aiNXYZ 

PRNTMDS 

XI 

aCTYPE 

aiNXYZ 

X2  ! 

aCTYPE 

aiNXYZ 

X3  : 

•CTYPE 

aiNXYZ 

XA  : 

*1 NLAYR 

XANG  i 

aiNLAYR 

TRECON 

XCGMP  : 

* COORD 

XI  : 

aCOORD 

CRAMER 

XL  : 

•DECOUP 

aPRNTMDS 

DMGFREQ 
•SI VIB2 

ERROR 

aFRQITER 

PPMODE 

XP  : 

aiNXYZ 

Y  : 

aDMGITER 
a PREMULT 

a DM G ORTH 
a RAN DOM 

aDX 

FRQITER 

apRELT 

Y1  : 

aCTYPE 

Y2  : 

•CTYPE 

Y3  : 

aCTYPE 

YA  : 

aiNLAYR 

YANG  : 

aiNLAYR 

TRECON 

YCOMP  : 

aCOORD 
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YOUNG M  : 

*INP03 

PREPAR 

YP  ; 

*INXYZ 

Z  : 

*DX 

*PREMULT 

ZA  s 

# INLAYR 

ZANG  : 

*INLAYR 

TRECON 

Z  COMP  : 

* COORD 

ZP 

* INXYZ 
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APPENDIX  B 

COMMON  BLOCK  /  SUBROUTINE  NAME 
CROSS-REFERENCE  TABLE 
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ALSTRS  s 

F'REF'AR 

TRQDOUT 

PRNTEL 

AMAXMINi 

AVERAGE 

INLAYR 

DMODE 

SCALE 

ANG  : 

INGNRL 

INLAYR 

AREA  i 

COORD 

LMF'ROD 

QDRLTL 

TRQDSTR 

ELSTIF 
PLMASS 
STIFFK 
UN I TEG 

BASEA  : 

AVERAGE 
GET BEST 
LAYCALC 
PRINT 
STRCON 

CURRENT 

INCONN 

LAYF'R 

PRNTEL 

UNTFRC 

END  i 

B0UND2 

RESTOR 

GETBC 

CONN  : 

ASEMBL 

LMPMAS 

TRQDOUT 

ASEMDK 

PUP 

DEFLMT  : 

DMODE 

INDSPL 

DK  : 

ASEMDK 

GENMSS 

DMG ITER 
GETDK 

DMG  : 

ASEMDK 

GENMSS 

INPT 

TRQDSTR 

CURRENT 
GETDK 
PPM ODE 
UNITED 

DMGERR  : 

DAMAGE 

DMGFREQ 

DUMMY  : 

DMODE 

UNITED 

EIGN  : 

ASEMBL 

DECOUP 

FRQITER 

ASEMDK 

DMGFREQ 

GENMSS 

QLSTRS 

STRCON 

STRESS 

EMODE 

I NCONN 

INGNRL 

TRECDN 

UNTFRC 

GETDK 

PLSTIF 

STRCON 

UNTFRC 

LAYCALC 

PREPAR 

STRESS 

LMPMAS 

PRNTEL 

TRQDOUT 

DMODE 

INDMGE 

LMSIZE 

SAVE 

WEIGHTS 

ELEMIN 

INLAYR 

POP 

SCALE 

EMODE 

INP03 

PPELEM 

STIFFK 

INPT 

PRNTDR 

REDUCE 

COORD 

PPELEM 

ELFORC 

PRNTEL 

INCONN 

QDRLTL 

INGNRL 

SCALE 

DMGORTH 

DX 

FRQITER 

DAMAGE 

INDMGE 

PRNTEL 

WEIGHTS 

DMGFREQ 

INGNRL 

PRNTMDS 

DMODE 

I NLOADS 
STRCON 

I NDMGE 

UNTFRC 

BOUNDS 

DMGORTH 

GETDK 

CONDNS 

ELSTIF 

INGNRL 

CURRENT 

ERROR 

PPMODE 
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F'REPAR 

SIVIB2 

EIGVEC  s 

DAMAGE 

FRQITER 

ELEM  : 

ASEMBL 

DMODE 

I NCGNN 
LMSIZE 
QDRLTL 
STRCON 
WEIGHTS 

ENERGY  s 

PREPAR 

TRQDSTR 

ESTRESSs 

PREPAR 

TRQDOUT 

FRDR  : 

ANALYZ 

INPT 

SCALE 

TRQDSTR 

GENM  : 

GENMSS 

GM  ! 

ASEMBL 

LAYENG  ! 

LMSIZE 

LAYERD  * 

AVERAGE 

INLAYR 

PREPAR 

LAYMIN  : 

INGNRL 

TRQDSTR 

LAYPRNT: 

LAYCALC 

LMASS  : 

INGNRL 

LMTEXCD: 

DMODE 

■> 


PRINT 

STIFFK 

PRNTMDS 

DECOUP 

GENMSS 

DMGFREQ 

PPMODE 

ASEMDK 

ELFORC 

INLAYR 

POP 

QLSTRS 

TRQDOUT 

AVERAGE 

EMODE 

LAYCALC 

PREPAR 

SAVE 

TRQDSTR 

PRNTEL 

QLSTRS 

PRNTEL 

TRQDSTR 

QLSTRS 

CURRENT 

PPDISP 

stiffk: 

UN I TEG 

DMODE 

PRNTDR 

STRCON 

UNTFRC 

PPMODE 

PRNTMDS 

BGUND2 

PREMULT 

TRQDSTR 

UNITES 

DMODE 

INPT 

PRINT 

EMODE 

LAYCALC 

PRNTEL 

INLAYR 

LAYCALC 

LAYPR 

INLPMSS 

INPT 

DOPROB 

INITIAL 

QDRLTL 

RANDOM 

DMGORTH 

PRNTMDS 

ERROR 

SIVIB2 

COORD 

GETBEST 

LAYPR 

PRINT 

SCALE 

UNITEG 

CURRENT 

GETDK 

LMPMAS 

PRNTEL 

STIFFK 

UNTFRC 

STRCON 

TRQDOUT 

STRCON 

STRESS 

GETDK 

PRNTEL 

STRESS 

WEIGHTS 

INLOADS 

QLSTRS 

TRQDOUT 

STIFFK 

GETBEST 

LAYPR 

SAVE 

INGNRL 

LMSIZE 

TRQDOUT 

LAYPR 

LMSIZE 

STIFFK 

WEIGHTS 

SCALE 

UNITEG 
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UNTFRC 

LOCAL  : 

COORD 

LMPROD 

CRAMER 

TRECON 

MASS  : 

ASEMBL 

LMPMAS 

ASEMDK 

LMPROD 

MAT  : 

DMODE 

EM  ODE 

MATAXIS: 

ELSTIC 

PLST IF 

N  : 

ASEMBL 

DECOUP 

DX 

FOR SUB 

I NGNRL 

ORTHOG 

PREMULT 

REDUCE 

TRNSFM 

ASEMDK 

DMGFREQ 

ELFORC 

FRQITER 

INLOADS 

POP 

PRINTK 
RE ST OR 
WEIGHTS 

NMAT  : 

ELEMIN 

INP03 

NODES  : 

COORD 

INXYZ 

PRNTDR 

DMGFREQ 

MESHG 

PRNTMDS 

NPROB  ; 

DAMOPT 

INITIAL 

OPT  : 

DOPROB 

LAYPR 

INDSPL 

SCALE 

POST PR  : 

I NGNRL 
PPNODE 

INITIAL 

PPTTL 

QUAD  s 

INITIAL 

QDRLTL 

SAVE  : 

AVERAGE 

PRNTEL 

GETBEST 

SAVE 

SK  : 

ASEMBL 

PRELT 

BACKSUB 

STIFFK 

ELFORC 

ELSTIF 

LMPMAS 

TRNSFM 

CONDNS 

ELSTIF 

GETDK 

PREPAR 

QDRLTL 

STIFFK 

INF' 03 

PREPAR 

WE  I GHTS 

PREPAR 

STRESS 

TRECON 

BACKSUB 

BOUND 2 

DAMAGE 

DMGITER 

DMGORTH 

DMODE 

ELSTIF 

ERR 

ERROR 

GENMSS 

GETBC 

I  NDSF'L 

INPT 

LLT 

LMPROD 

F'PDISP 

F'PMODE 

PRELT 

PRNTDR 

PRNTMDS 

RANDOM 

SCALE 

SIVIB2 

STIFFK 

PREPAR 

INDSPL 

INLOADS 

INPT 

F'PDISP 

F'PMODE 

PPNODE 

INITP 

WEIGHTS 

I NGNRL 

INITP 

INPT 

WEIGHTS 

F'PDISP 

PPELEM 

F'PHDR 

PRNTMDS 

QLSTRS 

LAYCALC 

LAYPR 

PRINT 

TRQDOUT 

B0UND2 

FORSLIB 

LLT 
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SKGM  : 

ASEMBL 

LLT 

STIFFK 

STIFF  : 

ASEMBL 

PREPAR 

UNTFRC 

STRESS  : 

F'RNTEL 

STRNENGs 

BMODE 

STRCON 

TEMP  ! 

ELSTIC 

PREPAR 

TITLE  : 

AVERAGE 

EMODE 

INDSPL 

1NLPMSS 

POP 

F'RNTEL 

WEIGHT  : 

DOPROB 

BACKSUB 

BOUND 2 

POP 

PRELT 

ASEMDK 

CONBNS 

QDRLTL 

QLSTRS 

QLSTRS 

STRCON 

EMODE 

INGNRL 

TRQDSTR 

UNITED 

ELSTIF 

LMPMAS 

F'RNTEL 

STRCON 

DAMAGE 

DMGFREQ 

GETBG 

GETBEST 

INGNRL 

INITP 

I  NF'03 

INXYZ 

PPMODE 

PRNTMDS 

F'PTTL 

INITP 

SAVE 

FORSUB 

PREMULT 

INITP 

F'RINTK 

ELSTIF 

STIFFK 

GETDK 

UN I TEG 

STRESS 

TRQDOUT 

SCALE 

STIFFK 

LMPROD 

TRQDOUT 

PLMASS 
UN I  TEG 

DMODE 

INCONN 

INLAYR 

LAYPR 

F'RINTK 

DOPROB 

INDMGE 

INLOADS 

NODCHCK 

PRNTDR 

WEIGHTS 
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APPENDIX  C 

SUBROUTINE  DESCRIPTIONS 
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ROUTINE  NAME 
PURPOSE 

CALL  SEQUENCE  - 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS  - 
FILE  NAMES 


NOTES 


Main  Program 
Initiate  Calculations 
N/A 
N/A 


N/A 

INITIAL,  D0PR0B 
NPR0B 

DP0ST ,  0UTPUT ,  TAPES,  TAPE 6 ,  TAPE 8 ,  TAPE 9 9 


The  call  to  INITIAL  sets  certain  constants. 

The  D0PR0B  call  is  in  a  loop  which  allows  several 
problems  to  be  run  at  once.  It  is  usually  more 
convenient  to  keep  the  data  sets  separate  and 
run  only  one  problem  at  a  time. 
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ROUTINE  NAME  - 
PURPOSE 

CALL  SEQUENCE  - 
ARGUMENTS 


CALLED  BY 
EXTERNALS 

COMMON  BLOCKS  - 
FILE  NAMES 


NOTES 


ANAL Y  Z 

Perform  Static  Analysis  of  Undamaged  Design 

CALL  ANALYZ 

None 


CURRENT 

BACKSUB,  BOUND 2 ,  FORSUB,  LLT,  REDUCE,  SECOND, 
STIFFK 

FRDR 

TAPE  6 


Calls  all  the  routines  which  form  the  mass  and 
stiffness  matrices,  imposes  the  boundary 
conditions,  decomposes  the  stiffness  matrix, 
and  solves  for  the  displacements. 
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ROUTINE  NAME  - 

ASEMBL 

PURPOSE 

Assemble  the  Global  Mass  and  Stiffness  Matrices 

CALL  SEQUENCE  - 

CALL  ASEMBL (L) 

ARGUMENTS 

L  =  Element  Number 

CALLED  BY 

STIFFK 

EXTERNALS 

PUTSK,  SK 

COMMON  BLOCKS  - 

C0NN,  EIGN,  ELEM,  GM,  MASS,  N,  SKGM,  STIFF 

FILE  NAMES 

None 

NOTES 

None 
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ROUTINE  NAME  - 

ASEMDK 

PURPOSE 

Assemble  the  Damaged  Mass  and  Stiffness  Matrices 

CALL  SEQUENCE  - 

CALL  ASEMDK (L,M,KDMG) 

ARGUMENTS 

L  =  Damaged  Element  Number 

M  =  Damage  Case 

KDMG  =  Index  of  Damaged  Element  in  DMFCTR  Array 

CALLED  BY 

GETDK 

EXTERNALS 

None 

COMMON  BLOCKS  - 

C0NN,  DK,  DMG,  EIGN,  ELEM,  MASS,  N,  STIFF 

FILE  NAMES 

None 

NOTES 

None 
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ROUTINE  NAME  - 

AVERAGE 

PURPOSE 

Average  Current  Design  with  Least  Weight  Design 

CALL  SEQUENCE  - 

CALL  AVERAGE 

ARGUMENTS 

None 

CALLED  BY 

D0PR0B 

EXTERNALS 

None 

COMMON  BLOCKS  - 

AMAXMIN,  BASEA,  ELEM,  LAYERD,  SAVE,  TITLE 

FILE  NAMES 

TAPE  6 

NOTES 

None 
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ROUTINE  NAME  - 

BACKSUB 

PURPOSE 

Backsubstitution 

CALL  SEQUENCE  - 

CALL  BACKSUB (W,U,I0,I1) 

ARGUMENTS 

W  =  Solution  Vector 

U  =  Right  Hand  Side 

10  -  First  Vector  to  be  Used 

11  =  Last  Vector  to  be  Used 

CALLED  BY 

ANALYZ ,  DMGITER,  DM0DE ,  FRQITER,  SIVIB2 

EXTERNALS 

None 

COMMON  BLOCKS  - 

N,  SKBM 

FILE  NAMES 

none 

NOTES 

T 

L  W  =  U  is  solved.  The  W  and  U  need  not  be 

in  separate  locations. 
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ROUTINE  NAME 
PURPOSE 

CALL  SEQUENCE 
ARGUMENTS 


BEAMKM 

Compute  Three-dimensional  Beam  Element  Mass  and 
Stiffness 

CALL  BEAMKM ( ITYPE, ITRAN, EM, CC, XL) 

ITYPE  =  0  For  Stiffness  Calculations 
>  0  For  Mass  Calculations 

ITRAN  =  0  For  3-D  Beam 

=  1  For  Plane  Beam 

EM  =  Element  Matrix 

CC  =  Axial  Strains 

XL  =  Length  of  Beam 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


GETDK,  STIFFK 
BMSYMM,  BMX 
BMFR0B,  TEMP 
None 


NOTES 


See  Vol.  I  for  discussion  of  beam  element. 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


NOTES 


BMGE0M 

Geometry  Calculation  for  Beam 
CALL  BMGE0M (TR) 

TR  =  Geometry  Matrix 


BMX 

None 

NV 

None 


None 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


NOTES 


BMSYMM 

Make  the  Beam  Element  Matrices  Symmetric 
CALL  BMSYMM ( EM, I0RDER) 

EM  =  Mass  or  Stiffness  Matrix 
I0RDER  =  Array  Containing  DOF  Ordering 


BEAM KM 
None 
None 
None 


None 
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ROUTINE  NAME 

PURPOSE 

CALL  SEQUENCE 

ARGUMENTS 

-  BMX 

Sort  the  Beam  DOF's 

-  CALL  BMX ( EM, I0RDER) 

EM  =  Mass  or  Stiffness  Matrix 
I0RDER  =  Array  Containing  DOF  Ordering 

CALLED  BY 

EXTERNALS 

COMMON  BLOCKS 

-  BEAM KM 

-  BMGE0M 

None 

NOTES 

None 
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ROUTINE  NAME 

-  BMXTRS 

PURPOSE 

-  Compute  Beam 

Stresses 

CALL  SEQUENCE 

-  CALL  BMXTRS (II, XL) 

ARGUMENTS 

-  II  =  Element 

XL  =  Length 

Number 

CALLED  BY 

-  PRNTEL 

EXTERNALS 

-  BMGE0M 

COMMON  BLOCKS 

-  C0NN,  ELEM,  FRDR,  MAT,  N,  N0DES,  TTBME 
AMAXMIN,  BASEA,  AFK 

FILE  NAMES 

None 

NOTES 


None 


ROUTINE  NAME  - 

B0UND2 

PURPOSE 

Apply  Boundary  Conditions  to  the  Mass  and  Stiffness 

CALL  SEQUENCE  - 

CALL  B0UND2 

ARGUMENTS 

None 

CALLED  BY 

ANALYZ 

EXTERNALS 

PUTSK,  SK 

COMMON  BLOCKS  - 

BND ,  EIGN,  GM,  N,  SKGM 

FILE  NAMES 

None 

NOTES 

The  rows  and  columns  in  the  mass  and  stiffness 

matrices  which  correspond  to  the  boundary  degrees 

of  freedom  are  eliminated. 
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ROUTINE  NAME  - 
PURPOSE 

CALL  SEQUENCE  - 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS  - 
PILE  NAMES 


NOTES 


CHANGE 

Perform  Row  and  Column  Interchanges  for 
Quadrilateral  Elements 

CALL  CHANGE (EKM, IX, IY) 

EKM  =  Element  Mass  or  Stiffness  Matrix 

IX,  IY  =  Row  and  Column  Numbers  to  be  Interchanged 


C0NDNS 

None 

None 

None 


The  rows  and  columns  of  a  quadrilateral  element 
mass  or  stiffness  matrix  must  be  interchanged 
so  that  the  degrees  of  freedom  are  in  ascending 
order  as  required  by  ASEMBL. 
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ROUTINE  NAME  - 
PURPOSE 

CALL  SEQUENCE  - 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS  - 
FILE  NAMES 


NOTES 


C0NDNS 

Condense  the  Quad.  Mass  and  Stiffness  Matrices 
CALL  C0NDNS (MB , MC ,MD , N0) 

MB,  MC,  MD  =  Node  Numbers 
N0  =  Control  Parameter 


QDRLTL 

CHANGE 

EIGN,  MASS,  STIFF 
None 


Since  the  quadrilateral  elements  are  formed  by 
assembly  of  four  triangles,  the  center  node 
degrees  of  freedom  must  be  reduced.  The  10  x  10 
mass  and  stiffness  matrices  are  reduced  to 
8  x  8  by  static  condensation. 
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ROUTINE  NAME  - 

C00RD 

PURPOSE 

Determine  Direction  Cosines  for  Element 
Transformation 

CALL  SEQUENCE  - 

CALL  COORD (L) 

ARGUMENTS 

L  =  Element  Number 

CALLED  BY 

GETDK ,  PRNTEL,  STIFFK,  STRC0N ,  UNTFRC 

EXTERNALS 

SQRT 

COMMON  BLOCKS  - 

AREA,  C0NN ,  ELEM,  L0CAL,  NODES 

FILE  NAMES 

None 

NOTES 

None 
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ROUTINE  NAME  - 

CRAMER 

PURPOSE 

Cramer's  Rule  for  Inversion  of  Triangle 
Coordinate  Matrix 

CALL  SEQUENCE  - 

CALL  CRAMERfA,  TRIANG ,  MA ,  MB  ,  MC) 

ARGUMENTS 

A  =  3  x  3  Matrix  which  Contains  the  Inverse 

TRIANG  =  Area 

MA,  MB,  MC  =  Node  Numbers 

CALLED  BY 

PLSTIF ,  STRESS 

EXTERNALS 

None 

COMMON  BLOCKS  - 

LOCAL 

FILE  NAMES 

None 

NOTES 

None 
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ROUTINE  NAME 
PURPOSE 

CALL  SEQUENCE  - 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAME 


NOTES 


CTYPE 

Transformations  to  Cartesian  Coordinates 

CALL  CTYPE ( ITYPE, Xl ,X2,X3) 

ITYPE=1  Cylindrical  coordinates  given 
-2  Spherical  coordinates  given 

X1,X2,X3  =  r,0,z  for  ITYPE=1 
=  for  ITYPE=2 


INXYZ 

None 

None 

None 


Other  transformations  can  be  supplied  by  the 
user.  On  return,  X1,X2,X3  are  the  x,y,z 
cartesian  coordinates . 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


NOTES 


CURRENT 

Evaluate  Current  Design  and  Call  Related  Routines 

CALL  CURRENT 

None 


D0PR0B,  LAYPR 

ANALYZ,  DAMAGE,  DMGFREQ,  SCALE,  SIVIB2,  WEIGHTS 

BASEA,  DMG ,  EIGN,  ELEM,  FRDR 

TAPES 


The  current  design  is  evaluated  from  a  stress 
and  deflection  standpoint.  If  required,  the 
natural  frequencies  and  modes  will  be  calculated. 
The  structural  weight  is  also  computed. 
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ROUTINE  NAME  - 
PURPOSE 

CALL  SEQUENCE  - 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS  - 
FILE  NAMES 


NOTES 


DAMAGE 

Perform  Static  Reanalysis  for  Damage  Cases 
CALL  DAMAGE (L,D) 

L  =  Load  Case  Counter 

D  =  Array  of  Deflections  for  Each  of  the 
Damage  Cases 


CURRENT,  DM0DE 
DMGITER,  ERR,  GETDK ,  SEC0ND 
DMG,  DMGERR,  EIGVEC,  N,  TITLE 
TAPE  6 


Main  iteration  loop  contains  the  Aitken's 
extrapolation  procedure.  Time  and  error 
summary  printed  at  the  end  of  the  loop. 
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ROUTINE  NAME  - 
PURPOSE 

CALL  SEQUENCE  - 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS- 
FILE  NAMES 


NOTES 


DEC0UP 

Decouple  the  Eigenvector  Matrix 

CALL  DEC0UP 

None 


SIVIB2 

SQRT 

EIGN,  EIGVEC ,  N 
None 


The  eigenvalues  are  also  sorted  in  descending 
order  of  magnitude.  The  two  most  recent 
eigenvectors  are  decoupled. 
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ROUTINE  NAME  - 
PURPOSE 

CALL  SEQUENCE  - 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS  - 
FILE  NAMES 


NOTES 


DMGFREQ 

Frequency  Reanalysis  for  Damage  Cases 

CALL  DMGFREQ 

None 


CURRENT 

DMG0RTH,  ERR,  FRQITER,  GETDK ,  SEC0ND 
DMG,  DMGERR ,  EIGN,  EIGVEC,  N,  N0DES ,  TITLE 
TAPE  6 


The  damaged  stiffness  and  mass  matrices  are 
obtained  and  the  vector  estimate  is  orthogonalized 
with  respect  to  the  modified  stiffness  matrix. 
Procedure  is  repeated  until  the  vector  error  in 
the  norm  is  small. 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


DMGITER 

Damage  Reanalysis  Iteration 
CALL  DMG ITER(X,YfD,L) 

X  =  Old  Iterate 
Y  =  New  Iterate 

D  =  Original  Displacement  Vectors 
L  =  Load  Case  Counter 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


DAMAGE 

BACKSUB,  DX,  F0RSUB 

DK,  N 

None 


NOTES 


Iteration  formula  is 


K  dx 
x 


new 

new 


aK  xold 
x  +  dx 


new 


response  in  damage  condition 


where  x  is  the  solution  of  the  undamaged 
problem  Kx  =  f.  If  loads  change  in  damage 
condition  f  =  f  -  f^,  then  modify  the  Y 
vectors  output  just  after  call  to  DX.  Set 
Y  =  Y  -  f ^  before  call  to  F0RSUB. 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
PILE  NAMES 


NOTES 


DMG0RTH 

Orthogonalize  Current  Eigenvector  Estimates 
CALL  DMG0RTH (T) 

T  =  Array  of  Vectors  to  be  Orthogonal i zed 


DMGFREQ,  FRQITER 
DX,  PRELT,  SQRT 
DK,  EIGN,  EIGVEC ,  N 
None 


The  current  eigenvector  estimates  for  the  damage 
cases  are  orthogonalized  with  respect  to 
K-dK  where  K  =  LLT. 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 

COMMON  BLOCKS 

FILE  NAMES 


NOTES 


DM0  DE 

Resize  in  the  Displacement  Mode 

CALL  DM0DE 

None 


D0PR0B 

BACKSUB,  DAMAGE,  F0RSUB ,  LMSIZE,  REDUCE,  REST0R, 
SQRT,  UN 

AMAXMIN,  BASEA,  DEFLMT ,  DMG ,  DUMMY,  ELEM,  FRDR, 
LAYERD,  LM 

TAPE  6 


A  check  is  first  made  to  see  if  there  are  any 
active  displacements  (i.e.  ones  sufficiently 
close  to  their  limits) .  Next  a  dummy  loads 
vector  is  formed  for  each  active  displacement. 
The  virtual  displacements  due  to  these  dummy 
loads  are  then  computed.  Static  reanalyses 
for  the  damage  cases  are  then  performed. 
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ROUTINE  NAME 
PURPOSE 

CALL  SEQUENCE  - 
ARGUMENTS 


CALLED  BY 
EXTERNALS 

COMMON  BLOCKS  - 


NOTES 


D0PR0B 

Overall  Sequence  Control 

CALL  D0PR0B 

None 


Main  Program 

INITP ,  CURRENT,  EM0DE ,  SAVE,  GETBEST,  DM0DE , 
AVERAGE,  PRINT 

LMTEXCD ,  0PT,  TITLE,  WEIGHT 


Routines  are  called  which  read  and  check  the 
input  data.  The  current  design  is  then  analyzed. 
Resizing  is  performed  in  the  energy  mode  and 
then  in  the  displacement  mode.  If  the  weight 
increases  during  the  next  energy  mode  resizing 
the  best  design  is  returned  and  the  program 
proceeds  with  the  next  step  in  the  optimization. 
If  the  weight  increases  during  the  displacement 
mode  resizing,  this  design  is  either  averaged 
with  the  current  least  weight  design  or  resizing 
is  terminated  if  the  weight  has  more  than 
doubled . 
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ROUTINE  NAME 

DX 

PURPOSE 

Premultiplication  by  Damaged  Mass  or  Stiffness 
Matrix 

CALL  SEQUENCE  - 

CALL  DX(Y,D,Z,N) 

ARGUMENTS 

Y  =  Vectors  Resulting  from  Premultiplication 

D  =  Damaged  Mass  or  Stiffness  Matrix 

Z  =  Vector  Premultiplied  by  D 

N  =  Number  of  Vectors 

CALLED  BY 

EXTERNALS 

COMMON  BLOCKS  - 

FILE  NAMES 

DMGITER,  DMG0RTH,  FRQITER,  GENMSS 

REDUCE,  REST0R 

DK,  N 

None 

NOTES 

None 
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ROUTINE  NAME 
PURPOSE 

CALL  SEQUENCE  - 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


NOTES 


ECH0 

Echo  the  Input  Data 

CALL  ECH0 

None 


INITIAL 

E0F 

None 

TAPES,  TAPE 6 


The  input  data  is  printed  at  the  beginning  of 
each  run  from  TAPES.  Fifty  lines  are  printed 
per  page  with  column  and  line  numbers . 
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ROUTINE  NAME 
PURPOSE 

CALL  SEQUENCE  - 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS  - 
FILE  NAMES 


NOTES 


ELEMIN 

Read  Element  Data 
CALL  ELEMIN 
None 


INPT 

INP03 , PPHDR 
BASEA,  NMAT 
TAPES 


Read  element  control  data  and  call  INP03 
which  reads  connectivity  and  materials  data. 
Several  key  variables  are  defined: 

NMAT  -  NMT  =  total  number  of  materials 

IS0TRN  =  NMAT-NC0MP  =  number  of  isotropic 

materials 

MEMBS  =  NELEM  -  number  of  elements 
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ROUTINE  NAME 

-  ELF0RC 

PURPOSE 

Get  Element  Displacements  and  Transform 
Local  System 

CALL  SEQUENCE 

-  CALL  ELF0RC (DR, EDR, L, LDS) 

ARGUMENTS 

DR  =  Global  Displacement  Array 

EDR  =  Local  Displacement  Array 

L  -  Element  Number 

LDS  =  Number  of  Load  Cases 

CALLED  BY 

-  PRNTEL,  STRC0N,  UNTFRC 

EXTERNALS 

None 

COMMON  BLOCKS 

-  C0NN ,  ELEM,  L0CAL,  N 

FILE  NAMES 

-  None 

NOTES 

None 
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ROUTINE  NAME 

-  ELSTIC 

PURPOSE 

Generate  Elastic  Constant  Matrix 

CALL  SEQUENCE 

-  CALL  ELSTIC 

ARGUMENTS 

-  None 

CALLED  BY 

-  GETDK, 

EXTERNALS 

None 

COMMON  BLOCKS 

-  MATAXIS 

FILE  NAMES 

None 

NOTES 

None 

STIFFK ,  TRQDSTR,  UNTFRC 

,  TEMP 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


NOTES 


ELSTIF 

Generate  Element  Mass  and  Stiffness  for  Bars 
CALL  ELSTIF (AE) 

AE  =  AE  for  Bars 


GETDK,  STIFFK 
LMPR0D 

AREA,  EIGN,  L0CAL,  MASS,  N,  STIFF,  TEMP 
None 


The  mass  terms  are  divided  by  386  for  dimen¬ 
sional  consistency  with  the  units  of  the 
material  density  that  are  input. 
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ROUTINE  NAME 

EM0DE 

PURPOSE 

Resize  in  the  Energy  Mode 

CALL  SEQUENCE  - 

CALL  EM0DE 

ARGUMENTS 

None 

CALLED  BY 

-  D0PR0B 

EXTERNALS 

-  LMSIZE,  SQRT 

COMMON  BLOCKS 

-  AMAXMIN,  BASEA,  ELEM,  LAYERD,  MAT,  STRNENG 
TITLE 

FILE  NAMES 

-  TAPE6 

NOTES 

Resize  according  to 

AE  =  AE  ,  j\|u/(RHO*AE  .  ,*L) 

new  old  N  old 

U  =  Strain  Energy 

RHO  =  Density 

L  =  Length 
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ROUTINE  NAME 

-  ERR 

PURPOSE 

Compute  an  Error  Estimate  Between  S  and  T 

CALL  SEQUENCE 

-  CALL  ERR { S ,  T , ERRMAX , L , IMAX) 

ARGUMENTS 

S,T  =  Vectors  to  be  Tested 

ERRMAX  -  Norm  of  Error  Vector 

L  =  Number  of  Vectors  to  be  Tested 

IMAX  =  Index  of  the  L  Vectors  which  has 

the  Largest  Error 

CALLED  BY 

-  DAMAGE,  DMGFREQ 

EXTERNALS 

-  SQRT 

COMMON  BLOCKS 

-  N 

FILE  NAMES 

None 

NOTES 

None 
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ROUTINE  NAME 

ERR0R 

PURPOSE 

Lock  Eigenvectors  that  have  Converged  within 
Tolerance 

CALL  SEQUENCE  - 

CALL  ERR0R 

ARGUMENTS 

None 

CALLED  BY 

-  SIVIB2 

EXTERNALS 

-  SQRT 

COMMON  BLOCKS 

-  EIGN,  EIGVEC ,  N 

FILE  NAMES 

None 

NOTES 

Eigenvector  predictions 

are 

locked  when 

the 

corresponding  error  is 

less 

than  TOLVEC 

(see 

input  instructions)  and 

higher  eigenvector 

have  also  been  locked. 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


NOTES 


F0RSUB 

Forward  Substitution  Routine 
CALL  F0RSUB (V,X, 10 , II) 

V  =  Solution  Vector 
X  =  Right  Hand  Side 

10  =  First  Vector  in  X  to  Use 

11  =  Last  Vector  in  X  to  Use 


ANALYZ ,  DMGITER,  DM0DE,  FRQITER,  SIVIB2 
SK 

N,  SKGM 
None 


None 


ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


NOTES 


FRQITER 

Dynamic  Reanalysis  Iteration 
CALL  FRQITER (IP1) 

IP1  =  Damage  Case  Counter 


DMGFREQ 

BACKSUB,  DMG0RTH,  DX,  F0RSUB,  PREMULT 

DK,  EIGN,  EIGVEC ,  N 

None 


See  Vol.  I  for  a  discussion  of  the  dynamic 
reanalysis  iteration. 
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ROUTINE  NAME  - 
PURPOSE 

CALL  SEQUENCE  - 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS  - 
FILE  NAMES 


NOTES 


GETBC 

Input  Boundary  Conditions 

CALL  GETBC 

None 


INPT 

None 

BND,  N ,  TITLE 
TAPE5,  TAPE 6 


See  input  instructions  for  a  discussion  of 
boundary  conditions . 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


NOTES 


GETBEST 

Returns  to  Previous  Best  Design 

CALL  GETBEST 

None 


D0PR0B 

None 

BASEA,  ELEM,  LAYERD,  SAVE,  TITLE 
TAPES 


In  the  energy  mode,  a  weight  increase  results 
in  the  lower  weight  design  returned.  If  there 
are  displacement  constraints,  GETBEST  returns 
the  lower  weight  design  if  weight  has  more 
than  doubled. 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


NOTES 


GENMSS 

Compute  Generalized  Masses 

CALL  GENMSS 

None 


PRNTMDS 

DX,  GETDK,  PREMULT 

DK,  DMG,  EIGN,  EIGVEC,  GENM,  N 

None 


The  generalized  masses  are  computed  in  both 
the  damaged  and  undamaged  conditions.  The 
units  are  converted  to  pound-force. 


109 


ROUTINE  NAME 

-  GETDK 

PURPOSE 

Set  Up  the  Damaged  Stiffness  and  Mass  Matrices 

CALL  SEQUENCE 

-  CALL  GETDK (M) 

ARGUMENTS 

M  =  Damage  Case  Counter 

CALLED  BY 

-  DAMAGE,  DMGFREQ,  GENMSS 

EXTERNALS 

-  ASEMDK,  C00RD,  ELSTIC,  ELSTIF,  LMPMAS , PLMASS , 
PLSTIF ,  PREPAR,  QDRLTL ,  TREC0N,  TRNSFM 

COMMON  BLOCKS 

-  AREA,  DK,  DMG ,  EIGN,  ELEM,  FRDR,  MASS,  STIFF 

FILE  NAMES 

None 

NOTES 

None 
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ROUTINE  NAME 
PURPOSE 

CALL  SEQUENCE  - 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS  - 
FILE  NAMES 


NOTES 


INC0NN 

Read  Connectivity  Data 

CALL  INC0NN 

None 


INP03 

N0DCHCK ,  PPELEM 

AMAXMIN ,  BASEA,  C0NN ,  ELEM,  TITLE 
TAPES,  TAPE 6 


Read  connectivity  data  and  generate  element 
connectivity  under  certain  circumstances 
{see  input  instructions) .  The  node  sequence 
is  checked  by  calling  N0DCHCK  and  the  element 
connectivity  is  printed. 
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ROUTINE  NAME  - 

INDMGE 

PURPOSE 

Read  Damage  Input 

CALL  SEQUENCE  - 

CALL  INDMGE 

ARGUMENTS 

None 

CALLED  BY 

INPT 

EXTERNALS 

None 

COMMON  BLOCKS  - 

BASEA,  DMG,  DMGERR,  TITLE 

FILE  NAMES 

TAPES,  TAPE 6 

NOTES 

Mass  and  stiffness  damage  factors  are  read 
for  each  element  in  each  damage  case. 
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ROUTINE  NAME  - 
PURPOSE 

CALL  SEQUENCE  - 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


NOTES 


INDSPL 

Read  Displacement  Constraint  Data 

CALL  INDSPL 

None 


INPT 

None 

DEFLMT,  N,  N0DES,  0PT,  TITLE 
TAPE 5,  TAPE 6 


Deflection  limits  for  all  nodes  or  selected 
nodes  are  input  and  printed. 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


NOTES 


INLAYR 

Read  Composite  Layer  Data 
CALL  INLAYR 
None 


INPT 

None 

AMAXMIN,  ANG,  BASEA,  ELEM,  LAYERD,  LAYMIN,  TITLE 
TAPES ,  TAPE 6 


The  fiber  directions  are  read,  initial  layer 
proportions  are  computed,  and  data  is  printed. 
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ROUTINE  NAME  - 
PURPOSE 

CALL  SEQUENCE  - 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS  - 
FILE  NAMES 


NOTES 


INL0ADS 

Read  Loads  Data 
CALL  INL0ADS 
None 


INPT 

None 

DMG,  DRDR,  N.  N0DES,  TITLE 
TAPES,  TAPE 6 


Applied  loads  are  input  and  printed.  Total 
forces  and  moments  are  computed  for  each 
load  condition. 
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ROUTINE  NAME 
PURPOSE 

CALL  SEQUENCE  - 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS  - 
FILE  NAMES 


NOTES 


INITIAL 

Sets  Constants 
CALL  INITIAL 
None 


Main  Program 
ECHO,  0PENSK 

LMTEXCD,  NPR0B,  QUAD,  P0STPR 
TAPES 


The  number  of  problems  is  read,  a  call  to  ECH0 
prints  the  data  listing,  and  the  constants 
NFIL,  MAXSK,  NACTIVE,  MAA ,  MBB ,  MCC  are  set. 
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ROUTINE  NAME 
PURPOSE 

CALL  SEQUENCE  - 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS  - 
FILE  NAMES 


NOTES 


INITP 

Problem  Initialization 

CALL  INITP 

None 


D0PR0B 
INPT ,  P0P 

NPROB,  0PT,  SKGM,  TITLE ,  WEIGHT 
TAPE  6 


This  input  initialization  routine  is  called  at 
the  beginning  of  each  problem.  Calls  are  made 
to  the  routine  which  reads  the  input  and  to  the 
routine  which  maps  the  stiffness  matrix. 
Variables  NPAGE ,  NECYCL,  NDCYCL,  and  WTLAST 
are  set. 
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ROUTINE  NAME 
PURPOSE 

CALL  SEQUENCE  - 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS  - 

FILE  NAMES 


NOTES 


INGNRL 

Read  General  Analysis  and  Optimization  Data 

CALL  INGNRL 

None 


INPT 

None 

AMAXMIN,  ANG,  DEFLMT,  DMG ,  EIGN,  LAYERD,  LAYMIN, 
LMASS ,  N,  0PT,  P0STPR,  STRNENG,  TITLE 

TAPE5,  TAPE 6 


The  general  input  is  read  which  controls  the 
program  options.  This  information  is  also 
printed  with  explanatory  headings. 
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ROUTINE  NAME  - 

INLPMSS 

PURPOSE 

Read  Lumped  Mass  Data 

CALL  SEQUENCE  - 

CALL  INLPMSS 

ARGUMENTS 

None 

CALLED  BY 

INPT 

EXTERNALS 

None 

COMMON  BLOCKS  - 

LMASS,  TITLE 

FILE  NAMES 

TAPES,  TAPE 6 

NOTES 

None 
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ROUTINE  NAME 
PURPOSE 

CALL  SEQUENCE  - 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS  - 
FILE  NAMES 


NOTES 


INP03 

Read  Element  Data 

CALL  INP03 

None 


ELEMIN 

PPHDR,  INC0NN 

BASEA,  MAT,  NMAT,  TITLE 

TAPE 5,  TAPE 6 


Read  materials  properties  and  allowables. 
Call  INC0NN  to  read  the  connectivity. 
Materials  properties  and  allowables  are 
printed. 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 

COMMON  BLOCKS  - 
FILE  NAMES 


NOTES 


INPT 

Calls  all  the  Input  Routines 

CALL  INPT 

None 


INIPT 

INGNRL,  INXYZ ,  ELEMIN ,  INLAYR,  GETBC ,  INL0ADS , 
INDSPL,  INDMGE,  INLPMSS 

BND,  DMG,  FRDR,  LAYERD,  LMASS ,  N,  N0DES ,  0PT 
None 


Various  input  routines  are  called,  depending  on 
the  input  read  in  the  general  input  routine 
INGNRL. 
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ROUTINE  NAME 
PURPOSE 

CALL  SEQUENCE  - 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS  - 
FILE  NAMES 


NOTES 


INXYZ 

Read  Nodal  Data 
CALL  INXYZ 
None 


INPT 

MESHG,  CTYPE ,  PPHDR,  PPNODE 
N0DES,  TITLE 
TAPES,  TAPE 6 


Nodal  data  is  read,  optional  mesh  generator 
called,  coordinates  transformed  to  cartesian 
form  from  cylindrical  or  spherical  if 
necessary  (call  to  CTYPE) ,  and  call  routines 
to  write  header  and  nodal  data  to  the  optional 
post-processor  file. 


122 


ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


NOTES 


LAYCALC 

Calculate  No.  of  Layers  in  Each  Fiber  Direction 
CALL  LAYCALC (L) 

L  =  Element  Number 


LAYPR 


None 

AREA,  BASEA,  ELEM,  LAYERD,  LAYMIN,  LAYPRNT,  SAVE 
None 


None 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 

FILE  NAMES 


NOTES 


LAYPR 

Print  #  Layers  Each  Direction,  Integerize  Design 

CALL  LAYPR 

None 


PRINT 

CURRENT,  LAYCALC 

BASEA,  ELEM,  LAYERD,  LAYMIN,  LAYPRNT ,  OPT,  SAVE, 
TITLE 

TAPE6 


None 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


NOTES 


LLT 

LLT  Decomposition  of  Stiffness  Matrix 

CALL  LLT 

None 


ANALYZ 

PRNTSK,  PUTSK,  SK,  SQRT 
N,  SKGM 
TAPE  6 


The  stiffness  matrix  is  decomposed  by  the 

T 

Choleski  method  into  the  form  K=LL  .  The 
skyline  storage  space  for  K  is  replaced  with 
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ROUTINE  NAME 

-  LMPMAS 

PURPOSE 

Compute  Local  Mass  Matrix  for  Triangles 
and  Quads 

CALL  SEQUENCE 

-  CALL  LMPMAS (L) 

ARGUMENTS 

L  =  Element  Number 

CALLED  BY 

-  GETDK 

,  STIFFK 

EXTERNALS 

None 

COMMON  BLOCKS 

-  AREA, 

C0NN,  ELEM,  L0CAL,  MASS,  TEMP 

FILE  NAMES 

None 

NOTES 

None 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


NOTES 


LMPR0D 

Compute  Mass  Matrix  for  Bar 
CALL  LMPR0D ( AE) 

AE  =  AE  for  Bar 


ELSTIF 

SQRT 

AREA,  L0CAL,  MASS,  N,  TEMP 
None 


None 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


NOTES 


LMSIZE 

Compute  Proportions  of  Fiber  Directions 

CALL  LMSIZE 

None 


DM0DE,  EM0DE 
SQRT 

BASEA,  ELEM,  LAYENG,  LAYERD,  LAYMIN 
None 


The  proportions  are  also  adjusted  to  satisfy 
minimum  size  contraints. 
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ROUTINE  NAME 
PURPOSE 

CALL  SEQUENCE  - 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS  - 
FILE  NAMES 


NOTES 


0PENSK 

Sets  Up  Mass  Storage  for  Stiffness  Array 

CALL  0PENSK 

NONE 


INITIAL 

0PENMS,  WRITMS 

SKBUF 

None 


A  mass  storage  file  is  set  up  to  be  used  as  a 
pseudo  array  for  the  stiffness  array  SK.  This 
array  is  divided  into  MAXREC  records  of  1000 
words  per  record.  Five  records  are  in  core  at 
any  one  time.  A  list  of  the  record  numbers 
for  these  five  records  is  contained  in  the 
INDXSK  array.  The  SK  buffer  SKBUF  contains  the 
last  five  SK  records.  This  is  done  to  minimize 
the  MS  reads  and  writes. 
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ROUTINE  NAME 

0RTH0G 

PURPOSE 

Orthonormal ize  W 

CALL  SEQUENCE  - 

CALL  0RTH0G (W , L0CK , LI , NTRIAL) 

ARGUMENTS 

W  =  Vectors  to  be  Orthonormalized 

L0CK  =  Number  of  Vectors  Locked 

Ll  =  Total  Number  of  Vectors  Less  L0CK 

NTRIAL  =  Number  of  Trial  Vectors 

CALLED  BY 

-  SIVIB2 

EXTERNALS 

-  SQRT 

COMMON  BLOCKS 

-  N 

FILE  NAMES 

-  None 

NOTES 


The  "unlocked”  trial  vectors  are  orthonormalized 
by  the  standard  Gram- Schmidt  process. 
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ROUTINE  NAME 
PURPOSE 

CALL  SEQUENCE  - 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS  - 
FILE  NAME 


NOTES 


MESHG 

Generate  Nodal  Data 

CALL  MESHG (NWORK,MGEN) 

NWQfRK,  work  space  dimension 

MGEN,  generator  option  (positive  number) 


INXYZ 
See  notes 
N0DES ,  see  notes 
See  notes 


This  is  a  user  written  subroutine.  In  the 
current  version  of  ADDRESS,  no  mesh  generator 
is  supplied.  The  routine  may  contain  additional 
common  blocks,  externals,  and  references  to 
files . 
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ROUTINE  NAME 
PURPOSE 

CALL  SEQUENCE  - 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
FILE  NAMES 


NOTES 


N0DCHCK 

Check  Connectivity  Data 

CALL  N0DCHCK ( IELN0 , ITYPE , N0D , D0NE , ERR) 

IELN0  =  element  number 

ITYPE  =  element  type 

N0D(4)  =  element  nodes  for  this  element 
D0NE  =  logical  error  flag  (warning) 

ERR  =  logical  error  flag  (fatal) 


INC0NN 
None 
TAPE  6 


Connectivity  data  checked  to  make  sure 
MA<MB  (bars) 

MA<MB<MC  (triangle) 

MA<min (MB,MC,MD)  (quadrilaterals) 

Data  is  corrected  for  the  bars  and  triangles. 
Error  flag  ERR  is  set  to  true  and  the  program 
stops  if  there  is  an  error  in  the  quadrilateral 
data. 
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ROUTINE  NAME 
PURPOSE 

CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


NOTES 


PLMASS 

Determine  Mass  Matrix  for  Triangular  Elements 
in  the  Local  Coordinates 

CALL  PLMASS {EEMM, AREA) 

EEMM  =  Mass  Matrix  Array 

AREA  =  Area  of  Element 


GETDK,  QDRLTL,  STIFFK 
None 

AREA,  TEMP 
None 


None 
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ROUTINE  NAME 
PURPOSE 

CALL  SEQUENCE 

ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


NOTES 


PLSTIF 

Determine  Stiffness  Matrix  for  Triangular 
Elements  in  the  Local  Coordinates 

CALL  PLSTIF { EEKK , AREA , MA , MB , MC , N0N0RM) 

EEKK  =  Stiffness  Matrix  Array 

AREA  =  Area  of  Element 

MA , MB , MC  =  Node  Numbers  for  Elements 

N0N0RM  =  0  for  membrane  elements 

=  1  for  shear  panel  calculations 


GETDK,  QDRLTL,  STIFFK ,  UNTFRC 
CRAMER 

AREA,  MATAXIS 
None 


None 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


NOTES 


PPDISP 

Write  Displacements  to  Post-processor  File  99 
CALL  PPDISP {IL) 

IL  =  Load  Case  Index 


PRNTDR 

None 

FRDR,  N,  N0DES ,  P0STPR 
NFIL 


None 
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ROUTINE  NAME 

PPELEM 

PURPOSE 

Write  Element 

CALL  SEQUENCE  - 

CALL  PPELEM 

ARGUMENTS 

NONE 

CALLED  BY 

- 

INC0NN 

EXTERNALS 

None 

FILE  NAMES 

NFIL 

NOTES 

None 

Connectivity  to  Post-processor  File 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS  - 
FILE  NAMES 


NOTES 


PPHDR 

Write  Header  to  Post-processor  File 


CALL  PPHDR {ITYPE, 11,12,13, 14, Rl) 


ITYPE=1 

geometry 

=2 

connectivity 

=3 

end  of  data 

=4 

nodal  displacements 

=5 

element  stress 

=6 

element  type 

=7 

element  data 

=8 

end  of  problem  trailer 

11=  number  of  nodes,  ITYPE=1 
=  element  type,  ITYPE=2,6 
=  load  case,  ITYPE=4 
=  number  of  element  types,  ITYPE=5 
=  element  number,  ITYPE=7 

12=  max.  no.  of  nodes/element,  ITYPE=2 
=  number  of  nodes,  ITYPE=4 
=  load  case  number,  ITYPE=5 
*  dimensionality,  ITYPE=6 
=  max.  number  of  nodes/element,  ITYPE=7 

13=  number  of  nodes,  ITYPE=2 
=  interpolation  code,  ITYPE=6 
=  integration  order,  ITYPE=7 

14=  number  of  element,  ITYPE=6 

Rl=  load  parameter,  ITYPE=4,5 


ELEMIN,  INP03,INXYZ, PRINT, PRNTDR 

None 

P0STPR 

NFIL 


None 
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ROUTINE  NAME 

-  PPM0DE 

PURPOSE 

-  Write  Mode  Shapes  to  Post-processor  File 

8 

CALL  SEQUENCE 

-  CALL  PPM0DE 

ARGUMENTS 

None 

CALLED  BY 

-  PRNTMDS 

EXTERNALS 

-  None 

COMMON  BLOCKS 

-  DMG ,  EIGN,  EIGVEC ,  GENM,  N,  N0DES ,  TITLE 

FILE  NAMES 

-  TAPE 8 

NOTES 

Mode  shapes  are  written  to  TAPE 8 .  It  is 
assumed  that  only  the  Z  component  of  the 

vector  is  of  interest  and  that  the  nodes 

numbered  in  accordance  with  the  standard 

are 

conventions  for  a  wing  structure.  Generalized 

masses,  frequencies,  and  (x,y)  locations 
the  grid  points  are  also  output  to  TAPE 8 . 

of 
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ROUTINE  NAME 

PPN0DE 

PURPOSE 

Write  Node  Coodinates  to  Post-processor 
File  99 

CALL  SEQUENCE  - 

CALL  PPN0DE 

ARGUMENTS 

None 

CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


INXYZ 

None 

N0DE,  P0STPR 
NFIL 


NOTES 


All  three  components  of  the  coordinates 
written  to  NFIL. 


are 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS  - 
FILE  NAME 


NOTES 


PPTTL 

Write  Title  to  Post-processor  File 

CALL  PPTTL 

None 


INGNRL 

None 

POSTPR,  TITLE 
NFIL 


The  title  that  is  read  as  input  is  written 
to  a  post-processor  file  NFIL  (set  to  99  in 
INITIAL)  for  possible  later  use. 
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ROUTINE  NAME 
PURPOSE 

CALL  SEQUENCE  - 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS  - 
FILE  NAMES 


NOTES 


P0P 

Map  the  Stiffness  Matrix 

CALL  P0P 

None 


INITP 


None 

BASEA,  C0NN ,  ELEM,  N,  SKGM,  TITLE 
TAPE  6 


The  gross  population  (total  number  of  elements 
in  the  upper  triangle)  and  the  apparent 
population  (total  number  of  elements  in  the 
skyline)  of  the  stiffness  matrix  are  computed 
for  the  unrestrained  structure.  The  apparent 
population  must  be  not  greater  than  MAXSK  which 
is  defined  in  INITIAL. 
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ROUTINE  NAME 

-  PRELT 

PURPOSE 

-  Compute  X  =  LT*Y 

CALL  SEQUENCE 

-  CALL  PRELT {X, Y,  N) 

ARGUMENTS 

X  =  Vectors  to  be  Computed 

Y  =  Given  Vectors 

N  =  Number  of  Vectors  to  be  Computed 

CALLED  BY 

-  DMG0RTH 

EXTERNALS 

-  SK 

COMMON  BLOCKS 

-  N,  SKGM 

FILE  NAMES 

None 

NOTES 


The  L  matrix  is  stored  in  skyline  form. 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


NOTES 


PREMULT 

Matrix  Vector  Multiplication:  Y  =  GM*Z 
CALL  PREMULT (Y, Z,N,L) 

Y  =  Vectors  to  be  Computed 
Z  =  Given  Vectors 

N  =  Number  of  Vectors  to  be  Computed 
L  =  Shift  Factor  for  Y 


FRQITER,  GENMSS ,  SIVIB2 
None 

GM,  N,  SKGM 
None 


The  mass  matrix  is  stored  in  skyline  form. 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 

FILE  NAMES 


NOTES 


PREPAR 

Compute  Moduli  and  Allowable  Stresses 
CALL  PREPAR (L, BA, BAE, INDEX, LDS) 

L  =  Element  Number 

BA  =  Reference  Thickness  Used  to  Establish 
Triangle  and  Quadrilateral  Element 
Thicknesses  and  Moduli 

BAE  =  Reference  Thickness  Used  to  Compute 
Allowable  Stresses 

INDEX  =  0,  Compute  Element  Thicknesses 
=  1,  Compute  Allowable  Stresses 

LDS  =  #  of  Load  Cases 


GETDK,  PRNTEL,  STIFFK ,  STRC0N,  UNTFRC 
None 

ALSTRS,  AREA,  EIGN,  ELEM,  ENERGY,  ESTRESS ,  LAYERD, 
MASS,  MAT,  MATAXIS,  NMAT ,  STIFF,  TEMP 

None 


BA  =  1.0  when  PREPAR  is  called  by  STIFFK, 
GETDK,  STRC0N,  and  UNTFRC 

=  BASEA  when  PREPAR  is  called  by  PRNTEL 

BAE  =  1.0  when  PREPAR  is  called  by  STIFFK, 
GETDK,  PRNTEL,  AND  UNTFRC 

=  BASEA  when  PREPAR  is  called  by  STRC0N 

INDEX  =  0  when  PREPAR  is  called  by  STIFFK,  GETDK, 
and  UNTFRC 

“  1  when  PREPAR  is  called  by  PRNTEL, 

STRC0N ,  and  UNTFRC 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


NOTES 


PRINT 

Scale  Thicknesses,  Call  Various  Print  Routines 

CALL  PRINT 

None 


D0PR0B 

LAYPR ,  PPHDR ,  PRNTDR ,  PRNTEL ,  PRNTMDS 
BASEA,  EIGN,  ELEM,  LAYERD,  SAVE 
None 


None 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


PRNTDR 

Print  Table  of  Node  Information 

CALL  PRNTDR 

None 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


PRINT 

PPDISP,  PPHDR 

BND,  FRDR,  N,  N0DES ,  TITLE 
TAPE  6 


NOTES 


Coordinates,  applied  loads,  and  displacements 
are  printed. 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


PRNTEL 

Print  Element  Information 

CALL  PRNTEL 

None 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 

FILE  NAMES 


PRINT 

C00RD,  ELF0RC,  PREPAR,  SK,  TRQD0UT,  TRQDSTR 

ALSTRS,  AREA,  BASEA,  C0NN,  DMG ,  ELEM,  ENERGY, 
ESTRESS,  FRDR,  LAYERD,  SAVE,  STRESS,  TEMP,  TITLE 

TAPE6 


NOTES 


Element  thicknesses,  connectivity,  stresses,  and 
strain  energies  are  printed. 
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ROUTINE  NAME 

-  PRNTSK 

PURPOSE 

Print  the  Structural  Stiffness  or  Mass  Matrix 
by  Rows 

CALL  SEQUENCE 

-  CALL  PRNTSK (SK) 

ARGUMENTS 

-  SK  =  Pseudo  Stiffness  Matrix 

CALLED  BY 

-  LLT 

EXTERNALS 

-  SK 

COMMON  BLOCKS 

-  N,  SKGM,  TITLE 

FILE  NAMES 

-  TAPE6 

NOTES 

-  This  replaces  the  routine  PRINTK  in  previous 

versions  of  this  program. 
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ROUTIN  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


PRNTMDS 

Print  Out  Frequencies  and  Modeshapes 

CALL  PRNTMDS 

None 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


PRINT 

ATAN,  GENMSS,  PPM0DE,  REST0R,  SQRT 

DMG,  EIGN,  EIGVEC ,  GENM,  N,  N0DES ,  P0STPR,  TITLE 

TAPE  6 


NOTES  -  Before  mode  shapes  are  printed  they  are 

normalized  so  that  the  maximum  component  is 

+  1.0. 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


NOTES 


PUTSK 

Buffer  Bookkeeping 
CALL  PUTSK (IW0RD, VALUE) 
See  Notes 


ASEMBL,  B0UNDS,  LLT,  STIFFK 
SK 

SKBUF 

None 


This  routine  loads  the  value  "VALUE"  into  the 
pseudo  array  SK;  i.e.  this  routine  replaces 
the  statement  SK(IW0RD)  =  VALUE. 
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ROUTINE  NAME 

QDRLTL 

PURPOSE 

Compute  Stiffness  and  Mass  Matrices  for  the 
Quadrilateral  and  Shear  Panel  Elements 

CALL  SEQUENCE  - 

CALL  QDRLTL (AREA, L,N0) 

ARGUMENTS 

AREA  =  Element  Area 

L  -  Element  Number 

N0  =  Control  Parameter  (See 

C0NDNS  Arguments) 

CALLED  BY 

EXTERNALS 

COMMON  BLOCKS 

FILE  NAMES 

-  GETDK,  STIFFK,  TRQDSTR,  UNTFRC 

-  C0NDNS,  PLMASS ,  PLSTIF,  SUM 

-  AREA,  C0NN,  EIGN,  ELEM,  MASS,  QUAD,  STIFF 

-  None 

NOTES 

This  routine  calls  the  triangular  element 
routines  for  each  of  the  four  elements  making 
up  the  quadrilateral  and  makes  sure  that  the 

results  are  properly  combined  into  one  mass 

and  one  stiffness  matrix. 
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ROUTINE  NAME 
PURPOSE 

CALL  SEQUENCE 
ARGUMENTS 


QLSTRS 

Compute  Stresses  for  4  Triangles  of  the 
Quadrilateral 

CALL  QLSTRS (L) 

L  =  Element  Number 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 

FILE  NAMES 


TRQDSTR 

STRESS 

ALSTRS,  ELEM,  ENERGY,  ESTRESS ,  FRDR,  QUAD,  STIFF, 
STRESS 

None 


NOTES 


None 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


NOTES 


RANDOM 

Enter  Random  Vectors  into  Certain  Columns  of  W 
CALL  RAND0M (W, IFLAG) 

W  =  Initial  Array  of  Eigenvectors 
IFLAG  =  See  Notes 


SIVIB2 
None 
EIGN,  N 
None 


If  IFLAG=10,  random  numbers  are  entered  in 
just  the  last  column  of  W.  When  IFLAG=0 , 
random  members  are  generated  for  all  columns 
of  W. 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


NOTES 


REDUCE 

Eliminate  Rows  from  Vector  for  Boundary  D0F 
CALL  REDUCE (F,L) 

F  =  Vector  to  be  Reduced 
L  =  Number  of  Vectors 


ANALY2 ,  DM0DE,  DX 

None 

BND ,  N 

None 


None 
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ROUTINE  NAME 
PURPOSE 

CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


NOTES 


REST0R 

Restore  the  Displacement  or  Force  Matrix  to  Full 
Size 

CALL  REST0R (D, L) 

D  =  Vector  to  be  Restored 
L  =  Number  of  Vectors 


DMIDE,  DX,  PRNTMDS ,  SCALE 

None 

BND ,  N 

None 


None 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


SAVE 

Save  Current  Least-weight  Design 

CALL  SAVE 

None 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


D0PR0B 

None 

BASEA,  ELEM,  LAYERD,  SAVE,  WEIGHT 
None 


NOTES 


The  design  is  saved  if  the  weight  decreased 
from  the  previous  cycle. 


ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 

FILE  NAMES 


NOTES 


SCALE 

Adjust  Scaling  Factor  to  Satisfy  Constraints 

CALL  SCALE 

None 


CURRENT 

None 

AMAXMIN,  BASEA,  DEFLMT,  ELEM,  FRDR,  LMTEXCD , N , 
0PT,  STRNE 

None 


Maximum-size  constraints  are  also  impared  by 
this  subroutine. 
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ROUTINE  NAME 

SIVIB2 

PURPOSE 

Compute  Eigenvalues  and  Eigenvectors 

CALL  SEQUENCE  - 

CALL  SIVIB2 

ARGUMENTS 

None 

CALLED  BY 

-  CURRENT 

EXTERNALS 

-  BACKSUB,  DEC0UP,  ERR0R,  F0RSUB,  0RTH0G,  PREMULT 
RAND0M ,  SEC0ND 

COMMON  BLOCKS 

-  EIGN,  EIGVEC ,  N 

FILE  NAMES 

-  TAPE 6 

NOTES 

This  is  the  main  calling  routine  for  the 
eigensolution.  See  Vol.  I  for  a  discussion 
of  the  simultaneous  iteration  techniques  used 

in  this  routine. 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


SK 

Bookkeeping  for  Stiffness  Matrix,  SK 
FUNCTI0N  SK(IW0RD) 

IW0RD  =  See  Notes 


CALLED  BY 

-  ASEMBL, 
PRNTEL, 

BACKSUB , 
STRC0N, 

B0UNDS,  F0RSUB ,  LLT ,  PRELT, 
PUTSK,  PRNTSK 

EXTERNALS 

-  READMS , 

WRITMS 

COMMON  BLOCKS 

-  SKBUF 

FILE  NAMES 

-  TAPES 

NOTES  -  This  function  returns  the  value  of  SK(IW0RD) 

from  the  pseudo  array  SK.  Since  SK  has  been 
replaced  by  a  MS  file,  this  routine  performs 
all  the  bookkeeping  required  to  access  the 
correct  word  on  the  correct  record. 
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ROUTINE  NAME 
PURPOSE 

CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 

COMMON  BLOCKS 

FILE  NAMES 

NOTES 


STIFFK 

Set  Up  Total  Stiffness  and  Mas  Matrices  for 
Structure 

CALL  STIFFK 

None 


ANALYZ 

ASEMBL,  C00RD ,  ELSTIC,  ELSTIF,  LMPMAS ,  PLMASS , 
PLSTIF,  PREPAR,  PUTSK,  QDRLTL,  TREC0N ,  TRNSFM 

AREA,  BASEA,  EIGN,  ELEM,  FRDR,  GM,  LMASS ,  MASS, 
N,  SKGM,  STIFF,  STRNENG 

None 


This  is  the  main  calling  routine  which  sets 
up  the  global  mass  and  stiffness  matrices. 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUNCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 

FILE  NAMES 


NOTES 


STRC0N 

Determine  Scale  Factor  Needed  to  Satisfy  Constraints 

CALL  STRC0N 

None 


SCALE 

C00RD,  ELF0RC ,  PREPAR,  SK,  TRQDSTR 

ALSTRS ,  AREA,  BASEA,  DMG,  ELEM,  ENERGY,  ESTRESS, 
FRDR,  STRESS,  STRNENG,  TEMP 

TAPES 


None 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


STRESS 

Compute  Strains  and  Stresses  for  Triangular  Element 

CALL  STRESS (UV,MA, MB ,MC, ENG, NND) 

UV  =  Local  Strains 

MA, MB ,MC  =  Node  Numbers 

ENG  =  Element  Strain  Energy 

NND  =  5  If  Triangular  Element  is  Part  of  a 
Shear  Panel 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


QLSTRS,  TRQDSTR 
CRAMER,  SQRT 

ALSTRS,  AREA,  ESTRESS ,  FRDR,  MATAXIS,  STRESS 
None 


NOTES 


None 
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ROUTINE  NAME 

-  SUM 

PURPOSE 

“  Assemble  a  Single  Matrix  from  4  Triangular 
Matrices 

CALL  SEQUENCE 

-  CALL  SUM (EKM, EEKM,MA, MB , MC) 

ARGUMENTS 

-  EKM  =  Mass  or  Stiffness  Matrix  to  be  Computed 

EE  KM  =*  Mass  or  Stiffness  Matrices  for  the 
Individual  Triangular  Elements 

MA , MB , MC  =  Node  Numbers  for  the  Triangular 
Element 

CALLED  BY 

-  QDRLTL 

EXTERNALS 

None 

COMMON  BLOCKS 

None 

FILE  NAMES 

None 

NOTES 

None 
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ROUTINE  NAME 

TREC0N 

PURPOSE 

Coordinate  Transformations:  Orthotropic 
Materials 

CALL  SEQUENCE  - 

CALL  TREC0N (L, IND} 

ARGUMENTS 

L  =  Element  Number 

IND  =  1  for  0°  calculations 

2  for  90°  calculations 

3  for  +45°  calculations 

4  for  -45°  calculations 

CALLED  BY 

GETDK,  STIFFK,  TRQDSTR,  UNTFRC 

EXTERNALS 

C0S ,  SIN,  SQRT 

COMMONG  BLOCKS  - 

■  ANG,  L0CAL,  MATAXIS 

FILE  NAMES 

None 

Notes 

None 
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ROUTINE  NAME 
PURPOSE 

CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 
FILE  NAMES 


NOTES 


TRNSFM 

Transform  a  Matrix  From  Local  to  Global 
Coordinates 

CALL  TRNSFM (EKM,CCC,NN0DES) 

EKM  =  Element  Matrix  in  Local  Coordinate 

CCC  =  Global  Matrix 

NN0DES  =  3  For  Triangular  Element 
=  4  For  Quad  Element 


GETDK ,  STIFFK 
None 

L0CAL,  N 
None 


None 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 

FILE  NAMES 


NOTES 


TRQD0UT 

Output — Triangles  and  Quadrilaterals 
CALL  TRQDOUT (L) 

L  =  Element  Number 


PRNTEL 

None 

ALSTRS,  AREA,  C0NN,  ELEM,  ENERGY,  ESTRESS, 
FRDR,  LAYERD,  S 

None 


Print  element  thickness,  stresses,  and  energy. 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 
COMMON  BLOCKS 

FILE  NAMES 


NOTES 


TRQDSTR 

Compute  Triangle  and  Quadrilateral  Stress,  Energy 
CALL  TRQDSTR (L, BA) 

L  =  Element  Number 

BA  =  1.0  for  stresses  in  element  of  unit 
thickness 

~  BASEA  for  stresses  in  current  design 


PRNTEL,  STRC0N 

ELSTIC,  QDRLTL,  QLSTRS ,  STRESS,  TREC0N 

AREA,  DMG,  ELEM,  ENERGY,  ESTRESS,  FRDR,  LAYENG, 
LAYMIN,  ST 

None 


Stresses  and  engergy  for  both  isotropic  and 
composite  elements  are  computed. 
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ROUTINE  NAME 

-  UNI TEG 

PURPOSE 

-  Determine  Energy  Due  to 

Dummy 

Unit  Loads 

CALL  SEQUENCE 

-  CALL 

UN I TEG (L, II, INDX) 

ARGUMENTS 

L  = 

Element 

Number 

II  = 

1,  0° 

composite 

larger 

calculation 

2,  90° 

composite 

larger 

calculation 

3,  +45° 

composite 

larger 

calculation 

4,  -45° 

composite 

larger 

calculation 

INDX 

=  0  if 

this  is  not 

a  composite  element 

CALLED  BY 

-  UNTFRC 

EXTERNALS 

None 

COMMON  BLOCKS 

-  AREA,  DMG,  DUMMY,  ELEM,  FRDR,  LAYENG,  LMTEXCD 

STIFF,  STRNENG 

FILE  NAMES 

-  None 

NOTES 

None 
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ROUTINE  NAME 
PURPOSE 
CALL  SEQUENCE 
ARGUMENTS 


CALLED  BY 
EXTERNALS 

COMMON  BLOCKS 

FILE  NAMES 


NOTES 


UNTFRC 

Determine  Stress  Due  to  Dummy  Unit  Loads 

CALL  UNTFRC 

None 


DM0DE 

C00RD,  ELF0RC,  ELSTIC,  PLSTIF,  PREPAR,  QDRLTL, 
TREC0N,  UNITEG 

ANG,  AREA,  BASEA,  DUMMY,  ELEM,  FRDR,  LMTEXCD , 
STIFF 


None 


None 
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ROUTINE  NAME 

-  WEIGHTS 

PURPOSE 

-  Determine  Weights  for  Structure,  Scale  Displacements 

CALL  SEQUENCE 

-  CALL  WEIGHTS 

ARGUMENTS 

None 

CALLED  BY 

EXTERNALS 

-  CURRENT 

None 

COMMON  BLOCKS 

FILE  NAMES 

-  BASEA,  EMG,  ELEM,  FRDR, 
0PT,  WEIGHT 

-  TAPE 6 

LMASS ,  : 

MAT,  N,  NPR0B, 

NOTES 

Total  structural  weight 

as  well 

as  weights 

of  all  the  elements  by 
printed. 

type  are 

computed  and 
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APPENDIX  D 

UPDATES  FOR  OUT  OF  CORE  SOLUTIONS 
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*TD  DELETE  1 
*D  RFTr.iAMDPT.52 

COMMON  /DM/  GM ( 1 ) 

*D  RFTDAM0PT.236 

COMMON  /GM/  GM ( 1 ) 

*D  RFTDAMOF'T.  2946 

COMMON  /GM/  GM(1) 
fl-’j  PF  1  D AMO PT.  3707 

COMMON  /GM/  GM ( 1 ) 

#.D  RFTDAMOF'T .  555 

i  OMMON  /EIGVEC/ 

*D  RFTDAMOF'T.  646 

COMMON  /EIGVEC/ 

*D  RFTDAMOF'T.  718 

COMMON  /EIGVEC/ 

RF rDAMOP T . 799 

COMMON  /EIGVEC/ 

*D  RF TDAMOPT.  1247 

I  OMMON  /EIGVEC/ 
s  D  R'P  TDAMOPT.  1305 

i' OMMON  /EIGVEC/ 

* D  RFTPAM OP T . 1 340 

i  OMMON  /EIGVEC/ 

<U  RE  1  DA MOP T .  1 887 
5i 40  FORMAT ( 3F10. 3) 

*D  RFTDAMOET. 2871 

OMMON  /EIGVEC/ 

T’t  RF  TDAMOPT.  3262 

r OMMON  /EIGVEC/ 

*D  RFT D AMORT. 3633 

COMMON  /EIGVEC/ 

*ID  R ( GSK 
•tf-D  RF  TDAMOPT.  3 

S'  T  AP  E.  99 = DP  OS  T ,  TAPE4  ) 

*■1  RFTDAMOF'T.  1836 
CALL  0 PENCK 
RFID AM OPT. 55 
RFT  DAMOPT  .21? 

WB  RFT DAMuPT. 238 
v  Cl  RE  1 1 1  AMORT.  1281 
frD  RE TDAMOPT. 2338 
an  RF  IDA  MART.  2927 
T.i  RF*  DAMOPT.  371 1 
:>  U  RF  1  DAMOPT.  J  06 

20  CAL L  PUTSK <JX, SK ( J X  >  +C ( I , J ) ) 

*D  RE TDAMOPT. 272 

CALL  P U T SK ( K  X , SK (KXX) ) 

*D  RFTDAMOF'T.  2343 

CALL  PUTSK ( 1 , SORT ( SK ( 1 ) ) ) 

*D  RF TDAMOPT. 2357 

4  CALL  PUTSK ( J+ 1 Q , EL /SK ( I D 1  AG  C  J ) > > 

*1  RFTDAMOF'T .  2339 
EXTERNAL  SK 


U  <  350  350  ,  6  )  ,  W  (  1 ,  1 ,  1  )  ,  XL  (  1 0 , 3  ) ,  V  (1  ,  1  ) 

U ( 350 , 6 ) ,  V  ( 350 , 6 ) , W ( 1 , 1 , 1 ) , XL  < 1 0 , 3  > , Y ( 1 , 1 ) 
U < 350 7  6  > , V 1 350, 6 ) , W  tl , 1 , 1 ) , XL ( 1 0 ,  3 ) , V  C 1 , 1 ) 
U ( 350 . 6 ) , V ( 350 , 6 ) 7  W ( 1 , 1 , 1 ) ,  XL ( 1 0*  3  J , V ( 1 , 1 ) 
U ( 350  7  6 ) , V ( 350 , 6 ) , W  < 1 , 1 , 1 ) , XL ( 1 0 , 3  > , Y  < 1 , 1 ) 
U ( 350, 6 ) 7  V  f 350,6 ) >  W ( 1 , 1 , 1 ) , XL< 10, 3 ) , Y 1 1 , 1 > 
U  <  350 , 6  > ,  V  ( 350 ,  6 ) ,  W  ( 1 ,  1 ,  1  > ,  XL  (  10 , 3 ) ,  Y  <T  ,  1  ) 

U ( 350 , 6 ) , V ( 350 , 6 ) >  W ( 1 , 1 , 1 > , XL ( 1 0 , 3) , Y ( 1 ,  1 > 
U ( 350,6 ) , V< 350, 6 ) , W ( 1 > 1 , 1 ) , XL ( 1 0, 3 ) , Y ( 1 , 1) 
U ( 350,6 ) , V ( 350, 6  > , W ( 1 ,  1 , 1 ) , XL ( 10, 3 ) , Y ( 1 ,  1  ) 
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a-L.i  RF  TDAMOF'T .  236'? 

i.  All  PRNTSMSK) 

■+1  D  KF  T PAMuF'T  .  37 20 
10  GALL  PUTSK(I,0. ) 

* r.i  RFTDAMOPT.  236 0 

CALL  PUT SK ( I DI AG k  I  )  » SQRT <  EL )  ) 
« 1  RFTDAMOPT. 4433 

SUBROUTINE  GPENSK 


r 

r 


i 


Tuts  ROUTINE  SETS  UP  A  MASS  STORAGE  FILE  TO  BE  USED 
FOR  T HI  STIFFNESS  ARRAY  SK.  SK  WILL  BE  DIVIDED  INTO 
OF  1000  WORDS  PER  RECORD.  5  RECORDS  WILL  BE  IN  CORE 
A  l  1ST  OF  THE  RECORD  NUMBERS  FOR  THESE  5  RECORDS  IS 
THE  ARRAY  INDXSK.  the  MAPPING  OF  SK  ONTO  THE  MASS  ST 


AS  A  PSEUDO  ARRAY 
MAX REG  RECORDS 
AT  ANY  ONE  TIME. 
CONTAINED  IN 
ORAGE  FILE  IS 


AS  FOLLOWS 


■  '-4  (  i  )  THROUGH  SK (  1000 )  RESIDES  ON  MS  RECORD  1 

L  SKtiOOl)  -  SK  <  2000 )  RESIDES  ON  MS  RECORD  2 

i 

C  AND  SO  ON 

t 


0  THE  SK  BUFFER  SKBUF  CONTAINS  THE  LAST  5  SK  RECORDS  THAT  HAVE  BEEN 

r  ACCESSED.  THIS  IS  DONE  SO  THE  AT  ,  HOP  FULLY  ,  THE  NUMBER  OF  MS  READS 

i  AND  WRITES  WILL  BE  MINIMIZED. 

J.  JENSEN  11/1/81 

C ,o MHO N  / SK BU F /  SK BUF  (1000,5),  I NDXSK  ( 5 )  >  MAXREC ,  I D I RSK  (  101  ) ,  NXTRF.C 

r., 

C 

MAX RFC  =  100 

CAL L  OF’ENMS  (  4 ,  I D I  RSK ,  MAXREC+ 1,0) 

C 

c  INITIALIZE  THE  MASS  STORAGE  FILE  TO  DEBUG  VALUES 

C 

DO  10  I  a  1, MAXREC 

1 0  CALL  WR I TMS ( 4 , SKBUF ( 1 , 1 ) , 1 000 ,1,-1) 

C 

C  SET  THE  BUFFER  INDEX  TO  ALL  ZEROES  (  NO  RECORDS  IN  THE  SK  BUFFER) 

f 

D'  i  20  I  =  1,5 
20  I NDXSK ( I )  =  0 

NXTREC  =  1 

C 

C  NXTREL.  POINTS  TO  THE  OLDEST  RECORD  IN  THE  BUFFER.  THIS  RECORD  WILL 
C  REPLACED  WHEN  IT  IS  NECESSARY  TO  LOAD  A  NEW  RECORD  INTO  THE  BUFFER. 

0 

RETURN 

END 

FUNCTION  SK(IWORD) 

C 
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i  THIS  FUNCTION  RETURNS  THE  VALUE  OF  SK < 1  WORD )  FROM  THE  PSEUDO  ARRAY  SK. 
i  SINCE  SP  HAS  BEEN  REPLACED  BY  A  MS  FILE,  THIS  ROUTINE  PERFORMS  ALL. 

'1  THE  POOF  I  EEP.INO  REQUIRES)  TO  ACCESS  THE  CORRECT  WORD  ON  THE  CORRECT 
c  Rl  T  URD 

i 

COMMON  /SK BUR /  SKBUF < 1000,5) , INDXSK <5> , MAXREC, IDIRSK (101), NXTREC 

i  i.OMpmfc  1  HE  RECORD  NUMBER  ON  WHICH  SK  <  1  WORD )  RESIDES 
c 

5  ITS  i  "  (i  WORD  -1)  /  1000  +  1 

f 

lHBTt  To  SEE  OF  THE  RECORD  IS  ALREADY  IN  THE  BUFFER 

L 

Mil  1  l.)  I  —  1,5 

IF*.]  N  B  XSK  < I )  .  NE . I REC  >  GO  T  0  1 0 

'  RE i  UKD  IF  IN  CURE,  RETURN  THE  APPROPRIATE  WORD 

1W  •-  IWORD  -  (  I  REC  -  1)  *  1000 
SE  ■=  SE  BIJF  (  1  W ,  I  ) 

Kp TURN 

1 0  t  i  in  i  I  Ml  IE 

i 

i  REi  ORD  I NO  1  IN  THE  BUFFER  SO  RETRIEVE  IT 

IMIF'K.  U  -  MAXREC  )  GO  TO  1 5 

WRITE*/-.,  11)  IWORD,  IREC, MAXREC 
1  FORMAT ( "  IN  SK  I  WORD, IREC, MAXREC=U , 31 14 ) 

STOP  "IREC  TO  BIG  IN  SK,  ERROR" 

1 5  CONTINUE 

i  j 

i  FIRST  WRITE  THE  RECORD  CURRENTLY  IN  THE  BUFFER  BACK  TO  DISK 

C 

1H1  N 1 1 X SK  ( N X TREC )  .  NE .  0 )  CALL  WR I TMS  (  4 ,  SKBUF  <  1 ,  N X TREC  ) ,  1 000 , 

X  INDXSK ( NXTREC) , -1 ) 

i  AND  READ  IN  THE  REQUIRED  RECORD 

C 

CALL  RE AD MS ( 4 , SKBUF  *  1 , NXTREC ) , 1 000 ,  I REC ) 

C 

i  AND  PUT  THE  RECORD  NUMBER  IN  THE  BUFFER  DIRECTORY 

C 

INDXSK  <  NXTREC )  =  IREC 
NXTREC  =  NXTREC  +  i 
LF (NXTREC. GT. 5)  NXTREC  =  1 

i  NOW  THAT  THE  RECORD  IS  IN  THE  BUFFER,  Gu  BACK  AND  RETRIEVE  THE 
i-  REQUESTED  WORK 

GO  TO  5 
END 

•  IJBRO U TINE  F'UT SK (  1 W 0 R D ,  V AL UE  ) 
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o  o 


C  THIS  SUBROUTINE  LOADS  THE  VALUE  "VALUE "  INTO  THE  PSEUDO  ARRAY  SK. 
THAT  IS  THIS  ROUTINE  REPLACES  THE  STATEMENT 
SKI  I WORD)  =  VALUE 

0 

c 

c OMMON  /SKBUF /  SKBUF ( 1 OOO , 5 ) , INDXSK 1 5  > , MAXREC , I D I RSK (101), NXTREC 

C 

C  F  IR'-.T  MAKE  SURE  SKI  I  WORD)  IS  IN  THE  BUFFER  BY  READING  SKI  I  WORD) 


i 

C 

r 


DUMMY  =  SKI  I WORD) 

F  I ND  WH 1 1 H  BUFFER  SK ( I WORD ) 


IS  IN 


1REL  =  (IWORD  -  1)  /  1000  +  1 
DO  10  I  =  1,5 

IF t INDXSK ( I ).NE. IREC)  GO  TO  10 
IW  =  IWORD  -  (IREC  -  1)  *  1000 
SKBUF (I W, I )  =  VALUE 
RETURN 

1"  CONTINUE 


L  J  t 

C. 


c 

c 

c 


NiU  FOUND,  AN  ERROR  HAS  OCCURED  SOMEWHERE 

STOP  "ERROR  IN  PUTSK,  RECORD  NO  IN  BUFFER" 

END 

SUBROUT  I NE  PRNTSK ( SK ) 

COMMON  /N/  MM , NM, NN 

■  OMMON  ,/SKGM/  I  COL  ( 270 ) ,  I D I  AG  1 270 )  ,  NONZRO 
COMMON  /TITLE/  NPAGE, TITLE ( 8 ) 

FXTERNAL  SK 

PRINT  THE  STRUCTURAL  STIFFNESS  OR  MASS  MATRIX 
BY  ROWS 


WR I TE ( 6 , 3000 >  TITLE, NPAGE 
NPAGE-NPAGE+ 1 
WRITE <6, 1000)  1 
WR [ TE ( 6* 2000 )  SK(l) 
un  SO  1=2, NM 

V X  =  I D I AG ( 1-1 )  +  l 
K Y-IDIAGI I ) 

WRITE  1 6, 1000)  I 
WR 1 TE I 6, 2000 )  I SK I K ) , K=KX , KY ) 
CONTINUE 

t *■'» K'  FORMAT  (2 X,  IS) 

20U0  F ORMAT ( 1 0  X , i OE 1 2 . 4 ) 

3000  FORMAT  1 1  HI /30X , 8A10, 1  OX , 5HPAGE  ,I3//30X, 
%  *TOTAL  STIFFNESS  OR  MASS  MATRIX*) 
RETURN 
END 


*D  RF  T  DA MOP T . 839 


COMMON  /MAT/  ALSTRSI 2265 ) , ELCNST 1 90S ) , POISON ( 453 ) , RH01 I  453 ) , 
X  YOUNGMI 453 ) 
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*D  RFTDAMOPT. 1 190 

COMMON  /MAT /  ALSTRS ( 2265 ) , ELCNST ( 906 ) , PO I SON ( 453 ) , RHO 1 ( 453 ) , 
X  YQUNGM ( 453) 

*D  RFTDAMOPT. 2115 

COMMON  /MAT /  ALSTRS ( 2265 ) ,  ELCNST < 906 )  ,  PO I SON ( 453 )  ,  RHO 1 ( 453 )  , 
X  YOUNGM ( 453 ) 

*D  RFTDAMOPT. 2981 

COMMON  /MAT /  ALSTRS <  2265 ) , ELCNST ( 906 ) , POISON < 453 ) , RHO 1 ( 453 ) , 
X  YOUNGM (453) 

*D  RFTDAMOPT. 4373 

COMMON  /MAT/  ALSTRS  <  2265 ) , ELC  NST ( 90 6 ) , PO I SON ( 453  > , RHO 1 ( 453 ) , 
X  YOUNGM (453) 

*D  RFTDAMOPT. 1918 

GO  T 0  (90,40, 50 , 60, 70, 80, 85 ) ,  L A M < I )  + 1 

*  INSERT  RFTDAMOPT. 1936 

GO  10  90 
85  AEX ( I ) =. 50 
AEY ( I >=. 10 
#U  RFTDAMOPT. 1835 

MAXSK  =  100000 
■*  I  RFTDAMOPT.  2000 

IF(KY.LE.O)  GO  TO  370 
*D  RFTDAMOPT. 214 

REAL  U( 270, 1 1 ) , W ( 270, 1 1 ) 

*D  RFTDAMOPT. 1283 

REAL  0(270, 1 1 ) , X (270, 1 1 ) 

*  D  RFTDAMOPT. 38 40,384 2 

SX (K)=E1*(EDRC2,K)-EDR< 1 ,K)  >/AL 

IF ( SX ( K ) . GE. 0. )  ESRTIO ( K ) =SK ( K ) /ALS ( 1 ) 

IF ( SX (K ) . LT. 0 . )  ESRT 1 0 ( K ) = ABS ( S  X ( K ) )/ALS(2) 

*D  RFTDAMOPT. 3225, 3227 

SX ( K ) =E i * ( EDR ( 2, K ) -EDR ( 1 , K) ) /AL 

IF ( SX ( K ) . GE. 0. )  ESRTIO ( K ) =SK ( K ) /ALS ( 1 ) 

IF ( SX ( K ) . LT. 0. )  ESRTIO (K ) =ABS ( SX ( K ) )/ALS(2) 
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APPENDIX  E 

TRUSS  PROGRAM  LISTING 
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PROGRAM  TRUSS 

COMMON  /ORDER/  N, NOPTION 

COMMON  / YML V /  E , F  <  S  > , OLDF ( 8 > , X  < 1 0 ) ,  X B ( 1 0 ) 

CHARACTER* 1  ANSWER 

C 

c - — — - 

c 

C  THIS  PROGRAM  DOES  A  DEFLECTION  AND  STRESS  ANALYSIS  ON 
C  A  "SIMPLE'*  TRUSS. 

c 

c  — ■ — ■ — - — — ■- — - - ■ — - - - ■ - ■ - 

c 

C  PR  I IMC I  PAL  VAR  I ABLES  = 


i" 

r 

C 

c 

I 


r 


F  ( 3 )  — — > 

X(  10)  — 

XB(iO>  - > 

EK (Hi  y  i  1  (.0  j? 

SK  (  S  i  S ) - > 

S(10,8)  - > 

R  (  8  )  y 

ST(8)  :> 


YOUNG'S  MODULUS 
APPLIED  LOAD  VECTOR 
UNMODIFIED  MEMBER  SIZES 
MODIFIED  MEMBER  SIZES 
ELEMENT  MATRICES 
SYSTEM  STIFFNESS  MATRIX 
STRESS  MATRIX 

SOLUTION  FOR  (UN) MODIFIED  STRUCTURE 
ELEMENT  STRESSES  FOR  STRUCTURE 


C  DATE!  MARCH  18,  1982  KEITH  MILLER,  PHONE: 229-4235 

C 

C  - — - 


C  OPEN  INPUT  &  OUTPUT  FILES  FOR  INTERACTIVE  PURPOSES. 

C 


c 

c 

C: 


C 

c 

c 


OPEN (5,  F ILE  =  '  INPUT-' ) 
OPEN  ( 6 ,  F I  LE=  *'  OUTPUT  -  ) 


KEEP  TRACK  OF  THE  CPU  SECONDS  USED. 


TZERO  =  SECOND M 


WRITE ( 6, 600 ) 

600  F  ORMAT ( 1 HO , '  NOPT I  ON ' , / ,  1 H  , 1 X , 7 ( 1H- ) , / , 1 HO , 3X , 

-•  1  DEFLECTION  AND  STRESS  ANALYSIS  ONLY „•',/,  1H0, 


&  3X,  - 

&  3X ,  •' 
S<  3X, 

&  3  X  ,  " 

& 

READ<5, 


2  DEFLECTION  AND  STRESS  ANALYSIS  0F',/,1H  , 

MODIFIED  STRUCTURE  USING  “TAYLOR  SERIES" . ' ,/, 1H0, 

3  DEFLECTION  AND  STRESS  ANALYSIS  0F"',/,1H  , 

MODIFIED  STRUCTURE  USING  SIMPLE  I TERAT ION .-',//,  1H0, 

NOPTION . i') 

*>  NOPTION 


READ  IN  THE.  ORDER  OF  THE  STIFFNESS  MATRIX. 


WRITE (6, 605) 

605  FORMAT (  1  HO,  ■'ORDER  OF  STIFFNESS  MATRIX . i') 
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READ (5,*)  N 
WRITE (6, 610) 

MO  FORMAT!  1H0,  '  YOUNGS  MODULUS . :'> 

READ ( 5»  * )  E 

WRITE (6- 620) 

620  FORMAT (  1  HO,  '  INPUT  LOAD  VECTOR . s  '  ) 

READ ( 5i # )  <F< I ) 7 1*1 ,N> 


COPV  ALL  LOADS  INTO  VECTOR  "OLDF" . 

THE  VECTOR  11  OLDF"  IS  NEEDED  ONLY  WHEN  THE  USER 
SELECTS  THE  3'RD  OPTION. 

DO  5  1  =  1  »N 
5  OLDF  (  I  )  =  F  ( I  ) 

WRITE <6, 630) 

&30  F ORMAT ( 1 HO , y INPUT  MEMBER  SIZES  OF , / , 1 H  , 

&  -UNMODIFIED  STRUCTURE . :') 

READ <  5>  # )  ( X  < I >, 1  =  1 , 10) 

I FLAG  DETERMINES  IF  THE  VECTOR  "X"  OR  "XB"  SHOULD 
BE  USED  FOR  THE  MEMBER  SIZES. 

I FLAG  =  1  - >  "X"  CONTAINS  THE  CURRENT  MEMBER  SIZES. 

I FLAG  =  2  - >  "XB"  CONTAINS  THE  CURRENT  MEMBER  SIZES. 

I  FLAG  =  1 

CALL  GENEK 


1 0  CALL  GENSK ( I FLAG ) 

GENERATE  THE  STRESS  MATRIX. 

CALL  GENS 

CHOLESKX  DECOMPOSITION  OF  MATRIX  "SK". 
CALL  LLT 

FORWARD  SUB'S T I  TUT  I  ON .  . 

CALL  FSUB 

BACKWARD  SUBSTITUTION. 

CALL  BSUB 

PRINT  OUT  SOLUTION  VECTOR 
CALL  F'RTANS 
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COMPUTE  ELEMENT  STRESSES. 


i_: 

C 

CALL  ESTRESS ( I FLAG ) 

C 

C  COMPUTE  CPU  TIME  THAT  HAS  ELAPSED. 

C 

L  f  =  SECOND () 

WRITE (6, 640)  ET-TZERO 

640  FORM  AT  ( 1H0,  'TIME  ELAPSED  FROM  PROGRAM  START  =  ' , F5 . 2 ,  "  CPU  SECS.') 
C 

C  DEPENDING  UPON  THE  VALUE  OF  "NuPTION"  PERFORM  A  DIFFERENT  TASK. 

C 

G 0 T  0  <20,30,4 0 ) ,  NOP  T I ON 

C  NUPTION  =  1  (DEFLECTION  &  STRESS  ANALYSIS  ONLY) 

20  CONTINUE 

c 

C  ASP  THE  USER  IF  SHE/HE  WISHES  TO  CHANGE  THE  STRUCTURE  SOMEWHAT. 

0 

WR I TE ( 6 , 650 > 

650  FORMAT l 1  HO, 'WOULD  YOU  LIKE  TO  CHANGE  THE  MEMBER  SIZES?  (Y/N)  -  > 
READ <5, 500)  ANSWER 
500  FORMAT (Al) 

I F ( ANSWER . NE . " Y ' . AND . ANSWER . NE . "N ' >  THEN 
WR I TE  <  6 , 655 ) 

655  FORMAT (  1  HO,  ■'SORRY.  .PLEASE  REPEAT  YOUR  REPLY.  .  .  '‘Y‘‘  =  YES,  NN"=lM0‘  ) 

GOTO  20 
END  IF 

IF  (ANSWER  .EQ.  "Y" )  THEN 
CALL  CHOS I ZE ( I  FLAG ) 
lFLAG  =  2 
GOTO  10 

ELSE 

STOP  'NOPTION  =  1  '' 

END  IF 

C  NOPTION  =  2  (DEFLECTION  &  STRESS  ANALYSIS  OF  MODIFIED  STRUCTURE 
C  USING  "TAYLOR  SERIES'.) 

iC 

30  CALL  CHGS I ZE ( I FLAG  > 

CALL  OPT 2 
ET  =  SECOND ( ) 

WRITE (6, 640 )  ET-TZERO 
35  WRITE (6, 650) 

READ (5, 500)  ANSWER 

IF ( ANSWER. NE. 'Y'  -AND.  ANSWER. NE. 'N" )  THEN 
WR I TE ( 6 , 655 ) 

GOTO  35 
END  IF 

IF ( ANSWER. EQ. 'N')  THEN 
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STOP  "NOPTION  =  2-' 

ELSE 

I  FLAG  =  2 
GOTO  30 
END  IF 

0 

C  NOPTION  =  3  (DEFLECTION  l<  STRESS  ANALYSIS  OF  MODIFIED  STRUCTURE  USING 
C  SIMPLE  ITERATION. > 

C 

40  CONTINUE 


C  OBTAIN  INFO  CONCERNING  THIS  OPTION. 

C 


CALL  OPTS (NITER) 

C 

C  CHANGE  CERTAIN  MEMBER  SIZES. 

C 

CALL  CHGSI ZE ( I  FLAG ) 

C 

CALL  I  TER ( N I  TER ) 

ET  =  SECOND O 
WRITE( 6,  640)  ET-TZERO 
4S  Wft I TE (6,  650 ) 

READ <5, 500)  ANSWER 

I F  ( ANSWER .  NE .  -  Y  "  .  AND .  ANSWER .  NE .  N '  )  THEN 
WRITE( 6j 655 ) 

GOTO  45 
END  IF 

IF  (ANSWER.  EC!.  "NO  THEN 
STOP  •"NOPTION  -  3"‘ 

ELSE 

I FLAG  =  2 
GOTO  40 
END  IF 
END 

C  ###  #**#***#*##*#**# ft******  ft#**#**#######  ft******#***##**#  ##■&* 

SUBROUTINE  GENEK 
COMMON  /ELEMAT/  EK(S,8,10) 

COMMON  /ORDER/  N> NOPTION 

C 

C  GENERATE  THE  ELEMENT  MATRICES. 

C 

C  INITIALIZE  ALL  ELEMENTS  TO  ZERO. 

C 

DO  5  1=1, N 
DO  5  J=1,N 
DO  5  K=l,10 
5  EK (I, J,K)  =  0.0 
L. 

C  NOW  GENERATE  EACH  NON-ZERO  ELEMENT  FOR  EACH  MATRIX. 

C 

C  MATRIX  Kl. 
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EK (1,1,1) 


=  1.0 

C 

C  MATRIX  K2. 

C 

EK(li  1)2)  =  1  *  0 

EK  (1,3,2)  —  "  1  ■  (j 
EK  (.3,1,2)  =  -1.0 
EK  (3,3,2)  =  1.0 

C 

C  MATRIX  K3. 

f. 

ER (7,7,3)  =  1.0 
C 

C  MATRIX  K4. 

C 

EK.  (5,5,4)  =  1.0 
ER  <  5 ,  7 , 4 )  =  -1.0 
EK (7,5,4)  =  -1.0 
EK (7,7,4;  =  1.0 

C 

C  MATRIX  R5. 

r 

E  h  (  2 ,  2 , 5 )  —  1  •  o 
E  K ( 2 , W , 5 )  —  —1.0 
EK (3,2,5)  =  -1.0 
ER (8,3,5)  =  1.0 

C:  MATRIX  K6 . 

C 

EK (4,4,6)  =  1.0 
EK.,  (  4 ,  6 ,  6/)  =  — 1.0 
EK.  (6,4,6)  =  -1.0 
EK.  (6,6,6)  —  1.0 
C: 

C  MATRIX  K7. 

C 

C  NOTE:  T  =  1 . / (2. *SQRT (2.0) ) 

C 

T  =  0.35355  339 
C 

EK (7,7,7)  =  T 
EK;  (7,8,7)  =  T 
EK  (£i,7,7)  =  T 
EK  <  8 , 8 , 7)  =  T 
C 

0  MATRIX  K8. 

C 

EK ( 1 , 1 , 8 )  =  T 
EK ( 1,2,3)  =  -T 
EK ( 2, 1,8)  =  -T 
EK (2,2,3)  =  T 
C 

C  MATRIX  K9. 
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EK ( 1 , 1 , 9 ) 

»  T 

EKt 1,2,9) 

=  T 

EK ( 2, 1,9) 

=  T 

EK (2,2,9) 

=  T 

EK ( 1 ,5,9) 

=  -T 

s£l 

Ju, 

W 

=  -T 

EK <2,5,9) 

=  -T 

EK (2,6,9) 

=  ~T 

EK (5,5,9) 

=  T 

EK (5,6,9) 

=  T 

EK (6,5,9) 

--  T 

EK (6,6,9) 

=  T 

EK<5, 1,9) 

=  -j 

EKt 5, 2,9) 

=  -T 

EKt 6, 1,9) 

=  -T 

EK (6, 2,9) 

=  -T 

MATRIX  K1Q. 

E 1'  (  3  i  3 ,  1 0 ) 

=  T 

EKt 3, 4, 10) 

=  -T 

EK  (  4 ,  j ,  10 ) 

=  -T 

EK  (  4 , 4 ,  1 0 ) 

=  T 

EK (3,7, 10) 

=  -T 

EKt 3, 8, 10) 

=  T 

EK ( 4 , 7 , 1 0 ) 

=  T 

EK (4,8,10) 

=  -T 

EK  (  7 , 7 ,  1 0 ) 

=  T 

EK (7,8, 1 0 ) 

=  -T 

EK (8,7,10) 

=  -T 

EK (8,8, 10) 

=  T 

EK (7,3, 10) 

=  -T 

EK (7,4, 10) 

=  T 

EK (8,3, 10) 

=  T 

EK (8,4, 10) 

=  -T 

RETURN 

END 

c  ****«*******#*#******#***********«********#**************.».****» 

SUBROUT  I NE  GENSK ( I FLAG ) 

COMMON  /ELEMAT/  EK<8,8, 10) 

COMMON  /ORDER/  N»NOPTION 
COMMON  / ST I FFM /  SK ( 8 , 8 ) , R ( 8 ) 

COMMON  /YMLW  E, F ( 8 ) , OLDF ( 8 ) , X ( 10 ) , XB ( 10) 

DIMENSION  XSIZE(IO) 

C; 

C  FORM  STIFFNESS  MATRIX; 

C 

C  INITIALIZE  STIFFNESS  MATRIX  TO  ZERO. 

f 

DO  10  1=1 iN 
DO  10  J=1,N 
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1 0  SF  (  I  ,  J ) 


0.  0 


t;  "XSI2E"  SHOULD  CONTAIN  THE  MODIFIED  MEMBER  SIZES, 

C  IF  I FLAG  =  1  THEN  THE  PROGRAM  IS  PERFORMING  THE  INITIAL  ANALYSIS 
C  ELSE  I FLAG  =  2  WHICH  INDICATES  THAT  THE  USER  HAS  MODIFIED  THE 
C  MEMBER  SIZES. 

C 

IF  ( 1FLAG  .EQ.  1)  THEN 
DO  15  1=1,10 
15  XSIZEI I )  =  X  < I ) 

EL  SE 

DO  18  1=1,10 
18  XSIZE(I)  =  X  B ( I ) 

END  IF 

C 

C  COMPUTE  STIFFNESS  MATRIX. 

C 

DU  20  1=1,10 
DO  20  J=  1 , N 
DO  20  K=1,N 

20  Sk(J,K)  =  SK(J,K)  +  XSIZE( I )*EK< J,K,  I  > 

C 

CONST  =  E/360. 

C 

DO  30  I=i,N 
DO  30  U= 1 , N 

30  SK  < 1 , J )  =  CONST *SK ( I , J ) 

C 

RETURN 

END 

SUBROUTINE  GENS 
COMMON  /ORDER/  N,NOPTION 

COMMON  / YMLV/  E, F ( 8 ) , OLDF < 8 ) , X ( 10 ) , XB ( 10) 

COMMON  /STRESS/  S ( 10, 3 ) , ST ( 10 > 

C 

C  GENERATE  STRESS  MATRIX. 


DO  1 0 

1=1 , 10 

DO  1 0 

Z 

i 

II 

10  S ( I , J  > 

=  0.0 

GENERATE 

NON-ZERO  ELEMENTS 

S  <  1  >  1  ) 

=  E/360.0 

3(2,1) 

=  -S (1,1) 

St  2, 3) 

=  -S (2, 1 ) 

S  (  3 ,  7  ) 

=  S( 1 ,  1 ) 

S(4, 5) 

=  SCI. 1 > 

S  ( 4 , 7 ) 

=  -S (1,1) 

'“'(Si  ) 

=  -S( 1,1) 

3(5,8) 

=  S(  1,  1 ) 

8(6,4) 

=  — S ( 1,1) 
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3(6,6)  =  S<  1  ,  1  ) 

S  <  7 »  7  >  =  0.5*S(1,1) 

S(  7,8)  =  S<7,7) 

8(8,1)  =  S(7,7> 

3(8,2)  =  -3(7,7) 

1 )  =  -8(7,7) 
ft <8,2)  =  -8(7,7) 

5 '  =  8(7,7) 

SC-C6)  =  8(7,7) 

S< 10,8)  =  8(7,7) 

M 10-4)  -  -8(7,7) 

St  10, 7)  =  -8(7,7) 

8 < 10,8)  -  8(7,7) 

C 

RFTURN 

END 

SUBROUTINE  ESTRE33 ( I FLAG ) 

COMMON  /ORDER/  N,NOPTIGN 
COMMON  /STIFFM/  SK(S,8),R(8) 

COMMON  /STRESS/  8 ( 1 0, 8) , ST ( 10 ) 

COMMON  /YMLV/  E, F(8) , OLDF (8 ) , X ( 10 ) , XB ( 10 ) 

DIMENSION  XSIZE(IO) 

i  COMPUTE  ELEMENT  STRESSES. 

DO  10  1=1,10 
ST ( I )  =  0.0 

DO  10  J= 1 , N 

ST ( I >  =  ST ( I )  +  S(I,J)#R(J) 

XF  ( I FLAG  .EQ.  1)  THEN 
DO  15  1=1,10 
XSIZE( I )  =  X ( I ) 

ELSE 

DO  18  1=1,10 
XSI ZE ( I )  =  XB (  I  ) 

END  IF 

WRITE (6, 600) 

FORMAT! 1 HO, /, 1 HO, 7X, "ELEMENT  STRESSES  AND  SIZES",/, 1H  ,7X,26(1H-)) 
DO  20  1=1,10 

WRITE <6, 610)  I ,ST( I ) » I, XSIZEf I ) 

FORMAT (1H  , "ST( ", 12, "  >  =  " , F10. 3, 2X , " X ( " , 12, '  )  «  ",F10.3) 
CONTINUE 

C  COMPUTE  WEIGHT. 

C 

WEIGHT  =  0.0 
WEIGHTD  =  0.0 

C 

DO  30  1=1,6 


i. 


10 


15 


t  E 


600 


610 

20 


l: 
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c 

c 

c 

c 

c 


1 , 

c 

r. 

I 


SO  WEIGHT  =  WEIGHT  +  XSIZE(I) 

DO  40  1=7,10 

40  WEIGHTS  =  WEIGHTD  +  XSIZE(I) 

WEIGHT  =0,1  *  (360.0  *  WEIGHT  +  360.0  *  SORT (2.0)  *  WEIGHTD) 

WRITE (6,620)  WEIGHT 
620  FORMAT ( 1 HO , ' WE I GHT  =  ' ,  F 1 0 .  3 ) 

RETURN 

END 

*  ■»*+******#*■»***##*#***#*#***#****#********#**********»******** 
SUBROUTINE  LLT 
COMMON  /ORDER/  N.NOPTION 

COMMON  /YMLV/  E, F ( S) , OLDF ( S ) , X ( 10 ) , XB ( 10 ) 
i  GMMuN  / ST  I FFM /  SK ( 8 , S ) , R (8 ) 

ROUTINE  "LLT"  PERFORMS  THE  CHOLESKI  DECOMPOSITION  OF 
THE  COEFFICIENT  MATRIX: 


C  SK  =  (L)  *  (LT) 

C  WHERE , 

i  L  =  LOWER  TRIANGULAR  MATRIX 

C  (STORED  IN  LOWER  TRIANGLE  OF  SK) 

i 

C  LT  =  TRANSPOSE  OF  L 

C 

DO  4  I = 1 , N 
DO  4  J=I,N 
S  =  SK  ( I ,  J  ) 

IF  cl  . EQ.  1)  GOTO  2 
11  =  1-1 
C 

DO  1  K= 1,11 

1  S  -  S  -  SK( I,K)*SK< J,K) 

2  IF  (J  . NE.  I)  GOTO  3 

IF  (S  .LE.  0.0)  STOP  'NOT  POSITIVE  DEFINITE  -  STOP  IN  LLT" 

T  =  SORT ( ABS ( S ) ) 

SK < 1,1)  =  T 
GOTO  4 

3  SM  J,  I  >  =  S/T 

4  CONTINUE 
C 

RETURN 

END 

C  ###  ■» -tHHJ- -fr -ft  tf -tt -fr -IHt  # -ft 

SUBROUTINE  FSUB 

COMMON  /ORDER/  N,NOPTION 

COMMON  /YMLV/  E, F(S) , GLDF(S) , X< 10) , XB( 10) 

COMMON  /ST I FFM/  SK(8,8),R(S) 

C 

C  ROUTINE  "FSUB"  PERFORMS  THE  FORWARD  SUBSTITUTION  TO  SOLVE 


# 


* 
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C  L*Y  a  B  WHERE  L  IS  STORED  IN  THE  LOWER  TRIANGLE  OF  A  AND 
C  Y  TS  STORED  IN  THE  X  VECTOR. 

C 

DO  10  I  =  1 , N 
S=F< I ) 

ini  =  i-i 

IF  (I  .EG.  1)  GOTO  10 

C 

DO  20  d=l,IMl 
20  S  =  S  -  SK<X.df*R(j) 

( 

10  R. I )  =  S/SK ( I » I ) 


return 

END 

SUBROUTINE  BSUB 
COMMON  /ORDER/  N,NGPTION 

COMMON  / YMLV/  E»F<8) iOLDF(S) , X < 10 ) >  XB< IQ) 
COMMON  /STIFFM/  SK<8,8),R(8> 


C  ROU'l  INE  "BSUB"  PERFORMS  THE  BACKWARD  SUBSTITUTION  TO  SOLVE 
C  L.r*R  =  Y  WHERE  L  IS  STORED  IN  THE  LOWER  TRIANGLE  OF  8K  AND 
V  IS  INITIALLY  STORED  IN  X. 

C 


C 

C 

c 

c. 


ON  RETURN  R  CONTAINS  THE  SOLUTION  OF  SK*R  =  F 


WHERE i 


SK  =  THE  ORIGINAL  SYMMETRIC  COEFFICIENT 


MATRIX. 


DO  3  I  =  1 ,  N 
Nil  -  N-I+l 
II  =  1-1 
S-  R  <  N II  ) 

IF  (I  .EQ.  1)  GOTO  2 
DO  1  J-1,11 

1  S  =  S  -  SK(NI 1+J, Nil >*R<NI 1+J> 

2  R < N 1 1 )=S/SK(NI 1 .Nil) 

3  CONTINUE 
RE  TURN 
END 

SUBROUTINE  PR TANS 
COMMON  /ORDER/  N.NOPTION 
COMMON  /STIFFM/  SKC8,8),R<8> 


C  ROUTINE  "PRTANS"  PRINTS  OUT  THE  SOLUTION  VECTOR  " X" 

C 

DO  10  I -1 » N 
10  WRITE ( 6, 600 )  I,R(I) 

600  FORMAT  (  1  HO, ’"RC 12, ''  )  =  ",F10.4) 
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c 

RETURN 

END 

C  HUHHHHHUHHHHUHHHHHHt#**#####****##**#****#*^^***#****##**#*#* 
SUBROUT  I NE  OPTS ( N I  TER ) 

COMMON  /EPS/  TOL 

COMMON  / YMLV/  E , F ( S ) , OLDF ( S ) , X ( 10 ) , XB (10 ) 

t" 

WRITE (6, 600) 

600  FORMAT  (  1  HO, ''WHAT  NUMBER  OF  ITERATIONS  WOULD  YOU  LIKE . :■) 

READ <“5,  *)  NITER 
C 

WR I TE  <  6, 610 ) 

610  FORMAT U HO,  WHAT  TOLERANCE  DO  YOU  REQUIRE  OF  THE', 

+  '* SOLLIT I  ON  VECTOR . .  '  ) 

REAB<5,*>  TOL 

t 

RETURN 

END 

SUBROUTINE  SETUP 

COMMON  /EL'EMAT  /  £K  (0,8,10) 

COMMON  /ORDER/  N,NOF'TION 
COMMON  /ST  I FFM /  SK (8,8) , R ( 8 ) 

I..OMMON  / YML.V/  E,  F(3) ,  OLDF (8)  ,  X  (  10) ,  XB(  10) 

D I  MENS  I ON  DEL TAK ( 8 , 8 ) , TEMP ( 8 ) 

O 

c.  COMPUTE  “DELTA  K"  MATRIX. 

CALL  GENDK (DELTAK ) 

i1 , 

C  COMPUTE  THE  RIGHT  HAND  SIDE  OF  EQ.(17) 

i.  FIRST  MULT.  C DELTA  K3  #  (R) 

i.: 

DO  30  I = 1 , N 
TEMP(I)  *  0.0 
DO  30  J=1,N 

30  TEMPI  I)  =  TEMP (I )  +  DELT  AK( I , U>  *  R(J) 

C 

DO  40  1=1, N 

40  F  < I )  ■  OLDF  < I )  -  TEMP ( I ) 

C 

RETURN 

END 

C  tt#*#***#*#**###*##*##**#**##**##***###***#.**##**#* *  •#■•#■■***** #**•*« 

SUBROUTINE  ITER ( NI TER ) 

COMMON  /EPS/  TOL 

COMMON  /ST  I  FFM/  SMS,  8),  R  (8) 

DIMENSION  OLDR ( 8 ) 

C 

C  this  ROUTINE  PERFORMS  THE  DEFLECTION  it  STRESS  ANALYSIS  VIA 
0  SIMPLE  ITERATION. 


-4 
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DO  10  1=1 i NITER 


SET  UP  THE  RHS  OF  EQ.(17> 
i.  ALL  SETUP 

SAVE  THE  OLD  SOLUTION  VECTOR  'H'  IN  •'OLDR'' . 

DC, i  20  U=  1 , 8 
20  OLDft(J)  =  R ( J ) 

SOLVE  THE  SYSTEM  FOR  A  NEW  "RHS". 

FIRST  CALL  THE  FORWARD  SUBSTITUTION  PART _ 

CALL  FSUB 

...NOW  USE  THE  BACKWARD  SUBSTITUTION. 

CALL  BSUB 

TEST  1U  DETERMINE  IF  MAGNITUDE  OF  "OLD"  SOLUTION  VECTOR 
IS  WITHIN  A  GIVEN  TOLERANCE  OF  THE  "NEW"  SOLUTION  VECTOR. 

DO  30  J- 1 , S 

D I FF  =  ABS ( OLDR ( J)  -  R < J > ) 

IF  <  DIFF  .LT.  TOL)  GOTO  30 
GOTO  10 
30  CONTINUE 

WRITE (6, 620)  I 

620  FORMAT  <  1  HO »  'THE  ANALYSIS  REQUIRED  ITERATIONS.') 

GOTO  40 

10  CONTINUE 

WRITE (6, 620)  NITER 

PRINT  OUT  THE  ANSWERS... 

40  CALL  PRTANS 

PRINT  OUT  THE  ELEMENT  STRESSES  AND  SIZES. 

CALL  ESTRESS(2 ) 

RETURN 

END 

tt  ft  ***■»■«  *******  *  *  #  #  *#**#**###*#*  -ft  *  #  *  *  ft  *  #■  #  #  ft 

SUBROUTINE  CHGSI ZE ( IFLAG ) 

COMMON  / YMLV/  E, F ( 8> , 0LDF(8) , X < 10 ) , XB<10 ) 

THIS  ROUTINE  ALLOWS  THE  USER  TO  MODIFY  THE  STRUCTURE. 
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5  WRITE ( 6>  600 ) 

600  FORMAT  (  1  HO, 'HOW  MANY  MEMBER  SIZES  WOULD  YOU  LIKE  TO  CHANGE  :  ■"  > 
READ  (  5 .  *  )  NCHGES 
IF  (NCHGES  .GT.  10)  GOTO  5 
IF  ( NCHGES  . LE.  0)  RETURN 
L 

IF ( IFLAG. EQ. 1 )  THEN 

C 

i  COF'V  EXISTING  MEMBER  SIZES  TO  MODIFIED  MEMBER. 

C 

no  s  1=1,10 

S  XB( I j  =  X ( I ) 

END  IF 

OU  to  1=1, NCHGES 
15  WRITE (6, 610) 

6  1 0  K  ORMAT  <  1  HO ,  '  ELEMENT  NUMBER . :  '  > 

READ  (  5 ,  *>  NIJMELEM 

r 

r  CHECH  TO  DETERMINE  IF  MEMBER  SIZE  VALID. 

IF  (NUMELEM  . LT.  0  .OR.  NUMELEM  .GT.  10)  THEN 
WR I TE ( 6 , 900 )  NUMELEM 

900  FORMAT <1  HO,  ^  ***  NO  SUCH  ELEMENT  NUMBER  — >  ',13,'  ***■') 

GOTO  15 
END  IF 

WRITE (6, 620) 

620  FORMAT (  1H  ,  -  MEMBER  SIZE . :  '  ) 

READ(5, *)  SIZEMEM 
10  XB (NUMELEM )  =  SIZEMEM 

RETURN 

END 

I.  ft-tf-tt-it  ft  1*^ 

SI  IBROUT I NE  GENDK  (  DELTAK ) 

COMMON  /ELEMAT/  EK(8,8, 10) 

COMMO N  /ORDER /  N , N OP T I  ON 

COMMON  / YMLV /  E , F ( S ) , OLDF ( S ) , X ( 1 0 ) , X  B ( 1 0 ) 

D I MENS  I ON  DELTAK (8,8) 

C 

C  INITIALIZE  THE  MATRIX  "DELTAK"  TO  ZERO. 

C 

DO  5  I = 1 , N 

DO  5  J= 1 , N 

5  DELTAK ( I , J)  =  0.0 

C 

C  COMPUTE  DELTA  K  (EQ.  (18)) 

C 

DO  10  1=1,10 
DIFFX  =  XB ( I )  -  X < I > 

C 

DO  20  J=1,N 

DO  20  K=1 , N 
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D I FF X  *  EK ( J 1K1 I )  +  DELTAK ( J , K > 


20  DELTAK ( J, K)  = 

10  CONTINUE 

C 

DO  23  IE1 i N 

DO  23  'J=1>N 
23  DELTAK ( I , J )  =  <  E/360.  )  *  DELTAK  <  I »  J ) 

C 

RETURN 

END 

C  »**■»*■****##***###»#****#*********•»****#**#*#***#***#***##*■» 

SUBROUTINE  OPT 2 
i COMMON /ELEMAT /  EK (8>8. 10  > 

COMMON/ YMLV/  E ,  F  <  8  >  , OLDF ( S ) ,  X  ( 1 0 > , XB ( 1 0 ) 
i  OMMON /ST I FFM /  SK (S,8>,R(8) 

D I MENS I ON  RQLD ( S ) , DELR ( 8  > 

L.  STORE  OLD  R  VECTOR 

i; 

DO  10  1=1,8 
10  ROl  D< I)«R< I ) 

t  INITIAL  I ZE  DELR=0. 
r 

DO  20  1=1,8 
20  DELR (I)  =0.0 
C 

t ;  L  OOP  F OR  SUMMAT I ON ,  EQ .  (15) 

C 

DO  200  J=  1 , 1 0 

C 

C  FORM  R I  OH T HAND  SIDE  OF  EQ.  (16) 

C 

DO  120  1=1,8 
F ( I >=0.0 
DO  110  L=1 , 8 

110  F ( I ) =F C I ) +EK ( I , L , J ) *ROLD ( L ) 

120  F ( I >=-F( I )*E/360. 

C 

C  SOLVE  EQ.  (16)  FOR  R 

C 

CALL  FSUB 
CALL  BSUB 
C 

C  SUM  DELR 

C 

DO  1 30  1=1,8 

130  DELR < I ) =DELR (I )  +  <XB(J)-X(J) )*R( I ) 

200  CONTINUE 

C 

C  SOLVE  EQ.  (15) 

C 

DO  210  1=1,8 
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210  ft ( I ) =R0LD ( I ) +DELR ( I ) 


r  PRINT  OUT  SOLUTION  VECTOR 

C 

CALL  PRTANS 

C 

C  COMPUTE  ELEMENT  STRESSES 

C 

CALL  ESTRESS (2) 

>■  REF’LACE  NEW  R  VECTOR  WITH  ORIGINAL 
C  (RETURN  TO  ORIGINAL  STRUCTURE) 

i 

DO  220  1=1,8 
220  R ( I ) =ROLD ( I ) 

RETURN 

END 


4 


l 
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